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CHAPTER 1 



Introduction 

1.1 Motivation 

Lepton-hadron deep inelastic scattering (DIS) is a very important tool to unravel the 
structure of hadrons. To the lowest order, this scattering takes place via the exchange 
of a virtual photon. When a very low mass virtual photon (Q 2 < lGev 2 ) scatters off a 
proton, the photon 'sees' only the total charge and magnetic moment of the proton and the 
proton appears to be a pointlike object. However, a higher virtual mass photon resolves 
the fine structure of the proton charge distributions and sees its elementary constituents. 
The cross section of such scattering can be expressed in terms of several functions which 
depend on two variables, Q 2 (momentum transfer square in the process) and v (energy of 
the virtual photon). These functions, to start with, are unknown and they depend on the 
hadron structure. They are called structure functions. Early SLAC experimental results 
showed that the structure functions show a phenomenon called scaling, which means that 
they depend only on the ratio of Q 2 and v, and do not depend on them separately. In 
the pre-QCD era, Bjorken used current algebra techniques to explain such behavior of 
structure functions. Later, Feynman proposed parton model idea according to which, 
the proton consists of pointlike, massless, noninteracting particles called partons and 
the scattering takes place incoherently from the partons. The parton model successfully 
explained scaling of the structure functions. However, substantial scaling violations were 
observed in later SLAC experiments and the structure functions were found to evolve with 
Q 2 . This violation of scaling occurs because of the interactions between the constituents of 
the proton, which are now known as the quarks and gluons. Their interaction is governed 
by QCD (quantum chromodynamics). 

The Q 2 dependence of the structure functions can be logarithmic or there may be 
power dependence. The contributions to the structure functions which are proportional 
to some power of ^ are commonly called 'higher twist' contributions. The most widely 
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used formal approach to DIS structure functions is based on Wilson's operator product 
expansion (OPE). There exists a more formal definition of twist in terms of OPE, which 
we shall discuss shortly. DIS is a high energy process, and Q 2 is very large. This means 
that higher and higher twist contributions are suppressed compared to the leading twist 
contributions. However, recent experiments indicate that the higher twist effects play 
a very important role in the SLAC kinematical range. Understanding the higher twist 
effects requires non-trivial non-perturbative information of the structure of the hadron. So 
far, a clear and physical picture of the higher twist contribution is lacking. To understand 
these effects, there is an urgent need to develop a non-perturbative technique to analyze 
DIS. For this, one needs theoretical tools which are based on physical intuitions and at 
the same time employs well-defined field theoretic calculational procedure. 

The usual field theoretical formulation is Lagrangian formulation, in which one starts 
from the local Lagrangian, gets the equations of motion and then quantizes the system by 
imposing suitable commutation or anti-commutation relations between the fields and their 
conjugate momenta on a space-like hypersurface. It is a manifestly covariant formalism. 
However, for interacting theory, the equations of motion involve non-trivial interacting 
fields and can only be solved approximately, using perturbative techniques. Perturbative 
methods cannot be applied to low energy QCD bound state problems since the coupling 
constant takes a large value in this energy range. Another point is that, in the covariant 
Lagrangian formulation, the intuitive picture of quantum mechanics is lost. There exists 
another formulation known as Hamiltonian formulation of quantum field theory in which 
one starts with the Hamiltonian, obtains the Hamiltonian equations of motion and sets 
the quantization conditions on a space-like hypersurface. The Hamiltonian gives the 
time evolution of the system. Though it is physically more intuitive, there are some 
disadvantages of Hamiltonian formulation. It is manifestly non-covariant since one has 
to choose a time axis and perturbation theory becomes more complicated. The latter is 
due to the fact that here one has to do old-fashioned time ordered perturbation theory 
and instead of Feynman diagrams, one has to calculate a larger number of time-ordered 
diagrams at each order. It has been shown that, light-front quantization (where the surface 
of quantization is tangential to the light cone) makes Hamiltonian formulation simpler. 
This is due to the fact that, the diagrams which produce particles from the vacuum are 
absent here (this point will be clarified in chapter 3). So one encounters lesser number 
of time ordered diagrams. In what follows, we have used Hamiltonian light-front QCD 

2 



framework (Hamiltonian formulation of QCD where the quantization surface is tangential 
to the light-cone) to analyze DIS structure functions. 

In order to get an intuitive picture of DIS, it is necessary to keep close contact with 
Feynman's parton model. As mentioned before, partons were originally introduced as 
collinear, massless, non-interacting constituents of the proton. It is necessary to generalize 
the parton ideas to introduce field theoretic partons which are massive, non-collinear and 
their interaction is governed by QCD. This goal is achieved in light-front Hamiltonian 
QCD analysis. This is essentially a non-perturbative method which allows us to explore 
the non-perturbative higher twist effects explicitly. In this approach, instead of calculating 
the moments of the structure functions as one does in the conventional OPE formalism, 
one directly calculates the structure functions themselves which involve hadron matrix 
elements of light-front bilocal operators. The structure functions can be calculated once 
the light-front wave functions are known. However, perturbative calculations are also 
possible, and both perturbative and non-perturbative aspects can be investigated within 
the same framework. 

In what follows, we have used light-front Hamiltonian QCD framework to study the 
twist four part of the longitudinal structure function F L in unpolarized DIS and the 
transverse polarized structure function g T in polarized DIS, which is a twist three con- 
tribution. Both of these structure functions involve non-trivial operators which brings in 
non-perturbative dynamics. It is well known that the evolution of F L and g T with the 
scale, Q 2 , are highly complicated. We have shown that the reason of this complexity is the 
fact that F L is related to the light-front QCD Hamiltonian density and qt is related to the 
transverse spin operator, both of which are complicated dynamical operators (they change 
the surface of quantization, see chapter 3). As a result, both these structure functions 
acquire off-diagonal contributions from the matrix elements, in contrast to the other two 
electroproduction structure functions, F 2 and g±, which involve diagonal matrix elements 
only. Another interesting point is, the operators involved acquire divergences in pertur- 
bation theory. The experimentally measured structure functions, of course, are finite 
quantities. Therefore, one must know how to remove these divergences (renormalization). 
We have addressed the issue of renormalization of these composite operators which, in 
light-front theory is highly non-trivial. As we mentioned before, a very interesting point 
is that, both perturbative and non-perturbative issues can be address in our approach. In 
the non-perturbative context, we have analyzed twist-four longitudinal structure function 
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for a meson in 1 + 1 dimensional QCD and for a positronium in 3 + 1 dimensional QED 
in weak coupling limit. 

Earlier, it has been shown that the integral of the structure function F 2 is related to 
the light-front longitudinal momentum density and the integral of the polarized structure 
function g\ is related to the helicity. In this work, we have shown that the integral of Fl is 
directly related to the light-front Hamiltonian density. Also the integral of the transverse 
polarized structure function is related to the transverse spin operator in light-front 
QCD. Thus our approach gives a unified picture of DIS structure functions and connects 
them to the light-front Poincare generators. We emphasize that such connections arise 
because our approach is different from the usual one and we get many new and interesting 
results. 

1.2 Plan of the Thesis 

The usual approach to analyze DIS structure functions is based on OPE. This mainly 
concentrates on the evolution of the structure functions in perturbative QCD. In these 
methods, the moments of the structure functions automatically arise. We briefly discuss 
the OPE method to analyze DIS structure functions in chapter 2. In this chapter, we also 
explain the term 'twist' in the context of OPE. It is also important to understand the 
origin and nature of the higher twist contributions. We discuss how these higher twist 
effects are addressed in the covariant formalism. For the sake of simplicity, we briefly 
outline the procedure and omit the mathematical details. The effect of the target mass 
gives a higher twist contribution, which, in the free theory assumption produces £ scaling. 
The dynamical higher twist effects are highly complicated. We point out that in the 
QCD improved parton model approach, these contributions are expressed as the Fourier 
transform of the target matrix elements of various bilocal operators. However, these 
matrix elements cannot be further investigated because of the non-perturbative hadronic 
state. 

We follow an alternative approach based on light-front Hamiltonian QCD. The light- 
front formulation of quantum field theory has many strikingly different features compared 
to the usual equal time formulation. Some of these features makes it more suitable for 
the analysis of DIS process. It is important to remember that the parton model was 
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formulated in the infinite momentum frame, which, in some sense is equivalent to light- 
front formulation. In chapter 3, we discuss the basic features of light-front field theory 
and light-front QCD. 

In chapter 4, we outline our approach. Bjorken- Johnson-Low (high energy) expansion 
of the Compton scattering amplitude gives the expressions of the structure functions as 
the Fourier transform of the matrix elements of light-front bilocal operators. We work in 
the light-front gauge. In this method, one gets the structure functions themselves and not 
their moments. By expanding the target state in Fock space and using truncation, the 
structure functions are expressed in terms of light cone wave functions. Splitting functions 
are obtained in perturbation theory simply by replacing the target by a dressed parton. 
In this chapter, we review the works done so far in this approach. We show that this 
method not only reproduces the well known results of the covariant formalism, like the 
Q 2 evolution of the twist two structure function F 2 , but also addresses some important 
issues, like the effect of quark mass in and the helicity sum rule of the proton, in a 
more transparent way. 

The actual work done in this research program has been discussed in detail in chapters 
5, 6, 7 and 8. These works are concerned with the higher twist effects F L and qt- All of 
these results are new and interesting, both from theoretical and phenomenological point 
of view. Since our approach is different from the usual one, such new outcomes are not 
unexpected. Also, our approach gives better insights into problems that are not easy to 
understand in the commonly used procedure. 

In chapter 5, we discuss a new sum rule that we have obtained which connects the twist 
four part of Fl to the fermionic part of light-front QCD Hamiltonian density. Introducing 
a gluonic part of Fl, we show that the sum rule relates the total Fl to the target mass. 
To our knowledge, this is the first sum rule for a higher twist structure function. In 
this chapter, we also perform one loop renormalization of Fl in light-front Hamiltonian 
perturbation theory. Replacing the target by a dressed quark, we show that the twist four 
part of Fl is directly related to the quark mass shift in perturbation theory. 

In chapter 6, we continue our investigation of Fl- In the non-perturbative context, 
we verify that the sum rule is obeyed for a meson in 1 + 1 dimensional QCD because 
of t' Hooft equation. An analysis for a bound state in QCD in 3 + 1 dimension is very 
complicated and there one has to use the recently developed similarity renormalization 
technique. In order to understand the calculational procedure, we consider a simpler case, 
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namely a positronium-like bound state in 3 + 1 dimensional QED in the weak coupling 
limit. We explicitly verify the sum rule in this case. 

In chapters 7 and 8, we investigate another most challenging problem, namely that 
of the spin operator for a composite system in an arbitrary reference frame in relativistic 
quantum field theory. We show that, this issue is very important in the context of polarized 
scattering, since the transverse polarized structure function g T is related to a part of the 
transverse spin operator. It is also known that g± is related to the helicity operator. In 
chapter 7, we discuss the problems associated with the definitions of the relativistic spin 
operators in equal-time theory and contrast with the light-front case. The light-front 
theory is more suitable for this purpose because boost is kinematical. We discuss the 
light-front spin operators for the massive and massless cases separately. The transverse 
spin operator contains the transverse rotation operator which is interaction dependent. 
We construct the transverse rotation operator in light-front QCD. In this chapter, we 
also show that though the light-front transverse spin operators cannot be separated into 
orbital and spin parts, there exists a decomposition of it which is physically interesting. 

In chapter 8, We explore the physical significances of this decomposition and compare 
and contrast it with the decomposition of the helicity operator into orbital and spin parts. 
We also perform one loop renormalization of the full transverse spin operator. Our results 
show that only one counterterm is needed and that has the same form as the linear mass 
counterterm in light-front Hamiltonian perturbation theory. We explicitly show that in 
an arbitrary reference frame, all the center of mass momentum dependence get canceled 
in the matrix element of the transverse spin operator. We emphasize that in the equal 
time case, this cancellation is almost impossible to prove because of the dynamical boost 
operators. 

Summary and conclusions are given in chapter 9. References are given at the end 
of each chapter. In order to clarify the notations and also the intermediate steps of 
calculations, we have provided several appendices. 
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CHAPTER 2 



Operator Product Expansion 
Approach to DIS 

The standard approach to DIS structure functions is based on Wilson's Operator 
Product Expansion (OPE). In this chapter, we introduce the idea of OPE and also the 
definition of twist of an operator. Then we describe briefly the OPE method of analyzing 
DIS. 

2.1 'Twist' of an Operator 

Products of fields at the same space-time point are called composite operators. The 
composite operator is not well-defined. The easiest way to see this is to calculate the 
vacuum expectation value of the propagator of free scalar field, A(x — y). In coordinate 
space, one clearly sees that A(x — y) diverges as x — > y. For free field theory, one then sub- 
tracts out the singularity by taking the normal ordered product. However, for interacting 
fields, this cannot be generalized in a straightforward manner. Wilson introduced the idea 
of operator product expansion (OPE) in which the singularities of operator products are 
expressed as a sum of non-singular operators with singular c-number coefficients. There 
are two types of OPE |J: 

2.1.1 Short Distance Expansion 

The product of two operators at short distance (x — yY — ^ is given by, 

A(x)B( y ) « J2 c i( x - v)Oi{\{x + y)); (x - yf - (2.1) 

where A, B and Oj are local operators. are singular c number functions which are called 
Wilson coefficients. The short distance behavior of the Wilson coefficients is expected to 
be given by naive dimensional counting. The most dominant term in the expansion at 
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x — > will be the term for which A ~ B is most singular, where di, (La and cLb are 
the mass dimensions of the operators. Higher is the dimension of the operator O;, less 
singular is the Wilson coefficient. 

2.1.2 Light Cone Expansion 

The product of two operators at short light cone distance x 2 — > is given by, 

^(f M-f) « Ec<(*Mf, -|). * 2 « °- ( 2 - 2 ) 

Oi(x,y) are regular bilocal operators which can be expanded in a Taylor series, 

°t\>-\) ,,,(")■ (2-3) 

So the OPE on the light cone becomes, in terms of local operators, 

A&B(~) ^^(^...^^(O) (2.4) 

where j is the maximum spin of the operators Oj£ w (0)- The light cone behavior of 
the Wilson coefficients is given by naive dimensional counting, and the most dominant 
singularity is for the term which has the lowest value of (d^ — j), i.e. the dimension of 
Ofil Mj . (0) minus the spin of 0^1^(0). The twist of an operator is defined as 

r = dimension — spin. (2.5) 

2.1.3 Spin of an Operator 

The representations of the homogeneous Lorentz group are labeled by the values of 
two positive integers and/or half integers, A, B. An operator that transforms like the 
(A, B) representation has components that rotate like objects of spin j, with 

j = A + B,A + B - 1, .... \A-B\. (2.6) 

A traceless symmetric tensor of rank 2A transforms like the (A, A) representation and has 
components of spin 2A,2A-1,...Q 0. In the definition of twist, one takes into account 
the maximum spin. 

Deep inelastic scattering is a light cone dominated process, which means that it gets 
major contributions near the light cone and the light cone expansion is applied here. 

In the next section, we shall describe briefly the conventional approach based on OPE 
to deep inelastic scattering. 
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2.2 A Brief Overview of Deep Inelastic Scattering in the Con- 
ventional Approach 



In this section, we describe the approach based on OPE to polarized and unpolarized 
deep inelastic scattering. Because of the complexities of the expressions, we shall avoid 
the explicit details but only outline the procedure. 

We begin with a brief review of the basic ingredients of lepton-nucleon deep inelastic 
scattering (DIS): 

e(k) + h(P) — > e(fc') + X(P + q) , (2.7) 

where we have specified the four momenta of the particles explicitly and q = k — k' is 
the momentum transfer in the process through the virtual photon. The inclusive cross 
section for the above scattering process is given by 

(2 - 8) 

Here E (£") is the energy of the incoming (outgoing) lepton, is the leptonic tensor, 

1 

s' 

= 2{k'^k u + k'^k^) — 2g flI/ k ■ k' — 2ie flI/pa q p s (T . (2.9) 

W pv is the hadronic tensor which contains all the hadronic dynamics involved in DIS 

process, 

W " U = ^J d ^ el9< (PS\[J"(0,J u m\PS), (2.10) 
where P and S are the target four-momentum and polarization vector respectively (P 2 = 
M 2 ,S 2 = —M 2 ,S ■ P = 0), q is the virtual-photon four momentum, and J^(x) = 
J2 a Za.ipct(. x )l^''Pa{.x) the electromagnetic current with quark field ip a (x) carrying the flavor 
index a and the charge e a . 

The above hadronic tensor can be decomposed into independent Lorentz invariant 
scalar functions: 

= (-g^ + ^-)w 1 (x,Q 2 ) + ( y P^-^)( y P l/ -^)w 2 (x,Q 2 ) 



S a W 3 (x,Q 2 )+P a S-qW A (x,Q 2 ) 

»n v , A . M 2 



= (^-^)(^l(o:,Q 2 )-^F 2 (x,Q 2 )) + [W-^(PV + PY 

+r ^_\ F^x^l _ te ^^\ S(TLgi ( X} Q 2 ) + SaTgT (x, Q 2 )} . (2.11) 



q 2 



v v 



The dimensionless functions 

Fl(x,Q 2 ) 



-W 1 + [M 2 - 



(p. q y 



Wo 



(2.12) 
(2.13) 



F 2 (x,Q 2 ) = vW 2 (x,Q 2 ) 

are the unpolarized structure functions measured from the unpolarized target and 

gi (x,Q 2 ) = u[w 3 (x,Q 2 )+uW 4 (x,Q 2 )} (2.14) 
g T (x,Q 2 ) = gi(x,Q 2 )+g 2 (x,Q 2 ) = uW 3 (x,Q 2 ) , (2.15) 

where gi(x, Q 2 ) and gr(x,Q 2 ) are known as the longitudinal and transverse polarized 
structure functions respectively. Here x = is the Bjorken variable, v = P ■ q is 
the energy transfer and Q 2 = — q 2 . The longitudinal and transverse polarization vector 
components are given by 

S-q 



Si, S, 



flT ) 



Sfj, - 



v 



(2.16) 



The hadronic tensor is related to the forward virtual-photon hadron Compton scat- 
tering amplitude as 

= — ImT^, (2.17) 
2n 

where T^ v is given by 

= if d\ e iq< {PS\T{J^), r{0))\PS) 



( - <T + ^)T 1 (x, Q 2 ) + - (f - ^)t 2 (x, Q 2 ) 



-le^qx 



S a T 3 (x,Q 2 ) + P a S-qT^x,Q< 



(2.18) 



It can be shown that the most dominant contribution to W^ v comes from the region near 
the light cone, £ 2 0, at high value of Q 2 . So one can expand the products of currents 
near the light cone. For the sake of clarity and simplicity we suppress the Lorentz indices 
and write 



T(z/,Q 2 )^ 



i,n 



E 2 ^"- 2 ^9(Syr / d 4 ^C[\e)(P I Of ~"» | P>. (2.19) 



d(iq 2 ) 



Here, in the last step we have replaced £ M s by ^-s and used the relation 
d 



d d n 



d(iq 2 
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trace terms. 



(2.20) 



The trace terms are neglected since they are suppressed. So we get 

T(v,Q 2 ) = £(^r% 1 ---^C'r(Q 2 )<P | Of -"-(O) I P>. (2.21) 

i,n 

Here C™(Q 2 ) are (Fourier transform of) the Wilson coefficients and Of "'^ n (0) are local 
operators. According to their most general Lorentz structure, the spin averaged matrix 
elements of the operators Oi can be written as || 

(P | Of -^(0) | P) = A?(p' 11 ...p Mn - mV'T-l (2-22) 

The terms proportional to g'^ 2 are trace terms which are subtracted to project out 
operators of definite spin. We then get 

2 

Ttl 

T(u,Q 2 ) =Y.CUQ 2 )x~ n A 1 t + 0[x- n+2 —]. (2.23) 

i,n ^ 

The neglected terms are either contributions from the trace terms (which are called target 
mass effects) or they arise from higher twist operators (in this case, the mass scale m 2 is 
not necessarily the target mass). 

The Wilson coefficients Cf{Q 2 ) depend on the renormalization scale \i and the coupling 
g(fi). They obey the renormalization group equation 

d d O 2 

■£ + 

where (3 is the QCD beta function and 70™ is the first non-trivial (one loop) order anoma- 
lous dimension of the operator Of [T|. The leading Q 2 dependence to all orders in a s 
can be obtained from the solution of the RG equation. One then uses optical theorem 
to connect T^ v to W^ v . This is not straightforward, since the physical region for DIS is 
< x < 1. But from Eq. (|2.23| ) one can see that T^ u diverges in the physical region. 
This difficulty is overcome by analytically continuing T^ u to complex x and get the result 
corresponding to the physical region as a limit of the analytic function in the unphysical 
region. In this way one automatically gets the x moments of the structure functions. We 
define T(x,Q 2 ) as a function of complex x which is analytic as | x |— > 00 and has a cut 
from —1 to +1. The coefficient of x~ n in Eq. ( |2.23|) is isolated by taking the x n ~ l moment 
and integrating along a closed contour in the complex x plane such that the cut lies inside 
the contour. We then get 

J c dxx n - l T{x, Q 2 ) = ]T CnQM?- (2-25) 
11 



0"— + - To?)CT(^, g(ji)) = (2.24) 



This is called Mellin transform. Shrinking the contour to the physical cut, and also using 
optical theorem, we get 



This equation uniquely predicts the Q 2 evolution of the DIS structure functions W(x, Q 2 ) 
provided we know the non-perturbative part, Af. These depend on the target hadron 
state and are related to the experimentally measured parton distribution functions. One 
takes the ratio of moments at different Q 2 to get rid of the A™s and compare moments 
at various scales of Q 2 . The moments can be inverted by what is called an inverse Mellin 
transformation to get the structure functions. 

So far, we considered only lowest twist contributions. Higher twist effects are sup- 
pressed by powers of compared to the lowest twist effects in DIS. They may become 
important at lower value of Q 2 where perturbative QCD is no longer applicable since the 
coupling constant becomes large. In this kinematical region, the effect of target mass 
also becomes important. Correction terms that appear when the target mass effects are 
included involve powers of and are often referred to as kinematical higher twist effects, 
whereas terms coming from higher twist operators are called dynamical higher twist ef- 
fects. 

Operators of higher dimensions contribute at higher twist. The detailed analysis of 
DIS using OPE keeping only the twist two contributions can be found in textbooks [[I], || . 
The analysis involving the higher twist operators is very complicated because there are a 
large number of operators. 

The target mass enters in the structure functions in the OPE method through the 
subtraction of the trace terms. The operators in OPE are made traceless to project out 
operators of definite spin. The trace terms are proportional to ^ which are subtracted. 
These are higher twist effects and in the leading twist analysis they are neglected. Target 
mass effects have been analyzed by Natchman by introducing a variable £ defined as, 



The basic assumption in this case is that the effective coupling constant is small such 
that free field OPE is a good approximation. Instead of scaling in terms of the Bjorken 
variable x, the structure functions scale in terms of the variable £. In other words, if the 




(2.26) 



2x 



(2.27) 



1 + (1 + 



4M 2 x 2 \ 2 
Q 2 > 
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structure functions F(x, Q 2 ) are expressed in terms of the variables £, such that, 

F=(^r 1 F, (2.28) 

the F(£, Q 2 ) is independent of the scale Q 2 in free theory. £ scaling was derived by 
Natchman || by expanding the products of currents in terms of operators of definite spin. 
He showed that, the light cone expansion of the product of currents in Mincowski space 
corresponds to the expansion in terms of O4 spherical harmonics. These are orthogonal 
Gegenbauer polynomials on the interval — iQ < v < iQ. The orthogonality relations 
among the spherical harmonics project out operators of definite spin n. An alternative 
way of deriving £ scaling is given in |J. Violations of £ scaling occurs due to interactions. 

£ scaling can be interpreted in the language of the parton model in the following way. 
We take the initial proton four-momentum as (p°,0,0,p 3 ) and the virtual photon four- 
momentum as (q°, 0, 0, q 3 ). Then taking both the initial and the final state of the struck 
quark to be on the mass shell, it can be shown that the expression for £ given by Eq. 
( P-27| ) can be written as 

£ = 4^4 ( 2 - 29 ) 

where pi is the initial momentum of the struck quark and we have assumed that p\, — 0, 
where pt is the part of pi transverse to p and q. 

As we mentioned before, higher twist contributions to DIS also come from the higher 
twist operators. These dynamical higher twist effects involve nontrivial quark-gluon in- 
teraction. There are two approaches to analyze the higher twist contributions: 

(1) OPE is used to relate the O(q^) corrections to the target matrix elements of certain 
classes of twist four local operators [[7]. 

(2) The other method is a graphical method to express O(qs) corrections in terms of 
certain quark-gluon-target forward scattering amplitudes [[J. 

Both of these methods are very complicated. The starting point is the assumption that 
the structure function can be factorized into two parts: the coefficient functions Cj which 
describes the short distance part of the interaction and can be calculated in perturbation 
theory, and the parton correlation functions fi which contain non-perturbative information 
about the distribution of partons inside the struck hadron. 

OPE deals with the moments of the structure functions. In the QCD improved parton 
model analysis one gets the parton correlation functions as the Fourier transform of the 
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matrix elements of light-front bilocal currents. One cannot explore these expressions fur- 
ther because of the non-trivial hadronic matrix elements. Here one has to parametrize the 
nonperturbative information and it is only the Q 2 evolution which is calculated pertur- 
batively. We have followed a recently developed alternative approach based on light-front 
Hamiltonian technique. It addresses the structure functions directly instead of its mo- 
ments, using the Bjorken- Johnson-Low (high energy) expansion of the Compton scattering 
amplitude. More important is the point that our approach gives a complete description of 
the structure functions because it has the potential of incorporating the non-perturbative 
contents of the structure functions. It is possible to investigate the structure functions for 
a bound state target by solving the light-front QCD bound state equation. At present, 
active research is going on for solving light-front QCD bound state equation using a tech- 
nique called similarity renormalization (see chapter 3 for references). Another advantage 
of our method is that the target and quark mass effects are automatically included and 
one does not have to introduce them by hand. The intrinsic picture of parton model is 
there but the partons are massive, non-collinear and interacting. 
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CHAPTER 3 



Basic Features of Light-Front Field 

Theory 

In this chapter, we introduce the main features of light-front field theory and light-front 
QCD. 

The light-front formulation of a dynamical system was originally introduced by Dirac 
in 1949 He called it 'Front Form'. Simply speaking, light-front dynamics is the 
description of the evolution of a relativistic system along a light-front direction |J. 

In the usual (equal-time) or instant formulation of field theory, the initial conditions 
are fixed on a spacelike hypersurface, x° = 0, which in the quantum theory becomes 
the surface of quantization. The time evolution of the system is given by the equal 
time Hamiltonian. In the front form, the initial conditions are fixed on the light like 
hypersurface, x + = x° + x 3 = 0, which is the surface of quantization. It is a surface 
tangential to the light cone. The light-front time is x + and the x + evolution of the system 
is governed by the light-front Hamiltonian, P~ . The light-front coordinates are, 

x ± = x°±x 3 , x ± = (x 1 ,x 2 ). (3.1) 

Light-front time derivative is, d~ = 2^-, the space derivatives are, d + = and 
d 1 - = — -Jp-. The metric and other useful relations are given in Appendix A. For a 
system of particles moving with relativistic speed, the particle world lines cluster near 
the light cone, in contrast to the non-relativistic situation where the particles move with 
a velocity much less than the velocity of light c and the world lines cluster near the time 
axis. If we choose x + as the new time direction, the structure of the time evolution of 
a relativistic system may look similar to a non-relativistic system in the ordinary space 
time coordinate system. As a result, relativistic dynamics is expected to become simple 
and more transparent || . 
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It may be mentioned here that, there exists another formulation of field theory, known 
as the infinite momentum frame formulation |3j, where the equal-time results are boosted 
to an infinite momentum frame. Although the Lorentz transformation is highly singular, 
the singularity cancels in the physical objects like Poincare generators. It can be shown 
that such a transformation basically produces the coordinate change, 



x 



± -x°±x\ x ± = (x\x 2 ) (3.2) 



which are the same as in light-front field theory. However, in light-front formulation, 
we start with these coordinates and no infinite boost is involved. One works in an arbi- 
trary reference frame in contrast to the infinite momentum frame formalism where one is 
restricted to one particular reference frame. 

Deep inelastic lepton-hadron scattering is a highly relativistic process. As described 
earlier, Feynman introduced the parton model in the infinite momentum frame which 
successfully described scaling of the structure functions. Light-front formulation of QCD 
gives an intuitive picture of DIS by keeping close contact with the parton ideas, however, 
the partons here are field theoretic partons and one obtains the expected scaling violations 
and many other interesting results. 

3.1 Some Special Features of Light-Front Field Theory 

3.1.1 Dispersion Relation 

Using the metric tensor given in Appendix A, we can see that the dot product of two 
4- vectors is given in terms of light-front components as, 

1 i 

a-b = = -a + b- + -a~b + - a 1 - ■ b ± . (3.3) 

2 2 

So the dispersion relation for an on-mass-shell particle is given by, 

k- = ( ^ + + '" 2 (3.4) 

where k~ is the light-front energy conjugate to light-front time x + and k + is the longitu- 
dinal momentum. The above dispersion relation is remarkable in the following aspects: 

(i) Even though it is a relativistic dispersion relation, there is no square root factor 
involved. This is helpful if one tries to solve eigenvalue equations like, H | ip) — E \ ip) . 

(ii) The dependence of the energy k~ on the transverse momentum k x is just like in 
the non-relativistic dispersion relation. 
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(iii) From the dispersion relation, we see that for an on-mass-shell particle, for k + 
positive (negative), k~ is also positive (negative). It so happens that the particles with 
negative energy which must have negative longitudinal momentum can be mapped into 
antiparticles with positive energy and positive longitudinal momentum. As a result, we 
always have k + > 0. This fact has several interesting consequences, particularly in the 
context of the vacuum structure in light-front theory, which we shall discuss shortly. 

(iv) The dependence of energy on k x and k + is multiplicative and large energy can 
result from large k x and/or small k + . This has drastic consequence for the renormalization 
aspects of light-front field theories ||. 

3.1.2 Poincare Generators in Light-front Formulation 

A relativistic dynamical system must be Lorentz invariant. The inhomogeneous Lorentz 
transformation is generated by ten operators, Hamiltonian, three momenta, three boosts 
and three angular momenta. These generators satisfy Poincare algebra, which is given in 
terms of 4-momenta P p and generalized angular momentum as, 

[P^ p v ] = 0, [M^, P p ] = i(-g pp P u + g up P p ), (3.5) 

[M pu , M pa \ = i{-g pp M ua + g vp M pa - g pa M pu + g ua M pp ). (3.6) 

In the usual equal time or instant form of dynamics, the rotation and boost operators 
are given by, = €ijkJ k and M 0t = K % respectively. Four out of the ten generators, 
namely, the Hamiltonian and the three boost generators, are dynamical (depend on the 
dynamics), the remaining six are kinematical (do not depend on the dynamics). The 
action of the dynamical generators change the surface of quantization whereas the action 
of the kinematical generators do not. In field theory, the ten Poincare generators are 
constructed from the symmetric energy momentum density 9^: 

P^ = J d 3 x6 0p , (3.7) 

M pv = J d 3 x[x p 6 0u - x u 9°% (3.8) 

In light-front field theory, the ten Poincare generators are obtained in the same way, from 
the energy momentum density; 

P ,M = -J dx-d 2 x ± 9 +p , (3.9) 
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M = -J dx-cPx 1 - \x»6 +v - x u 9 +f1 } . (3. 10) 

Here, one has to remember that x + is the light-front time. P~ is the light-front Hamil- 
tonian, this produces the x + evolution of a system. P + is the longitudinal momentum 
and P x are the transverse momenta. M + ~ = 2K 3 and M +t = E l are light-front boosts. 
M 12 = J 3 and M~ % = F l are the three rotation operators. 

A very important aspect of the light-front theory is that, here, the boost operators 
are kinematical. Moreover, the longitudinal boost generator K 3 produces just a scale 
transformation and the transverse boosts behave like the non-relativistic Galilean boosts 
H]. On the other hand, the transverse rotation operators, F l are dynamical which means 
that they change the surface of quantization, x + = 0. So the kinematical subgroup of the 
Poincare group is enlarged in the light-front theory and instead of four dynamical gener- 
ators we have only three, which are the light-front Hamiltonian and the two transverse 
rotation generators. Because of this different nature of the Poincare generators, certain 
problems are easier to tackle in the light-front formulation of field theory. For example, 
in any practical calculation using the Fock space expansion of the state, one has to do 
particle number truncation (called Tamm-Dancoff truncation). In the equal time theory, 
such a truncation has a major difficulty since this violates boost invariance. In the light- 
front theory, on the other hand, such a truncation violates rotational invariance. It is 
argued that the restoration of rotational invariance is easier than the restoration of boost 
invariance because the rotational group is compact in contrast to the boosts which are 
non-compact |J. This makes the light-front Hamiltonian technique more advantageous 
than the equal time Hamiltonian techniques. Because of this, the highly complex QCD 
bound states may be represented by a few boost invariant multi-particle wave functions 
in the Fock space expansion, as we shall see later in this section. 

3.1.3 Simplicity of Light-front Vacuum 

One of the most remarkable features of the light-front field theory is the so called 
triviality of the vacuum. In the equal time formulation of QCD, vacuum is highly non- 
trivial and is the source of various non-perturbative effects like chiral symmetry breaking, 
confinement etc. In light-front theory, a particle's momentum is separated into two com- 
ponents, longitudinal component which is along the light-front time direction x + and 
transverse components which are perpendicular to x + direction. As we said before, in the 
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light-front theory, longitudinal momenta k + > 0. Now, vacuum is an eigenstate of the 
longitudinal momentum, P + | 0) = 0. So the vacuum of a light-front theory can only have 
particles with zero longitudinal momenta, since the total P + is given by, P + = 
where kf are the longitudinal momenta of individual particles in the state. At this point, 
the vacuum is still nontrivial. To construct hadronic bound states in the Fock space of 
quarks and gluons, one needs to simplify the vacuum. If we impose a small longitudinal 
momentum cutoff e in the theory, light-front vacuum becomes trivial, i.e. it does not con- 
tain any particle. It is important to understand that, this trivial vacuum is the eigenstate 
of the full interacting Hamiltonian. In the equal time case, there is no such restriction 
on the longitudinal momenta k 3 which can be positive, zero or negative. So even if we 
impose a cutoff k 3 > e, the total longitudinal momentum P 3 = J2ikf, can have zero eigen- 
value (vacuum) for non-zero kf and we do not get a trivial vacuum devoid of any particle. 
Because of this triviality of the vacuum in light-front theory, Hamiltonian perturbation 
theory becomes simpler since we no longer have the graphs that generate particles from 
the vacuum (Z graphs and vacuum diagrams). All such diagrams require negative kf 
which is not possible. In the case of QCD, the effects associated with the non-trivial 
vacuum in the equal time formulation are expected to be recovered while removing the 
cutoff dependence from the theory, i.e. by the process of renormalization which produces 
new interactions in the theory. The wave function of a relativistic bound state can be 
expressed as an ordinary Fock space expansion using the trivial vacuum, 



For QCD, f{a)) should consist of quark, antiquark and gluon creation operators with 
non-zero longitudinal momenta and it must be a color singlet operator. 

3.2 Light Front Bound State Equation 

Due to the simplicity of the light-front vacuum, the hadronic bound states can be 
expanded in terms of the Fock space components. Explicitly, a hadronic bound state of 
momentum P and helicity A can be written as J7|, 






n,\i 
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Here n represents n constituents in the Fock state | n,x i P + ,XiP ± + K^,Xi), Aj is the 
helicity of the ith constituent, /' denotes integration over the space, 

5> f = i, £«+ = o, ( 3 - 13 ) 

i i 

where Xi is the fraction of the longitudinal momentum carried by the ith constituent, kj- 
is its transverse momentum relative to the center of mass frame, 

*i=fq:, K i=pt-x i P\ (3.14) 

Here pf and pf~ being the longitudinal and transverse momenta of the zth constituent, 
nj-, Xi) is the amplitude of the Fock state | n, XiP + , XiP 1 - + nf, Xi) and satisfies the 
normalization condition, 

Y, f dxid 2 Kj\ <5> x n (x u ^,Xi) | 2 = I. (3.15) 

n,Xi 

The states obey the relativistic Schroedinger-like equation, 

(P^) 2 + M 2 

P-\P,X}= { -^ \P,X). (3.16) 

P~ is the light-front Hamiltonian. One remarkable property of the light front is that 
boost does not involve dynamics, so if one can solve the bound state equation in the rest 
frame, one can understand the particle structure in any frame. Solving the bound state 
equation for the full Fock space is, of course, impossible, but one can restrict oneself to 
a few particle sector (Tamm-Dancoff truncation) || and get useful information of the 
theory. It is important to mention here that Tamm-Dancoff truncation in the usual equal 
time theory encounters a major problem since it violates boost invariance. But in the 
light-front theory, since boost is kinematical, it breaks only rotational invariance. It is 
argued that the restoration of rotational invariance is easier than the restoration of boost 
invariance since rotational group is compact whereas the boosts are non-compact 0. 

The feasibility of describing hadronic bound states in terms of Fock states gives a hope 
of reconciling constituent quark model, which has been so successful in describing hadron 
spectra, with QCD. We shall see in the later chapters, how this particularly helps us to 
investigate DIS structure functions. 
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3.3 Light-Front QCD 



QCD (Quantum Chromodynamics) has now been taken as the fundamental theory of 
strong interactions. Light-front QCD, which means the theory quantized on the surface 
x + = 0, is basically used to explore non-perturbative effects of QCD and is based on 
Hamiltonian theory. The technique was first developed by Kogut and Soper || for QED. 
Because of the triviality of the vacuum and kinematical boost, hadronic bound states 
can be studied in a truncated Fock space of quarks and gluons. The low energy hadron 
structures like the parton distribution functions, fragmentation functions and various 
hadronic form factors can be addressed in this framework. Perturbative calculations are 
also possible using the old fashioned light-front Hamiltonian perturbation theory, which, as 
described before, is simpler than the equal-time Hamiltonian perturbation theory because 
of less number of graphs at each order. 

Before the advent of QCD, Feynman proposed the parton model in the infinite momen- 
tum frame. Since the infinite momentum frame formulation is equivalent to the light-front 
formulation, the natural framework to investigate parton phenomena in high energy scat- 
tering process is the formulation of QCD on the light-front in light-front gauge. Tomboulis 
TO first proposed Yang-Mills theory on the light-front in 1973. A systematic formulation 



of light-front QCD was given by Casher []TTJ and Bardeen et al \T2\. Thorn [I3| studied 
various aspects of light-front QCD including asymptotic freedom for pure Yang Mills the- 
ory. Leapage and Brodsky [|l4j] applied it to the study of exclusive processes in the early 
eighties. However, only very recently there has been a lot of interesting investigation and 
development in light- front QCD |15[]. In particular, active research is going on in the 



area of recently developed similarity renormalization (16 
The QCD Lagrangian is given by 

£ = ~Tr[F^F^\ + ^(i^D» - ™)V> (3.17) 

where F^ = d»A v - d u A fi - ig[A' x ,A u ], N 1 = T,aK Ta is the S luon field color matrix 
and T a are the generators of the SU(3) color group. The field variable ip describes quarks 
with three colors and Nf flavors, = <9 M — igA^ is the covariant derivative, m is the 
quark mass (matrix). The equations of motion are : 

g^p^a + g F abc AlF^ c + gi>YT a ^ = 0, (3.18) 

(i^-m + g^A^^O. (3.19) 
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Dynamical and Constrained Fields: The fermion spinor in light-front coordinates 
can be split into if) = ip + + ip~ , if)^ = A ±r ip where A ± = |7°7 ± (see Appendix A). It 
turns out that, the fields ip~ and A~ ( minus component of A^) are both constrained, the 
equations of constraints are given by, in the light-front gauge A + = 0, 

id+if)~ = [a 1 - ■ {id x +gA L ) + 7 °m]^ + , (3.20) 
-d + A~ a = d l A ia + gf abc ^-(A ib d + A ic ) + 2g^-(f ] T a i\ (3.21) 

where d + = It is a unique feature of light-front formulation of any field theory 

that some of the fields which are dynamical in the ET (equal time) formulation become 
constrained here, as a result the equations of motion of the dynamical fields ip + and A 1 - 
are first order in light-front time. In order to solve the constraint equations, it is needed 
to define the operator Assuming antisymmetric boundary condition, it is defined as 

m 



^:/0O = \ /"J dy~e(x- - y-)f(y-). (3.22) 

Choosing the antisymmetric boundary condition completely fixes the residual gauge free- 
dom. 

The Hamiltonian can be obtained from the Lagrangian. Using the constraint equa- 
tions, the constrained field can be removed from the Hamiltonian and it can be written 
entirely in terms of dynamical variables in light-front gauge. The light-front QCD Hamil- 
tonian in terms of the dynamical variables is given explicitly in Appendix C. There is 
another method of constructing light-front QCD Hamiltonian which is based on light- 
front power counting ||, we shall not discuss it here. 

Quantization: In order to quantize such a constrained system, one has to follow 
Dirac procedure [0 or the phase space quantization method [15, 17 1. The non- vanishing 



commutators between the dynamical fields are given by, 

{if) + (x),if)i + (y)} x+=y+ = A + S 3 (x-y), (3.23) 
[Ai(x), A{{y)\ x+=y+ = -i8 ab 8v\e{x- - y~)5\x ± - (3.24) 

where 5 3 (x — y) — 5(x~ — y~)5 2 (x ± — y L ). It is to be noted that these commutators are 
evaluated in equal light-front time, x + = y + . The ET (equal-time) commutators between 
two field variables must vanish because of causality since the separation is spacelike, 
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however, the light-front commutators are non-vanishing because here the two fields are 
separated by a light like distance, and a signal traveling with the velocity of light can still 
communicate between them. 

Non-locality: Another interesting point to note is that, the above commutators and 
anticommutators are non-local in the longitudinal (x~) direction. This non- locality creeps 
in since in light front coordinates, x 2 = x + x~ — (x -1 ) 2 , and x + = still allows non-local 
x~. This non-locality in the longitudinal direction is also seen in the definition of and 



has far reaching consequences ||19|| . 



Two Component Formulation: By eliminating the constrained fields completely 
in light-front gauge, tht light-front QCD Hamiltonian can be written completely in terms 
of the dynamical field components, ip + and A- 1 . A very interesting point is that, in a 
particular representation of the gamma matrices, (see Appendix A) the four component 
fermion field can be reduced to a two component field, i. e. the dynamical field ip + can 



be written as [17 







(3.25) 



where £ is the two component fermion field. This is one of the major advantages of the 
light-front theory that simplify the relativistic fermionic structure. The two component 
fermionic field satisfy the anticommutaion relation, 

{^£}} x+=y+ =5\x-y). (3.26) 

In the two component formalism, the Fock space expansion of the dynamical fields become, 

AI{X) = ^2^ [ela{q ' X)el9X + kX - ] > (3 ' 27) 



2(2tt)V 

= E xx J XX P r + [b{p > xyQX + rft °°' X)e ~ iP% (3 ' 28) 



dp + d p 



where q = g+ and p = p+ — . A can be 1 or —1 for the gluon and ^ or — ~ for the 
quark. The gluon polarization vectors are, 

4 = i=(l,z), e i 1 = _L(i,_,). (3.29) 

The quark two-component spinors are 

XI = (1,0), X _i=(0,l). (3.30) 
24 



[a(q, A), a\q\ A')] = 2(2tt)V * 3 (? - q')5 xx ,, 
{b(p, A), b\p>, A')} = {d(p, A), S(p', A')} = 2(2vr) V5 3 (g - q'). (3.31) 

As described before, the constrained fields are eliminated from the theory using the con- 
straint equations and everything can now be written in terms of the two-component 
dynamical fields, £ and A x . This makes computations easier. 

Power Counting and Renormalization: Power counting in light-front theory is 
different from the usual equal time case. In order to understand this, let us look into the 
dispersion relation, 

k~ = {kl) l + m2 (3.32) 

which means that k x scales as m and the product of k + k~ scales as m 2 , but k + and 
k~ individually can have different scaling behavior. Among the light-front coordinates, 
only the transverse directions x x carry inverse mass dimensions, while the longitudinal 
coordinate x~ does not carry any mass dimension ||. This makes the crucial difference 
from the instant form where all directions carry the same mass dimension. 

In perturbative calculations using field theory, one often encounters divergences. In 
light-front field theory, one has to treat transverse and longitudinal directions separately in 
determining the superficial degree of divergence of a divergent integral by power counting. 
Also since light-front theory is manifestly non-covariant and gauge fixed, one gets new 
types of divergences. However, the physical quantities themselves are finite. One has 
to regulate these divergences and remove them, this process is called renormalization. 
Renormalization in light-front field theory is completely different from the usual instant 
form of field theory. One has to add non-local non-covariant counterterms to restore all the 
invariances in the renormalized theory. Various investigations regarding renormalization 
of light-front theories are in progress. Recently a special renormalization scheme, called 
similarity renormalization has been proposed |2(| which is suitable for the bound state 



studies on the light-front. In this work, we have used ultraviolet transverse momentum 
cutoff and a small infrared longitudinal momentum cutoff. 
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CHAPTER 4 



Light-Front Hamiltonian QCD 
Approach to DIS 

A large number of recent experiments on polarized and unpolarized deep inelastic 
scattering are providing valuable information on the so called 'higher twist' or power 
suppressed contributions to DIS. A theoretical understanding of these effects require non- 
perturbative information on the structure of the hadron. There is an urgent need to 
develop a non-perturbative approach to DIS which is preferably based on physical intu- 
itions but at the same time employs well defined field theoretical calculational procedure. 
This goal is achieved in light-front Hamiltonian QCD framework. We follow a recently 
developed method of calculating the structure functions which is a combination of co- 
ordinate space approach based on light-front current algebraic techniques and momentum 
space approach based on Fock space expansion method in light-front theory in Hamilto- 
nian QCD. 

The most intuitive approach to DIS is the parton model || proposed by Feynman. 
Partons were originally introduced as collinear, massless, non-interacting pointlike con- 
stituents of the proton. However, the QCD governed interacting partons need not be 
collinear and massless. Can one generalize this concept to introduce field theoretic par- 
tons, which are non-collinear and massive, but still on-mass-shell objects in an interacting 
field theory? The answer is yes, and this is achieved in light-front Hamiltonian QCD by 
introducing many body (or multi-parton) wave functions. 

In order to introduce these many body wave functions in the description of DIS struc- 
ture functions, we use Bj or ken- Johnson-Low (BJL) expansion of the scattering amplitude 
@, which is essentially a non-perturbative approach where the expansion parameter is 
the inverse of the light-front energy of the probe (virtual photon). The leading term in 
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this expansion in the high energy limit contains the matrix element of the equal light- 
front time current commutator which can be evaluated, knowing the light-front commu- 
tators/anticommutators between the field variables. Using optical theorem (see Chapter 
2) one obtains the expressions of the structure functions as the Fourier transform of the 
matrix elements of light-front bilocal vector and axial vector currents. The bilocality is 
only in the longitudinal (minus) direction and the matrix element is between the target 
hadron state. In the next step, one uses the Fock space expansion of the target state using 
Tamm-Dancoff truncation. This introduces multi-parton wave functions in the expression 
of the structure functions. The structure functions, then, can be evaluated once these 
wave functions renormalized at scale Q are known. At present major efforts are under 
way to calculate these wave functions. It is important to emphasize here that, such an 
approach is possible because of certain special properties of light-front formulation, in 
particular, the fact that light-front boost is kinematical (see chapter 3). The advantage 
of this technique is that both perturbative and non-perturbative issues can be addressed 
within the same framework. So, a unified description of perturbative and non-perturbative 
QCD underlying the deep inelastic structure functions can be realized. 

In this section, we briefly outline our approach and summarize the previous works 
done using this formalism. 

4.1 Structure Functions in the BJL Limit 

The hadronic tensor W^ v is given in terms of hadronic matrix elements of the current 
commutator as in chapter 2. However, it is not equal-time or equal-x + commutator. The 
DIS structure functions are related to the equal x + current commutators through the BJL 
(high energy) limit of the forward virtual photon-hadron Compton scattering amplitude 

T*" larg = 9 " -— f dCd 2 ^e iq <(PS\[J^{0, J u {0)y =0 \PS) + ... (4.1) 
q J 

where q~ is the light-front energy of the virtual photon, and (£ + , £",£*) are the light- 
front space-time coordinates. For large Q 2 and large v limits in DIS, without any loss of 
generality we can always select a Lorentz frame such that the light-front energy q~ of the 
virtual photon becomes very large. Explicitly, in terms of light-front variables, we can 
choose q + to be negative and finite for the virtual photon (since the photon is virtual, 
longitudinal momentum can be negative). Also, keeping q % to be finite, one can get large 
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space-like q 2 (Q 2 — > oo) by taking q~ — > oo. In this case, x = — is positive and finite, 
since q + is negative while both q + and P + are finite. So by taking the limit q~ — > oo, we 
do not go beyond the physical region of DIS. Also, in the large Q 2 region, we take only 
the first term in the above expansion since the higher and higher terms are suppressed. 

The structure functions are related to the equal light-front time current commutator 
by using optical theorem and dispersion relation (see chapter 2). The light-front current 
commutator can be calculated exactly from light-front QCD in the light-front gauge. 

It is worthwhile to point out that there exists an equal-time BJL expansion of the 
scattering amplitude where the expansion parameter is -\. But taking q° — > oo gives 
timelike q 2 which is unphysical for DIS. In order to overcome this difficulty, one has 
to go over to the complex q° plane and take iq° — ► oo. This increases the complexity. 
Another important point is, equal-time BJL expansion gives an infinite set of relations 
connecting the moments of the structure functions to the corresponding term in the BJL 
expansion ||. So in order to get the structure functions, one has to invert these equations. 
Light-cone BJL expansion instead gives the structure functions themselves in terms of the 
light-cone current commutators. 

The various expressions of the DIS structure functions that are obtained in the BJL 
limit and inserting the current commutator are given below: 

F2iX,Q2) 1 f dve-^J2el(PS\^ a (rh + M0) 

a 

-^ a (0h + MC)\PS) (4.2) 

a 

-^a(oyAMn\ps) ■ (4.3) 



x AttP + 



AttP] 



Eq. (4.4) may be reduced to the same expression obtained by the collinear expansion 
in the Feynman diagrammatic method up to the order twist-four [||. But it is obtained 
directly here in the leading order in the l/q~ expansion without involving the concept of 
twist expansion. The polarized structure functions come out to be, 

9i(x,Q 2 ) = — / dr ] e-^Y.el(PS\^ a (Ch + ^MO) + h.c.\PS), (4.5) 
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9t(x,Q 2 ) = ^/^e-^^e^P^|^(r)(7 J -^7 + )75^(0) 

+ h.c.\PS). (4.6) 

We refer to |TJ for a detailed derivation of the above results. In our work, these are 
the starting expressions. The detailed derivation of the structure function Fl taking 
into account the flavor structure is given in the next chapter. The above results are 
derived without recourse to perturbation theory, and also without the use of the concept 
of collinear and massless partons. Also, they are not derived using the twist expansion 
and these are the leading terms in the -4r expansion and not the leading terms in terms of 
twist, the target here is in arbitrary Lorentz frame. In the rest frame, these expressions 
reduce to those obtained by Jaffe and Ji in the impulse approximation || and also by 
Efremov et. al. in QCD field theoretic model |J. Also, the above expressions are valid 
only in the light-front gauge, A + = 0, otherwise the bilocal expressions should involve a 
path ordered exponential to ensure gauge invariance. Since the bilocality is only in the 
longitudinal direction, the exponential factor is unity. 

In the above expressions, the bilocal operators arise from the light-front current com- 
mutator. It can be seen from Eq. (|4.2| ) and ( |4.3|) that the plus (involving 7 + ) and 
perp (involving 7- 1 -) components of the bilocal current are related to the same structure 
function F%, whereas the minus component (involving 7") is related to the longitudinal 
structure function Fl. Explicit calculation shows that the minus component involves the 
constrained field, if>~ , which, when eliminated using constraint equation, gives rise to in- 
teraction dependent terms in the operator. This interaction dependence is absent in the 
plus component. Similarly in the polarized structure function g 1; the operator involved is 
the plus component of the axial vector current which does not depend upon interactions 
whereas in transverse polarized structure function g?, the perp component is involved, 
which contains t/>~, elimination of which brings in interaction dependence. The evolution 
of the structure functions F 2 and g\ with Q 2 are generated entirely from the target state 
in the leading log approximation. In these cases, all the complexities are buried in the 
target hadronic state and the operators involved have simple structure. The situation 
is completely different for the other two structure functions, namely Fl and gx, where 
complexity is already there in the operator structure. The interaction dependence of the 
operator gives off-diagonal contributions to the hadron matrix elements. 



31 



The evaluation of the hadronic matrix elements is straightforward in this approach and 
this clarifies the physical picture of DIS. In this framework, we do not use dimensional 
regularization and the presence of 75 does not produce any problem. Also, in contrast 
to the OPE formalism, the quark mass effects can be included without any difficulty. In 
this approach, we deal with probability amplitudes rather than probability densities and 
interference effects are easy to handle. 

In the next section, we briefly review the previous works done in this approach. 

4.2 Review of Previous Works Done in This Approach 

4.2.1 F 2 Prom Transverse Component of Bilocal Current 

Earlier in this section, we have written F 2 in terms of both plus and perp components 
of bilocal current. That the perp component also gives the same F 2 structure function is 
a new result and comes out directly from our approach based on BJL limit and light-front 
Hamiltonian QCD. This is different from [7[], according to which, the plus component is 
twist two and perp component is twist three. So the perp component should not have 
partonic interpretation. Explicit calculation shows that the perp component involves the 
constrained field ip~ . When it is eliminated using the constraint equation, the expression 
for F 2 becomes, 

^ = / dy-e-i p+ y-*(P I e{y)[O m + <V + O,]£(0) I P) + h.c. (4.7) 

h.c. is the Hermitian conjugate and, 

1 



io + 

1 



O g = g-^[A 1 + i a 3 A% (4.8) 

There is an interaction term O g in the operator. However, contribution from this term 
to the matrix element vanishes and one gets the same structure function F 2 as from the 
plus component and it carries the same parton interpretation ||. 

4.2.2 Factorization of Hard and Soft Dynamics 

Any general cross section in hadron physics contains both short distance and long 
distance behavior and hence is not accessible to perturbative QCD. Factorization theorem 
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allows one to separate the two behaviors in a systematic way. According to factorization 
theorem, the cross section can be separated into two parts, the non-perturbative soft 
dynamics factors out of the hard part M. So the hard (short distance) part of the 
cross section is calculated using perturbative QCD and the soft (long distance) part is 
given in terms of a set of operator matrix elements which have to be calculated non- 
perturbatively using some non-perturbative techniques like lattice gauge theory. This 
factorization theorem has not been proved for higher twist contributions to the cross 
section. Recently, a factorization scheme is presented in our approach based on light- 
front Hamiltonian QCD and Fock space expansion I]. The light-front Hamiltonian is 
separated into three parts, depending on the hard and soft transverse momenta k L , 

p- = (P~) h + (P~) m + (P-) 1 (4.9) 

where the soft part {P~) 1 contains only those momenta for which {k^-) 2 < fi 2 and hard 
part {P~) h contains Q 2 > {k^) 2 > fi 2 where fi 2 is some intermediate scale separating the 
hard and soft dynamics. This is called factorization scale. Q 2 is the scale of the virtual 
photon. The low energy effective Hamiltonian can, in principle, be obtained by integrating 
out all modes with {k^) 2 > fi 2 from the canonical QCD Hamiltonian. The mixed part 
contains interaction and mixes the hard and soft partons. The structure function can be 
written as, 

Fi(x, Q 2 ) « / drjexp-^ £ e 2 (PS | i(r)r,i(0) ± h.c. \ PS) (4.10) 
where Tj denote the gamma matrices. The target bound state can be expressed as, 

\PS) = U h \PS,fi 2 ) (4.11) 
where U h = T+exp(-| / "°° dx + [(P-) h + (P~) m ]) and, 

(P-) 1 | PS,^ 2 ) = £ L± | PS, fi 2 ). (4.12) 



p+ 



Then one can write, 



Q 2 ) w / drjexp-^ E E ( PS ' ^ I n i) 

a ni,ri2 

(n 2 | PS^ 2 )( ni | U^$ a (C)Tiip a (Q) ± h.c.]U h | n 2 ) (4.13) 



where | ni), | n 2 ) are complete sets of quark and gluon Fock states with momenta k 2 < fi 2 . 
The hard contribution is described by the matrix element, (m | U^ l [%[! a {^~)Yiip a {Gi) ± 
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h.c.]Uh | n 2 ) which is evaluated in light-front Hamiltonian perturbation theory. The soft 
contribution is characterized by an overlap of the multiparton wave functions, (PS, /i 2 | 
n i)( n 2 I PS,fi 2 )(rii |. For F 2 and gi, \ n\) =| n 2 ) —\ n), only one parton actively takes 
part in the scattering process, the remaining ones are spectators. Such a factorization 
scheme is valid for the F 2 and g\ structure functions. The structure function can be 
written as, 

F,(x,Q 2 ) ^$> 2 [ 'dyP pp ,(y,x,%)q ai (y,^ 2 ) (4.14) 

where P pp \y,x, %) are the leading hard contributions to the structure functions, called 
splitting functions, and q a i{y,l^ 2 ) is called the distribution function which is to be calcu- 
lated by non-perturbative light-front QCD approaches to hadronic bound states. Thus, 
a unified treatment of both perturbative and non-perturbative aspects of DIS is possible 
in this framework. 

4.2.3 F 2 Structure Function, Q 2 Evolution and Sum Rule 

The factorization scheme shows that the dressed parton structure functions are impor- 
tant, since they give informations on the splitting functions and scaling violations. The 
F 2 structure function has been calculated for dressed quark and gluon target states [T^ . 



The state is expanded in Fock space for a dressed quark in terms of bare states of quark 
and quark plus gluon (upto 0(g 2 )) as, 

\P,a) = <pib\P,a) I 0) 



f dktd 2 kj I dkjd 2 kj / 2(27r)3 p + ^3 ( p _ k _ h) 

£m J p{2^fkt J p^fkt V 
<h(P, ° I h, <ri\ h, X 2 )b\h, a 1 )a^(k 2 , A 2 ) I 0) (4.15) 

where a is the helicity, <p 2 is the probability amplitude to find a bare quark of momentum 
k\, helicity <7i and a bare gluon of momentum k 2 and helicity A2 in the dressed quark. 
Explicit evaluation leads to, upto 0(at s ), 

F2(x,Q 2 ) cz s Q 2 1 -|- x 2 3 

—r~ = 6{1 ~ x) + 2-. Cfln ^ [ jr^x)i + 2 6{1 - x)] - (416) 

Here we have taken the state to be normalized to unity. The first term in RHS is the 
parton model result and the next term is the QCD correction. InQ? in the second term 
arises from the transverse momentum integration of the field theoretic partons and it 
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gives rise to scaling violations. The collinear (x = 1) singularity is canceled between the 
real and virtual gluon emission, which is indicated by the plus prescription in the above 
expression. The Altarelli-Parisi splitting function can be obtained from this expression: 

It is known that the F 2 structure function gives the longitudinal momentum sum rule. It 
has been shown that the integral of F% over x is related to the hadronic matrix element of 
the fermionic part of the light-front QCD longitudinal momentum density 9 ++ . Thus, it 
gives the total fraction of the hadronic longitudinal momenta carried by the quarks and 
antiquarks. The structure function for the gluons F$(x, Q 2 ) has been calculated and it 
has been shown explicitly that the longitudinal momentum sum rule is obeyed, 



f 1 dx(F$(x) + FHx)) = 1. (4.18) 
Jo 



This sum rule is verified upto 0(a s ) in light-front Hamiltonian perturbation theory. 

In order to understand the factorization of hard and soft dynamics in this picture, F2 
has been calculated for a meson-like bound state by expanding the target state in Fock 
space |I0| . Upto 0(a s ), one considers the bare states consisting of a quark-antiquark pair 
and quark-antiquark-gluon. Q 2 evolution of F 2 has been calculated and it has been shown 
that it factorizes into soft (0 < k L < /i) and hard (/i < k L < A), where A is the cutoff on 
transverse momenta. 

4.2.4 Polarized Structure Function: gi(x,Q 2 ) and Helicity Sum Rule 

The first moment of the polarized structure function gi(x, Q 2 ) gives a measure of the 
total intrinsic helicity contributions of quark and antiquark in the nucleon. The EMC data 
in the late 80's produced a lot of interest in polarized structure function by confirming 
that Ellis- Jaffe sum rule, according to which about 50-60 percent of the nucleon helicity 
should be carried by quarks and antiquarks, is violated. The EMC data showed that 
only a very small fraction of the nucleon helicity is carried by the charged constituents. 
This gave rise to a lot of confusion in understanding the spin contribution of the nucleon 
and this is what is known as the spin crisis. Later, it has been shown by calculating the 
triangle graph anomaly that the gluon coming from the anomaly graph interacts with 
the virtual photon. This leads to a small quark/ antiquark helicity contribution because 



of certain cancellations ||11|| . However, in the total helicity of the nucleon, in addition to 
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the quark, antiquark and gluon intrinsic helicity, one has to take into account the orbital 
angular momenta of the constituents. In fact, the orbital angular momentum of the 
constituents play a very important role in the nucleon helicity. The role of orbital helicity 
has been investigated recently in our approach [12|. A detailed analysis of light-front 
helicity operator is needed for this purpose. Light-front helicity operator J3 (in the frame 
P 1 - = 0, J 3 is the helicity operator; see chapter 7) can be written in terms of manifestly 
symmetric gauge invariant energy-momentum tensor 9^ u in QCD as, 

J3 = i J dx~d 2 x^\xH+ 2 - x 2 9 +1 \. (4.19) 



Recently it has been shown explicitly JTSJ that J 3 constructed in light-front gauge A + 



0, after eliminating the constrained variables ip~ and A~ using constraint equations, is 
independent of interaction and its form agrees with the naive canonical form, provided 
we assume that the fields vanish at the boundary. J 3 thus constructed can be separated 
into four parts, 

J 3 = J% + Jfo + J 3 9 o + 4 (4.20) 



where Jj i is the fermion intrinsic part, Jj Q is the fermion orbital part, J 3 is the gluon 
orbital part and J gi is the gluon intrinsic part. It has been shown that J% is directly 
related to the polarized structure function gi(x, Q 2 ) measured in longitudinally polarized 
DIS. The helicity sum rule can be written as 

(PS I J 3 I PS) = (PS I J% + J% + J 3 go + J% I PS) = ±\, (4.21) 

that is, the sum of the orbital and intrinsic helicities of the constituents gives the nucleon 
helicity. This has been verified in perturbation theory for a dressed quark a and dressed 
gluon explicitly. 

4.2.5 Polarized Structure Function: g2(x,Q 2 ) and the Violation of Wandzura 
Wilczek Relation in pQCD 

The polarized structure function g2(x,Q 2 ) is the least well known structure function 
in DIS. The contribution of g2(x,Q 2 ) to the cross section is suppressed by a factor of 
compared to the leading contribution in the longitudinally polarized scattering. So it is 
commonly known as higher twist effect. Wandzura and Wilczek showed that |13| g-2,{x, Q 2 ) 
can be separated into two parts, 

g 2 (x, Q 2 ) = g™(x, Q 2 ) + U*> Q 2 ) ( 4 - 22 ) 
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where g2(x, Q 2 ) is a twist three piece and the twist two part is related to gi(x, Q 2 ) through 
the Wandzura-Wilczek relation: 

gr(x, Q 2 ) = -m(x t Q 2 ) + ^ %!(x, Q 2 ). (4.23) 

Jx y 

The twist three part was negligible in their model calculation and the quark mass was 
neglected. Later, there has been a lot of interest in g2(x,Q 2 ). The effect of the quark 



mass has been investigated [14| and it was found that the twist three part is a direct 



quark-gluon interaction effect and plays a very important role fL5 |. The importance of 
quark mass in g2(x,Q 2 ) has been investigated recently in the light-front QCD approach 
fl6[l . The pure transversely polarized structure function is defined as, 

9t(x, Q 2 ) = 9l {x, Q 2 ) + g 2 (x, Q 2 ). (4.24) 

The expression for in terms of the hadronic matrix element of light-front bilocal current 
has been given earlier. After one eliminates the constrained fields using the constraint 
equations, the expression of gx becomes, 

9t(x,Q 2 ) = — L. / dr t e- ir *(PS\(O m + O kx + O g ) + h.c.\PS) (4.25) 

where 

O m = m ^-)Q 2 1± (-L)^ + (0l (4.26) 

O k± = -^( r )g2 (7± _L 7 ± . d ± + 2^) 75 ^ + ( ), (4.27) 

O g = #UDQ 2 (7 ± • A ± (C)^ T 1± ■ 75^ + (0). (4.28) 

qt has been calculated for a transversely polarized quark state dressed with one gluon. 
This state can be expressed in terms of helicity states as, 

I k + , k ± , s 1 } = -±=(| k + , k 1 , T)± | k + } At 1 , I)) (4.29) 

with s 1 = ±ms, and mjj is the renormalized quark mass. For this state, 

9tM >) = - x) + | C/; „|[(1±^ ]} (4 .30) 

which is independent of quark mass. However, it is important to start from massive theory 
and renormalize the quark mass in light-front Hamiltonian perturbation theory. However, 
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starting from a massless theory one gets wrong result for g T . The structure function g x 
can also be calculated for a dressed quark in the same method and it can be verified that 
Wandzura-Wilczek relation is not obeyed in perturbative QCD. It is also easy to verify 
that the Burkhardt-Cottingham sum rule is obeyed, that is, 

f 1 dxg 2 (x,Q 2 ) = 0. (4.31) 
J o 

In the next few chapters, we discuss the works that we have done in this research 
program using this framework. 
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CHAPTER 5 



Twist Four Longitudinal Structure 
Function in Light-Front QCD 

In chapter 4, we have outlined our approach based on light-front Hamiltonian QCD 
to DIS. In this chapter and in the next three chapters we shall discuss the problems that 
we have investigated in this research program. 

In this chapter, we shall analyze twist four longitudinal structure function. 

5.1 Motivation for Studying Twist Four Fl 

The higher twist contributions or power corrections to DIS structure functions play 
a very important role in the SLAC kinematical range |]. The longitudinal structure 
function Fl involves coefficient functions which are zero at zeroth order and the lowest 
twist contributions in this case involves calculations to the same order in a s as the NLO 
results for F 2 @. Thus, twist four contributions to Fl are extremely important. These 
are the first non-perturbative contributions to Fl- The higher twist contributions to the 
unpolarized structure functions F 2 and Fl have been analyzed previously in operator 
product expansion method || and also in the Feynman diagram approach M. Qiu has 
given an alternate method based on special propagators utilizing some unique features 
of light-front coordinates ||. At twist four level, in the OPE analysis, there appears a 
proliferation of operators. The operators are not all independent and they are related 
through the equation of motion and the analysis is highly complex. On the other hand, 
in the Feynman diagram approach, it has been shown that one can make contact with 
light-front current algebra analysis, namely, twist four part of F L can be expressed as 
the Fourier transform of the hadronic matrix element of the minus component of light- 
front bilocal current. It is important to go beyond phenomenological parametrizations 
for a proper understanding of the higher twist effects. However, even after many years 
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of investigations, an intuitive physical picture of the interaction dependence of twist four 
longitudinal structure function is lacking. In our analysis we get an intuitive picture of 
F L . 

Another important problem of current interest is the perturbative aspect of the twist 
four matrix element. Simple power counting indicates that in the bare theory the twist 
four matrix element contains quadratic divergences. Understanding the origin and nature 
of these divergences will be quite helpful in finding procedures to remove them (the process 
of renormalization). 

5.2 Plan of This Chapter 

In this chapter, we show that using light-front Hamiltonian QCD framework, one can 
resolve outstanding issues associated with the twist four contribution to the longitudinal 
structure function. Our starting point is the Bj or ken- Johnson-Low (BJL) expansion for 
the forward virtual photon-hadron Compton scattering amplitude. This leads us to the 
commutator of currents which we present in detail for arbitrary flavors in SU(3). Next 
we consider the specific case of electromagnetic currents and arrive at expressions for the 
twist two part of F 2 and twist four part of F L in terms of specific flavor dependent form 
factors. We identify the integral of Fl ^ with the fermionic part of the light-front QCD 
Hamiltonian density. The consideration of mixing in the flavor singlet channel leads us 
to the definition of the twist four longitudinal gluon structure function and then we find 
a sum rule, free from radiative corrections. 

The sum rule which the physical structure function has to satisfy involves the physical 
mass of the hadron which is a finite quantity. A theoretical evaluation of the sum rule 
which starts with the bare theory, on the other hand, will contain various divergences 
depending on the regulator employed. In order to compare with the physical answer re- 
sulting from the measurement, we need to renormalize the result by adding counterterms. 
For the dressed parton target, for example, these counterterms are dictated by mass coun- 
terterms in the light-front Hamiltonian perturbation theory. For a dressed gluon target, 
calculations show that quadratic divergences are generated and one does not automati- 
cally get the result expected for a massless target. The divergence generated is shown to 
be directly related to the gluon mass shift in old fashioned perturbation theory. To a given 
order in perturbation theory, counterterms have to be added to the calculated structure 
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function. The precise selection of counterterms is dictated entirely by the regularization 
and renormalization of the light-front QCD Hamiltonian. The choice of counterterms in 
the Hamiltonian, in turn, determines the counterterms to be added to the longitudinal 
structure function which results in a theoretical prediction of the physical longitudinal 
structure function. We recall that in Hamiltonian perturbation theory we cannot auto- 
matically generate a massless gluon by the use of dimensional regularization. The point 
we emphasize is that the twist four longitudinal structure function is one to one related 
to the Hamiltonian density and that there is no arbitrary freedom in this relationship. 

To understand the nature of quadratic divergences, we evaluate the twist four longi- 
tudinal structure functions for quark and gluon target each dressed through lowest order 
in perturbation theory. The sum rule allows us to relate these divergences to quark and 
gluon mass corrections in QCD in time-ordered light-front perturbation theory. Finally 
we discuss the relevance of our results for the problem of the partitioning of hadron masses 
in QCD. 

The plan of this chapter is as follows. In Sec. 5.3 we derive the expressions for the 
twist two structure function F 2 and the twist four longitudinal structure function Fl using 
the BJL expansion and equal time (x + ) current algebra and also taking into account the 
flavor structure. The flavor singlet expressions are given in Sec. 5.4. The sum rule for Fl 
is derived in Sec. 5.5. In Sec. 5.6 we evaluate F L for a massless quark dressed through 
lowest order in perturbation theory and explicitly verify the sum rule. Then we evaluate 
Fl in perturbation theory for dressed quark (massive) and dressed gluon target. In Sec. 
5.7 we discuss the issue of the breakup of hadron mass in QCD in the context of our sum 
rule. Finally, in Sec. 5.8, we discuss the implication of our results. Our notations and 
conventions have been clarified in Appendix B. We also refer to our original papers, || [7J. 

5.3 Longitudinal Structure Function: Flavor Structure 

In this section we present the expressions for structure functions for arbitrary flavors 
in SU(3) which follow from the use of the Bjorken- Johnson-Low expansion and light-front 
current algebra. In terms of the flavor current J£ (x) = ip(x) , y^^Yijj(x), the hadron tensor 
relevant for deep inelastic scattering is given by 

W2? = ^f*t e^<(P\[J^),^m\P)- (5-1) 
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The forward virtual photon-hadron Compton scattering amplitude is given by 

d^<(p\T{j^)jm)\p)- 



rpflU 



We have 



Tab ( x > 



°° dq)+ W%(x',Q*) 



qt+ _ q-\ 



(5.2) 
(5.3) 



Using the BJL expansion ||, we have 

T£ = ~f dC<P£L<P*(P\[J£{t), ^(o)] e+=0 |P) + ... (5.4) 

q J 

where ... represents higher order terms in the expansion which we ignore in the following. 
In the limit of large q~ , from Eq. ( |5.1| ), we have 



W. 



ab 



L(ab) 



+ 



XV 



2(a6)j 



with x = -g 2 - and v — P.q. On the other hand, from Eq. 



Limit n 



+oo J- ab 



~ I drd 2 ^<(p\[j:(o,J b -(o)y=o\P)- 

(J J 



(5.5) 



(5.6) 



The components of the flavor current </„ (x) obey the equal- x + canonical commutation 
relation (to be specific, we consider SU(3) of flavors) 



x = y 



Zifabc ip(x) 7 y^(x) 5 2 (x ± - y 1 ) S(x -y~ 



- - d + 
2 x 



e(x - y ) if abc V c (x \ y) + id abc V c (x \ y) d 2 (x 1 - y 1 ) 



S'ix^-y 1 -) V c (x\y) - Al(x\y) 



+ - idabc e(x~ 



5 2 {x ± - y L ) V c (x | y) + e« A>(x | y) 



- y-) 91 

In deriving the above relations, use has been made of the relation 

A a Af, = ifabc^c + d a b c \ c . 

Here a, 6, c run from to 8. We have defined the bilocal currents as follows. 



(5.7) 



(5. 
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y7^(y) + ^(2/)y7^(x) 
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(5.9) 



Further, we introduce the bilocal form factors 



(PI 


v^(eio)|p) = 




(5.10) 


(PI 


v:(eio)|p) = 




(5.11) 


From Eqs (|5.6|) and ( 


57|), we get 








Limit q-^oo q 


- T+- = -2if ahc p-T c 





+V / dCe^e(r) f abc (P I V-(£ | 0) | P) + d abc (P | V c "(£ | 0) | P) 



q 
"2 



dfe* 9 r e(D fabc (P I V c (£ | 0) | P) + d afcc (P | V c (£ | 0) | P) . (5.12) 



It is to be noted that matrix elements of A^(x \ y) do not contribute to unpolarized 
scattering. Using the dispersion relation given in Eq. (|5.3| ), together with Eqs. (|5.5| ) and 
( |5.12| ) and comparing the coefficient of q l on both sides, we get 



2(ab) 



X 



X 



(5.13) 



Comparing the coefficients of q + on both sides, we get 
1 i (q+) 2 



L{ab) 



Q 2 7T P^ 

{P L ? 



d V e^ x [fabc{P I V;{t | 0) | P) + d abc (P | V c (£ | 0) | P) 
drje- i ^\f ahc (P\V+^\0)\P) 



-x 



Q 2 ixP 
+ d abc (P\VUt\0)\P) 



(5.14) 



We have introduced r\ = |P + £~. 

Note that our result for Fl differs from the one given in the literature ||. The 
difference can be traced to the expression for F L that one employs. It is customary [|] 
to ignore target mass M 2 in the expression for Fl (see Appendix B). This leads to an 
incorrect expression for Fl which in turn will lead to an incorrect sum rule (see section 
5.5). 

The electromagnetic current 

1 



J"(x) = J£(x) + -j=JS(x). 



(5.15) 



From the flavor structure of electromagnetic current, we observe that, only d abc contributes 
to the structure functions in deep-inelastic electron-hadron scattering. Explicitly, we have, 
F 2 (x) i 



x 



2ttP- 



J dr]e- inx (P | V + (£ | 0) | P). 



(5.16) 
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The longitudinal structure function is given by 

Fl(x) = l±(£ljdT,e-**(P\Tr(t\0)\P) 



" 2 ^"^ 2 / dr l e ~ inX ( P I I 0) I P )- ( 5 - 17 ) 



Q 2 ixP 

We have defined the functions 



V (e 1 0) = UJ V ° (e 1 0) + 3 V ^^ 1 0) + 37I V ^ 1 0) ' (5 ' 18) 
In arriving at our final results we have used explicit values of the structure constants 
of SU(3), 

11 l~2 

^338 — ^888 = — ~^=> ^330 = ^880 = y 3' (5.19) 

Vq is the flavor singlet component of the fermion bilocal vector current. 

5.4 Flavor Singlet Case 

In the flavor singlet channel, for simplicity, we do not explicitly write the flavor charge 
dependence. The BJL limit, together with light-front current algebra || in A + = gauge, 
the tools used in the pre-QCD era, lead to the twist four part of the fermionic contribution 
to the longitudinal structure function 

FffiW = ^^V^ / dy-e-i p+ «-*(P | J-(y | 0) | P) - 4^xF 2(/) (x) (5.20) 

where / represents a quark (q) or anti-quark (q) or both depending on the target | P), 
and where the twist two contribution to the i*2(/) structure function 

X 

The bilocal current operator 



%l d y~ e ~ l2P+v ~ x ( p l J + (v l °) l p )- ( 5 - 21 ) 



j\v i o) = \-. U>(y)rm - mn>(y)} ■ (5.22) 

Note that in the case of quark (anti-quark) contributions, the second (first) term in the 
expression for the bilocal current in Eq. (|5.22[) vanishes, because of a momentum non- 
conserving delta function. We have 

F[^(x) = M 1 + M 2 , (5.23) 
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with 



Mi 



and 



M, 



J_ x 2 (P+) 2 

q 2 2^ 



(P x ) 2 1 x 2 (P H 



dy" e -i p+ ^- (P | V(r)7^(0) | P), 



tfy- e -i p+ ^ * (P | ^(^) 7 +^(0) | P). 



(5.24) 



(5.25) 



(P+) 2 Q 2 2tt 

Mi contains the 'bad' (minus) component of bilocal current operator. This involves the 
constrained field ip~ and one has to use the equation of constraint to eliminate it. 
Using ijj{y~)^~ip{fS) = 2ip~' (y~)tp~ (0) , where 



*l>~(z) = j % I dy-e{z- - y~) .{id 1 - + gA 1 ) + 7 °m ^ + {y ), 



(5.26) 



and e(x ) - 
Mi 

and 



— ±p f dw 6 2^ x we arrive at 

n J 



hl dre 



a ± .[^ ± + gA ± {y)\ + 7°ml [^.[id 1 + gA J 



+ 7°m 



V> + (0) | P), (5.27) 



1 



-x 



dy-e—* p+y X (P | ij + \y-)^j + (0) \ P). 



(5.28) 



Q 2 2tt 

We see that the elimination of the constrained field ip~ brings in interaction (g) de- 
pendent terms in the operator in Mi. In fact, this feature distinguishes higher twist 
contributions from the leading twist, since it directly involves quark-gluon interactions. 
Thus we have obtained an expression for the twist four part of the fermionic component 
of the longitudinal structure function. 

After the establishment of QCD as the underlying theory of strong interactions, the 
twist four part of the quark contributions to the longitudinal structure function has been 
given in the limits of vanishing target transverse momentum and massless quark and in 
the A + = gauge using the QCD factorization method [HI. 9 



F 7 L{q j[x) 



vrQ 2 



dy-e-2 p+y ~ x (P | ij + \y-) 



+t/ 



[id 1 + gA ± {y)}.a ± a ± .[id ± + ^(0)]^ + (0) | P). (5.29) 



Using Eqs. Q5.23 ), ( 5.27 ), and ( |5.28 ), and taking the limit of vanishing target trans- 
verse momentum and massless quark, our result given in Eq. ([5.20 ) reduces to that 
obtained via the QCD factorization method. 
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Note that the expression for FJT^ given in Eq. (|5.20|) appears to violate transverse 
boost invariance since it involves P x . But, we exhibit below with explicit calculations in 
Sec. 5.6 that the P 1 - dependence cancels between Eqs. ( |5.27| ) and ( |5.28| ) so that the full 



^t(q) is indeed boost invariant. 



5.5 Sum Rule for the Twist Four Longitudinal Structure Func- 
tion 



From Eq. ( ggTJ ) it follows that F 2 



(/)» 



-x) 



2(/) 



i X 



and from Eq. ( |5.20| ) we explicitly 



find that F[,£(—x) = —F[^(x). Consider the integral 
dx- 



X 



2 / dx 

Jo X 



00 F kf){ x ) 



/•+ 00 dx 
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{xP + f I 


/— oo X 


[v 


7T J 



(P\J (y\0)\P) 



dy e 

(p^y 



\P+y-x 



(P | J + (y | 0) | P) 



(5.30) 



The part containing J can be written as 
(P+) 2 



2ttQ 2 



dx 



dy -^-(e^ p+ y~)(^)(p | I p) 

oy~ iP + 



dy 







e 2 



P+y- 



)(-^-)(P\m7-^(y-)\P), 



(5.31) 



dy~ v /y iP+' 

where in the first term, we have changed y~ to — y~ and used translational invariance of 
the matrix element. Interchanging the orders of x and y~ integrations and carrying out 
the x integration, we get, from the above expression, 



2zP H 



dy-(P | ^(0) 7 -9+V(0) | P)5( l -P + y-) 



Q 



2 (P | ^(0)7^9+^(0) | P). 



(5.32) 



Similar manipulations can be done for the part containing J and we arrive at |10[ , 



x 



4 

~Q 2 



(P | iip-f d + -01(0) 
(P^ 2 



(P|^ 7 + 9 + 0| { o)|P) 



(5.33) 
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Identifying iip r y~d + ip = the fermionic part of the light-front QCD Hamiltonian 

density and iip^ + d + ip = 0t + , the fermionic part of the light-front QCD longitudinal 
momentum density, (see Eqs. ( |5.36| ) and ( |5.38| ) below), we have: 

jf d* m *' Q2) = ^[{p I «T(0) | p) - p^(P I I P)l (5-34) 

where we have used the fact that the physical structure function vanishes for x > 1. 

The integral of L< - f) is therefore related to the hadron matrix element of the (gauge 
invariant) fermionic part of the light-front Hamiltonian density. This result manifests the 
physical content and the non-perturbative nature of the twist-four part of the longitudinal 
structure function. 

The fermionic operator matrix elements appearing in Eq. ( |5.34j ) change with Q 2 as a 
result of the mixing of quark and gluon operators in QCD under renormalization. Next 
we analyze the operator mixing and derive a new sum rule at the twist four level. 

The symmetric, gauge-invariant energy-momentum tensor in QCD is given by 

9^ = ^i[^D- + j u D^-F^ a F\ a +^{F Xaa ) 2 

(h X D x - m) ip. (5.35) 



The last term vanishes using the equation of motion. Here F^ va = d^A va — d v A^ a + 
gfabcA^A^. Formally, we split the energy momentum tensor into a "fermionic" part 9^ u 
representing the first term in Eq. (|5.35f) and a "gauge bosonic" part 6^ v representing the 
second and third terms in Eq. (|5.35| ). To be consistent with the study of deep inelastic 
structure function formulated in the A + = gauge, we shall work in the same gauge. We 
have, for the fermionic part of the longitudinal momentum density, 

e++ = hp^+d+tfj. (5.36) 

For the fermionic part of the Hamiltonian density, we have 

d+- = iip +j d~ip + + gip +i A-ip + + ^~ f d + <0~. (5.37) 

Using the Dirac equation for the fermion, we find that the sum of the first two terms 
equals the third term in the above equation. Therefore, 

0+- = i^-d+ip = 2iip- 1 d + ip~ (5.38) 

= 2^[a ± .(id ± +gA x ) + 7 Q m] — [^.(i^ + gA x ) +7°mU+. (5.39) 

The gauge boson part of the Hamiltonian density is given by 
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= —F +Xa F Xa + \g + -(F Xaa y = \{d + A' a ) 2 + \f^F«. 

= (d*Ai) 2 + 2gr bc A\A{d l A{ + A\A{A\A{ 

+ 2gd i A i a (i-) (/ a6c ^9 + ^ c + 2(^ + ) t T> + ) 

+ <? 2 (^) (/ abc ^9 + ^ c + 2ty + ) t TV) (ip) (/ ade ^9 + ^ e 
+ 2(^+) t T> + ) 



(5.40) 



where we have used the equation of constraint for the gauge field. 
We define the twist four longitudinal gluon structure function 



1 xP^ 



Q 2 2tt 



dy e 



-^P+y-x 



(P | (-)F +Aa (y-)F Aa (0) + -g + -F^ a (y-)F Xaa (0) + (y~ <-> 0) | P) 



(P^l 

(p+y 



(P | F +Xa (y-)F+ a (0) + (y <-> 0) | P) 



(5.41) 



Then we have 

r 1 dx 
Jo x 



771T=4 , pT=4 



1 ^ z?r =4 
X 



2 

Q~ 2 

(p^r 



(P | e + -(o) I P) - 

(P I 6 ++ (0) I P) 



(5.42) 



(P+) 2 

Neglect of M 2 in the expression of Fl (see Appendix B) will lead to (P -1 ) 2 + M 2 instead 
of (P -1 ) 2 in the right hand side of the above equation and this would spoil the correct 
sum rule. 
We have 

(P | e + -(0) | P) = 2P+P- = 2(M 2 + (P x ) 2 ) and (P | # ++ (0) | P) = 2(P+) 2 , (5.43) 

where M is the invariant mass of the hadron. Thus we arrive at the new sum rule for the 
twist four part of the longitudinal structure function 



L 



dx F[=4 = 4 M^ 



(5.44) 



To our knowledge, this is the first sum rule at the twist four level of deep inelastic 
scattering or for QCD in general. 
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5.6 Dressed Parton Calculations 
5.6.1 Dressed Quark with Zero Mass 



Next, we investigate the implications of Eq. Q5.34Q for quadratic divergences in FJl q y 
For simplicity, we select a dressed quark target and evaluate the structure functions to 
order g 2 . That is, we take the state | P) to be a dressed quark consisting of bare states 
of a quark and a quark plus a gluon: 

\P,a) = <f)itf(P,a) | 0) 

d k-^ r dJ\i>2 d k 



E / 7^'\ / ^p= p{2irfpn\p -h- h) 



<p 2 (P,a | k 1 ,a 1 ;k 2 ,X 2 )b^(k 1 ,a 1 )a\k 2 ,X 2 ) I 0). (5.45) 

Here we have truncated the expansion at the quark-gluon level. We introduce the 
Jacobi momenta, (xi,K,j-), where, 

k+ = Xi P + , kf = + %iP L (5.46) 

so that 

5> = 1, E^ = 0- M7) 
The amplitudes (fii,4> 2 are related to the corresponding boost invariant amplitudes as 

0! = Vi, Vp+M*?, k i) = M*, *t). (5.48) 

Using the notation, and Eq. ( |5.47| ) we have 



[M 2 



m 2 + ( K X)2 ( K X)2 1 



1-x 



— T a ^_ -2 <T - zmtr 1 ef*. 5.49 

2(2vr) 3 v 7 !^ 1-^ * 

Here M and m are the masses of the dressed and bare quarks respectively. The state is 
normalized to one, and from the normalization condition, we obtain 



| iM a = 1 - £.C, ^ dx\±^ln% (5.50) 
2n Je 1 — x fi z 

within order a s . Here e is a small cutoff on x. 
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Mi = I dy-e-i p+ y-*(P \ V> + V) 



First we work in the massless limits for the dressed and bare quarks. 
1 

^Q 2 

[id 1 - + gA ± (y)].oc L ac L .[id L + ^ ± (0)]^ + (0) | P>, (5.51) 

and 

M * = - A ^r xF m( x )- ( 5 - 52 ) 

As mentioned before, in this case, it is sufficient to keep only the first term of the bilocal 
operator. Note that the matrix elements appearing in Eqs. (|5.27| ), fl5.28Q , ( |5.39|) , and 
( |5.40| ) involve products of operators. We treat them as normal ordered which is sufficient 
for our purposes here. However, in doing so, the terms we drop are quadratically divergent 
and they will affect only the counterterm structure. 

We express the operators in terms of two component fermion field, £, and use the 
Fock space expansion of the fields £ and A 1 . First we evaluate the contribution M 2 given 
in Eq. ( |5.52| ). The operator in this case involves the 'good' (plus) component of the 
bilocal current which is independent of interaction and the matrix element receives only 
diagonal contribution from the single particle and two particle Fock space sectors. M 2 is 
proportional to F 2 . A detailed derivation of the structure function F 2 for a dressed quark 
state is given in We obtain 



M 2 = -AC f [ —fx 2 \5(l-x) + -^\nA 2 [--^--5(l-x) dy-^-]] (5.53) 
U z 1 In 1 — x J L — y J 



Q 1 

where we have cutoff the transverse momentum integral at A. Note that the result in Eq. 
( |5.53| ) violates transverse boost invariance. In the above expression, we have used Eqs. 
( jT5g ) and (Egg) . Here, C f = ^§=1 for SU(N). 

Mi involves the 'bad' (minus) component of bilocal current and the operator depends 
upon interaction. We have, from Eq. ( |5.51| ), 

-Mi = / dy-e~^ p+ y-%P | V^VX^WO) | P ) 

+ 9-7^ I dy-e-^ p+y ~ x (P | tfj + \y-)id ± .a ± a ± .A ± (0)^ + (0) | P> 

TlQ J 

+ 9^ I dy-e'^ p+ y-%P | ^ + \y~)a ± .A ± (y)zd ± .a ± ^ + (0) | P) 



+ 9 2 ^Q2-j dy-e-^v-^P \ ^ + \y-)A^{y)A^{U) \ P) (5.54) 
= M\ + M\ + M\ + Mf ■ (5.55) 
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Here, the operator in M.\ only does not depend explicitly upon interaction and will have 
diagonal contributions. Ai\ and M.\ explicitly depend upon g and the contributions will 
come due to the interference between the one particle and two particle sectors. This 
type of off-diagonal contributions are absent in the structure function F 2 and they involve 
quark-gluon dynamics. Since the operators in Eq. ( |5.51| ) are taken to be normal ordered, 
the contribution of A4f vanishes to order g 2 . 

Explicit calculation leads to the diagonal Fock basis contributions 



(.Mi)** = M\ = 4C f (^lx 2 \d(l -x) + ^lnA 2 [i±^ - 5(1 - x) f dy 1 -^- 

Q z 1 2tt 1 — x J I — y 

g 2 A 2 l+x 2 

+ ^2^^i— 7- (5 ' 56) 

The first term here explicitly cancels the term Ai 2 given in Eq. ( |5.53| ), thus the P x 
dependence goes away. Here, we have used the expression of ip 2 given by Eq. ( |5.49| ) and 
Eq. (|5T50|) . 

Off-diagonal contributions will come from Ai\ and Ail and these matrix elements 
will contain ip^ipi and its Hermitian conjugate. The following results are obtained after 
substituting ip 2 : 

(M^^ = Ai\ + M\ = -Cf^^Y^' (5 ' 57) 



Adding all the contributions, we have 

Fli*{x) = -C f j- 2 — {! + ,). (O.os 



As anticipated from power counting, we have generated quadratic divergences for Fl 
in the bare theory. This quadratic divergence arises from k x integration. Now we are 
faced with two issues: (a) What is the principle for adding counterterms? (b) What 
determines the finite parts of the counterterms? 

Further, since FJT$ is directly related to the physical longitudinal cross section, we 
expect FJT^ to be positive definite (see, for example, Ref. ||). From our answers we see 
that we get a negative answer which is free from end point singularities but is quadratically 
divergent. Note that it is easy to show from our expressions in Eqs. (|5.27|) and ( |5.28| ) for 
a free quark of mass m, F£^ = 4^-5(1 — 2), a well-known result. 

Now, the question is whether we can invoke the sum rule to get insight into the answers 
we obtained in the bare theory and also determine the counterterm structures to be added 

52 



for renormalization. To investigate whether our sum rule can throw some light on this 
problem first we have to check the validity of the sum rule to order g 2 in perturbative 
QCD. Towards this goal we have to evaluate the expectation value of the H — component 

of the energy-momentum tensor which we consider in the following. 

FTrUx) 

Now we show that the sum rule relating the integral of — ^ — to the fermionic part 
of the light-front Hamiltonian density helps us to understand the results we obtained 
for Fl for a dressed quark in the bare theory and indicates how to add counterterms to 
renormalize the theory. 

A straightforward evaluation leads to 

(P I 0r(0) | P) nondiag = -Cf^A 2 / ^ T ^, (5.59) 



(p I e+-(o) | P) diag - |^£(P | 6f+(0) | P) dia9 = c f \^j (5-60) 

Adding the diagonal and off-diagonal contributions from the fermionic part of the 
Hamiltonian density, we arrive at 

2 r ( P^ 2 i n 2 A 2 r dr 

Q~ 2 [(P l CP) l P ) ~ l ^ + (0) l ^)] = ~ C fh^ 2 J + ^ (5 - 6l) 

In all the x integrations we have taken a small x cutoff e. Comparison of Eqs. (|5.58|) 
and (|5.61| ) immediately shows that the relation given in Eq. (|5.34j) has been successfully 
tested to order g 2 in perturbative QCD in the bare theory. 

It is now clear that the quadratic divergences we have generated in FJT^ are directly 
related to the fermion mass shift due to the fermionic part of the light-front Hamiltonian 
density in perturbation theory. Since we want the renormalized fermion mass to be zero, 
we need to add counterterms to the Hamiltonian density. This precisely dictates the 
counterterms to be added to F£Z\ in order to renormalize via the relation given in Eq. 

To complete the discussion of the fermion mass shift, we now consider the contributions 
from the gluonic part of the Hamiltonian density. 
Take the off-diagonal contribution arising from 

M 3 = -(P | 4#^(«9M w )^ +t T> + | P). (5.62) 
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A straightforward evaluation leads to 



M, = -d 



9 tA 2 



2vr 2 



dx r 2x 2 x 
+ 



x 



1-x) 



1 — X 



(5.63) 



Thus from Eqs. ( |5.59|) and ( |5.63j ) we have 



k + -(o)+e-(°)l ip> 



nondiag 



-c 



A 2 



2vr 2 



da; 1 + x 2 
s (1 — a;) 5 



(5.64) 



Next let us consider diagonal contributions. From the gauge boson part of the Hamil- 
tonian density, we have, 



<P | 0+-(O) | P> 



diag 



■ (Pie + (0)|P) 



diag 



1 47T 2 



dx 1 + x 2 
1 — a; 1 — x 



(5.65) 



Adding the diagonal contributions from the fermion and gauge boson parts, i.e., Eqs. 
Q5.60Q and (|5.65|) we arrive at 



(P I + "(O) | P) diag - ^-(P | 6++(0) | P) diag = C f £- 2 A*J 



dx 1 + a; 2 



(P 



x 



(5.66) 



Thus the total contribution (quark and gluons) to the expectation value of the Hamilto- 
nian density is given by 



(P | 0+-(O) | P) 



(P 



(P | 9 ++ {0) | P) = -C 



4vr 2 



A 



dx 1 + x 2 
x (1 — x) 2 



(5.67) 



This result is directly related to the mass shift of the fermion to order g 2 in light-front 
perturbation theory. In the x + ordered perturbation theory, upto order g 2 , contribution to 
the mass shift of the quark comes from three diagrams, (a) one gluon exchange interaction, 
(b) instantaneous fermion interaction, (c) instantaneous gluon interaction. It can be 
shown that, the contribution to the fermion mass shift due to (a) in the limit of zero bare 
quark mass is, (see Eq. (4.10) in Ref. [|T2]j ), 

dx 1 + x 2 



2P 



1 ~ 2 
Cf 



9 A2 



47T Z 



x (1 — x) 



(5.68) 



Note that in the massless limit we encountered only quadratic divergences in the twist 
four part of the longitudinal structure function. The case of a massive quark is discussed 
in the next section. 
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5.6.2 Dressed Quark with Non-zero Mass 

We select the target to be a massive dressed quark and evaluate the structure functions 
to order g 2 . That is, we take the state | P) to be a dressed quark consisting of bare states 
of a quark and a quark plus a gluon, as before 



P,a) = ^(P,a)\0) 

/ 



^2 a fc 2 I 2 (2kYP+5\P -h-k 2 



4> 2 {P, cr I ki, a\\ k 2 , A 2 )& t (fci, ai)a\k 2 , A 2 ) | 0) 



(5.69) 



In the last section, we have shown that the twist four longitudinal structure function 
has quadratic divergences in perturbation theory. In this section, we show that for a mas- 
sive quark, in addition to quadratic divergences, logarithmic divergences are generated. 
The state is normalized to one. As before, we introduce the Jacobi momenta. We have 



F L = M.i + M 2 



(5.70) 



where 



Mi 



vrQ 2 



+t /-„.-> 



a 1 . [id 1 - + gA ± (y)} + 7 °m] [a 1 . [id 1 + ^(0)] + 7 °m]^+(0) 
+ h.c. | P), 



(5.71) 



and 



M 2 = -A^-xF 2{q) (x) 



(5.72) 



where | P) now has a mass M and m is the bare quark mass. 

In the case of quark contributions, the second term in the expression for the bilocal 
current in Eq. ( [5.9D vanishes. First we evaluate the contribution Ai 2 given in Eq. ( [5.721 ). 
We obtain 

M 2 = - A^^x 2 \d(l - x) + -f^Cf (i \/ 2 /v 



where C 



f 



N 2 ~l 
2N 



Q 2 

- 5(1 - x) I dyd 2 k 



for SU(N). 



l+x 2 u2 
l-x 



k'j + (1 — x) 3 m 2 



7^1 + (1 - y) 



[m 2 (l - x) 2 + k\} 2 

3 m 2 



[m 2 (l-y) 2 + k 2 ± } 2 



(5.73) 



55 



Here we have presented the result without working out the transverse integration to 
maintain a greater degree of transparency. 

It is to be mentioned that, to calculate all the matrix elements involved in Fl, one has 
to use the expression of ip2 given in Eq. ( 5.49 ). It can be easily seen from Eq. ( 5.49|) that 
the term which in linear in quark mass m causes helicity flip of the quark. Linear mass 
terms do not contribute in F^. However, in the two-particle sector, the diagonal part of 
the matrix element involves ^2^2 an d contribution comes from the terms quadratic in m. 

The contribution from Mi is split into four parts with additional contributions coming 
from quark mass terms and can be written as follows. 

1 r . _ _id+.,-»,„ , _ N 



Mi 



ttQ 2 



+ 9 



+ 9 



dy-e—* P y ~ x (P I ^ \y~){ - (c^) 2 + m 2 )^+(0) | P) 



dy-e-$ p+y ~ x (P I ^(y^iid^.a 1 + j°m)a ± .A ± (0)i/; + (0) | P) 



nQ 2 
1 

nQ 2 

1 



dy-e—* p+y ~ x (P I tP +, (y-)a ± .A ± (y)(id ± .a ± + 7 u m)^+(0) | P) 



-.+t 



dy~e-i p+y X (P I 4> +, (y-)A ± (y).A ± (Q)4> + (0) \ P) 



,+t 



vrQ 2 

Ml + M\ + Mi + M( 



(5.74) 
(5.75) 



Since the operators in Eq. ( p.71| ) are taken to be normal ordered, the contribution of Mf 
vanishes to order g 2 . 

In this case also, M\ gives diagonal contribution and M\ and M\ give off-diagonal 
contribution. Explicit calculation leads to the diagonal Fock basis contributions 



(M 



1 ) diag 



Ml 



Q 2 



-X 



g- , r 2 ^k'{ + (1 - xfm 2 
8^ Cf (J 6 k± [m 2 (l — x) 2 + kj] 2 



— 5(1 — x) I dyd 2 k 



Am 2 . 
+—6{1 -x) 



3 m 2 



T^fcj + (1 - y) 
[m 2 (l -y) 2 + k\} 2 



Q 2 

Q 2 8tt 3 



g 1 



i-c f 

d 2 k ± (k]_ + m 2 



dyd 2 k\ 



7^1 + (l-y) 



3 m 2 



8tt 3 7 " M \m 2 (l-y) 2 + k\} 2 
l -^k\ + (1 - xfm 2 
\m 2 (l - x) 2 + k 2 ± } 2 ' 



(5.76) 



The first term here explicitly cancels the term M2 given in Eq. ( |5.73| ). 
Off-diagonal contributions 

C> Ar,i ^ /" J -. J 2,_ m 2 (l-y) 

Q 2 7T 3 



{Ml)nondiag = M.\ + M\ 



5(1 - x) / dyd 2 k 



[m 2 (l - y) 2 + k\] 
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The linear in m terms in Ai\ and Ai\ receives contributions from the mass dependent 
term in ip2. As a result, in the final expressions, these are quadratic in m. 
Adding all the contributions, we have 

4m 2 4C f a 2 r f , , , ^k 2 ± + (1 - xfm 2 

^2 + (1 _ yfm 2 



-5(1 -x)m 2 J dyd 2 k : 1_ " 



[m 2 (l-?/) 2 + & 2 ] 2 



2 



cTfc 



g 2 7r 3 Ly "(i-x)[m 2 (i-x) 2 + £; 



2 ' 



, m 2 (l-y) 



-8{l-x)J dyd k ± [m2(1 _ y)2 + kl] \. (5.78) 

Here we have used M = m, since the difference that it entails is higher order in the 
coupling. Note that we are getting back the free quark answer once we switch off the 
interaction. Also, the dressed mass-less quark answer can be easily regenerated by putting 
M = m = 0. Note that the fc^-integration now produces logarithmic divergences with 
the expected quadratic one. 

To check the sum rule explicitly, we evaluate the right hand side of Eq. ( |5.34j ) next. 
A straightforward evaluation leads to 

(P | «T(0) | P)^ 9 = -C,£- S J fed V 1 + "I 1 - Xf _ 1 (5.79) 



(p | e*~(o) | p) dm „ - jj^(p I «, ++ (o) I p>*« s 



Adding the diagonal and off-diagonal contributions from the fermionic part of the Hamil- 
tonian density and multiplying it by one obtains the right hand side of the sum rule. 



Comparing it with the integral of L , where Fl is given in Eq. ( 5.78 ), one easily sees 
that the sum rule is verified. 

To see the connection of Fl with the fermionic mass shift, we calculate the contribution 
of the gluonic part of the energy momentum tensor 9 + ~ to the sum rule for the total Fl- 
Explicit calculation gives 

(p I er (o) i pu*. = -c,i-Ji ^\^ [m2{1 _\ )2 + kl] (5.8D 
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T=4 



(p i &i~ (o) i p) dia9 - |^(p i ^ 9 ++ (o) i p) diag 



q 2 r k 2 ^k 2 +(l-x) 3 m 2 
2C f^ / dxd 2 k ± /i 1 . ^ ^ \ 9 ^— . (5.82) 



Thus, we get 



(P|0+-(O) + 0+-(O) \P) nmdiag = 

q 2 f 9 (i±^ifc 2 + (1 - x) 3 m 2 1 

-C f ^r / dxrf 2 fc x ^ —. - - TT , (5.83) 

f 2ir 3 J x(l-x) [m 2 (l - x) 2 + klY K ' 



(p i or (o) + e+-(o) i p) dia9 - |^(p i r + + ^ ++ (o) i p) diag 



9 



2 



^^kl + (1 - xfm 2 1 



8^] x(l — x) [m 2 (l — x) 2 + /qj 

Adding diagonal and off-diagonal contributions, we get, 

(^l^(0)l^)-7S^I^ + (0)l^) = 



,2 ^g^i + (1 - a;) 3 m 2 



- C/^-g / tfcrf 2 fc ± '~ x - v " " y - % 2 ■ fr 2T ( 5 - 85 ) 

An 6 J x{l — x) [m 2 {l — x) 2 + k 2 ± \ 

Note that this result is connected to the full fermion mass shift 5p{ in second order 



perturbation theory. We have (see Eq. (4.10)) in Ref. [12], in the massive case 



1 q 2 r , k 2 + (1 - xfm 2 1 

2P+ J An 6 J x(l — x) [m 2 (l — x) 2 + k^_\ 

5.6.3 Dressed Gluon 

In this section we check the sum rule explicitly for a dressed gluon target. We consider 
the gluon to be composed of a bare gluon and a quark anti-quark pair. 

\P,a) = ia t(P,A) |0) 



dk-^ d /l-j^ f dk<2 d k 



+ y, , - / -2 J 2 (2tt) 3 P+5 3 (P -h- k 2 ) 

<^2 J ^2(27r) 3 kf J ^2{2Tifkt 

4> 2 {P,o | ki,ui; k 2 , a 2 )b ] {k 1 ,a 1 )S{k 2 , a 2 ) | 0). (5.87) 

0i is the probability amplitude of finding one bare gluon and <j) 2 is the probability 
amplitude of finding one quark- ant iquark pair in the dressed gluon. We introduce the 
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Jacobi momenta x, k x . The boost invariant amplitude ip2 is written in terms of ipi from 
the light-front bound state equation (see chapter 3): 



9 x(l - x) f 

2 -^cr' 



From this, we get 



[o ■ k ) a • ( <t • /-r 



.r 



1 — x 



spins 



(2tt) s 



[x 2 + (l-x) 2 ]|^i I 2 - 



(5.89) 



The target gluon and the bare quark and anti-quark masses are taken to be zero. Note 
that, to the order g 2 , there will be a contribution from the two-gluon Fock sector due to 
the non-abelian nature of the gauge coupling. For simplicity, we exclude that contribution. 
It is easy to incorporate that contribution by trivially extending our calculation presented 
here. 

Fl can be written in terms of M.\ and M.2 given in Eqs. ( |5.71| - |5TT^ ), where | P) now 
stands for the dressed gluon represented by Eq. ( |5.87| ). In this case, contribution comes 
from both components of bilocal current. It can be easily seen that Fl for a free gluon is 
zero. Explicit calculation gives 



-4 



■pdressed— gluon 

^-^Tfix 2 + (1 - *)>A 2 . 



(5.90) 



Here Tf = | and Nf is the number of flavors. 

M i is again divided into four parts as in Eq. (|5.75| ) and explicit calculation in this 
case gives the following. 

x 2 (P x ) 2 g 2 



(diag) 



Mi (a 



(a) 



Q 2 ^N f T f [x 2 + {l-xY]\nV 
+ ^N f T f [ X 2 +(l-xn 



(5.91) 



Ml(off-diag) = Ml(b)+Ml( c ) 

A 2 q 2 

= -^JWti-z). 

In the above results, we have used Eq. ( |5.88| ) and Eq. ( p.89|) . Thus, we get 



A 2 2 

Fl = T^NfTjUx 2 + (1 - x) 2 - 2(1 - x)}. 



Q 



71 



(5.92) 



(5.93) 
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On the other hand, we get 



P I I P ><^ - 7pT^( P I e \ ; + (0) I i 3 ) 



diag 



and 

2 2 / 1 \ 2 

(P I C (0) I P)off- d ia 9 = -^N f T f ^- 2 I dx f +^~*> ]. (5.95) 
Adding diagonal and off-diagonal contributions, we get 

(P I C (°) I P ) - l^fJ P I C(0) I P) = ~ A2N f T S I ^ 1 (5 ' 96) 
Note that this result is connected to the gluonic mass shift 5q 2 due to pair production, 
since the contribution from the gluonic part of the energy-momentum tensor in this 
case vanishes. In the massless limit, we have (see Eq. (4.40) in Ref. (||12||), gluonic mass 
shift due to pair production is 

5q 2 = l —k 2 N f T f ^- I dx f^ 1 '^ 2 }. (5.97) 

H2 2P+ 1 Un 2 J 1 x(l -x) J V ; 

From Eq. ( |5.93| ) we compute J dx^-. Since x-integration is from to 1, it can be written 

in the following form. 

Comparing Eq. ( |5.96|) and Eq. ( |5.98| ) , one explicitly verifies the sum rule for a dressed 
gluon target. 

As we have emphasized, in the bare theory, the twist four longitudinal structure func- 
tion is afflicted with divergences. We have to add counterterms to carry out the renor- 
malization procedure so that we have physical answers. The sum rule for the bare theory 
clearly shows that the quadratic divergences generated are directly related to the gluon 
mass shift in second order light-front perturbation theory arising from an intermediate 
quark anti-quark pair. As we mentioned earlier, in light-front Hamiltonian perturbation 
theory one cannot automatically generate a massless gluon by dimensional regularization. 
In order to ensure a massless gluon in second order perturbation theory, we have to add 
the negative of the shift as a counterterm. After adding the counterterm, the gluon mass 
shift in second order perturbation theory is zero and the twist four longitudinal structure 
function for a massless gluon becomes zero. Thus, after renormalization, the sum rule is 
satisfied, with a trivial (i.e., zero) gluon longitudinal structure function. 
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5.7 Partition of the Hadron Mass in QCD 



As is well-known, experiments that measure the twist-two part of the F2 structure 
function yield information on the fraction of longitudinal momenta carried by the charged 
parton constituents of the hadron (quarks and anti-quarks). The sum rule we have derived 
yields other useful information about the hadron structure. Namely, our sum rule shows 
that experiments to measure the twist four part of the longitudinal structure function will 
directly reveal the fraction of the hadron mass carried by charged parton components of 
the hadron. This sum rule is richer in content since it involves dynamics. The light-front 
Hamiltonian provides theoretical insight into this fraction as follows. 

According to our analysis, the twist four part of the longitudinal structure function is 
directly related to the fermionic part of the light-front QCD Hamiltonian density 0j~~ in 
the gauge = 0. explicitly we have 

1 



01 



2ip 



+t 



a 



id+ 



a. (id- 1 + gA ) + Ym i> 



(id 1 - + gA 1 ) + 7°m 
Thus we have the fermion kinetic energy contribution given by 

and the interaction dependent part given by 



(5.99) 



(5.100) 



^ q(int) 



2# 



+t 



a x .A J 



1 {a 1 . id 1 + 7°m) + (a^.id 1 + 7°m)ia ± .A J 



id+ 



id+ 



4>- 



+ 2g 2 ^ +j a ± .A ± a ± .A ± ^ + . 



(5.101) 



Note that the fermion kinetic energy constitutes only a part of the total contribution from 
fermions. Any theoretical estimate of the fermionic part of the longitudinal structure 
function necessarily has to involve off-diagonal contributions from Fock states differing in 
the number of gluons by one and two. 

It is important to emphasize the difference between equal time and light-front Hamil- 
tonians in the context of our calculations. The equal-time Hamiltonian contains the scalar 
density term (tpip) accompanying the quark mass m. In contrast, the quark mass appears 
quadratically in the free part of the light-front Hamiltonian. Recently the question of 



the partition of hadron masses in QCD has been addressed by Ji [13] in the context of 
the equal-time Hamiltonian and in terms of twist-two and twist-three observables. In 
his analysis, the extraction of the fraction of the hadron mass carried by the fermion 



61 



constituents is not straightforward because of the presence of the scalar density term. 
The hadron expectation value of the strange quark scalar density remains unknown (ex- 
perimentally). Our analysis, however, shows that the twist four longitudinal structure 
function, once extracted experimentally, directly yields the fraction of the hadron mass 
carried by fermionic constituents. 

5.8 Discussions 

Our results indicate that the experiments to measure the twist four longitudinal struc- 
ture function reveal the fraction of the hadron mass carried by the charged parton compo- 
nents. Thus these experiments play a complementary role to the longitudinal momentum 
and helicity distribution information obtained at the twist two level. It is of interest 
to investigate the feasibility of the direct measurement of the twist four gluon structure 
function in high energy experiments. Recent work of Qiu, Sterman and collaborators have 
shown that semi-inclusive single jet production in deep inelastic scattering and direct 



photon production in hadron nucleus scattering [[DJ provide direct measurement of twist 
four gluon matrix elements. 

We also note that in the pre-QCD era, there have been discussions about a possible 
S(x) function contribution to the longitudinal structure function which may appear to in- 
validate the sum rule derived ignoring such subtleties. In two-dimensional QCD Burkardt 
has shown [IB| that ^ has a delta function contribution which comes from surface terms 



and he has discussed implications of this for the sum rule for Obviously, — will not 
be affected by such a singular contribution and in the next chapter, we show explicitly 
that the sum rule is verified in two-dimensional QCD by virtue of the 't Hooft equation. 

In the next chapter, we present a non-perturbative calculation of the longitudinal 
structure function for a meson-like bound state in 1 + 1 dimensional QCD and for a 
positronium-like target in 3 + 1 dimensional QED in weak coupling limit. 
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CHAPTER 6 



Twist Four for Light-Front Bound 

States 

Higher twist or power suppressed contributions to the DIS structure functions involve 
non-trivial non-perturbative information on the structure of hadrons. As we mentioned 
before, the most important point of our approach based on light-front Hamiltonian QCD 
is that both perturbative and non-perturbative calculations are possible within the same 
framework [|IJ. We have calculated the twist four part of the longitudinal structure func- 
tion Fl in the previous chapter for a dressed quark and a dressed gluon state in perturba- 
tion theory. When the target is a bound state in QCD like a nucleon or a meson, one can 
use the same technique of Fock space expansion of the target state, but then has to know 
the bound state wave function. The analysis of QCD bound states in 3 + 1 dimension 
in light-front Hamiltonian theory is highly complicated and it requires the recently de- 
veloped similarity renormalization technique 0. The spectra of the heavy quark bound 
states like charmonium and bottomonium have been investigated using this technique 
H]. However, the analysis of the structure functions requires knowledge about the bound 
state wave functions, the analytic form of which is not obtained so far. This is because, 
the similarity renormalization technique generates an additional confining interaction in 
the effective Hamiltonian in 0(g 2 ) which makes it impossible to solve the effective Hamil- 
tonian analytically even in the leading order in bound state Hamiltonian perturbation 
theory |§. 

In this chapter, we investigate the twist four part of Fl for a meson in 1 + 1 dimensional 
QCD. We explicitly show that the sum rule which we derived in the previous chapter is 
obeyed in this case by virtue of t'Hooft equation. In order to understand the calculational 
procedure in 3 + 1 dimension, we perform a simpler but interesting analysis. We calculate 
the twist-four part of the longitudinal structure function for a positronium-like bound 
state in light-front QED in the weak coupling limit. In this limit, QCD results are not 
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expected to differ much from the QED results. As a result, this analysis is important 
since it tests and illustrates the approach in QCD. The advantage here is that, in the 
weak coupling limit, the bound state equation can be solved analytically and the wave 
function is known. This allows an analytic understanding of the problem. We refer to our 
original works, |]4]. |5[. 

The plan of this chapter is as follows. In chapter 6.1, we calculate twist four part of 
Fl for a meson in 1 + 1 dimensional QCD and explicitly verify the sum rule. In section 
6.2, we calculate the twist four longitudinal structure function for a positronium-like state 
in light-front QED in the weak coupling limit and verify sum rule for F[ =4 in this case. 
We also show that it reduces to a relation connecting the kinetic and potential energies 
to the binding energy of positronium. In section 6.3, we analytically calculate F[ =4 using 
the wave function for positronium. Discussions and conclusions are given in section 6.4. 
The bound state equation for positronium in light front QED in the weak coupling limit 
is derived in Appendix C. 

6.1 1 + 1 Dimensional QCD: Explicit Calculations 

In the previous chapter, we have derived a sum rule for the twist four part of the longi- 
tudinal structure function Fl- We have shown that the sum rule is obeyed in perturbative 
QCD for dressed quark and gluon targets. In this section, we turn to two-dimensional 
QCD in order to test the sum rule in a non-perturbative context. In 1+1 dimensions, in 
A + = gauge, we have, 

L - pT ^ {x) = qi {p 1 K~(°) + ^ (0) 1 1 p) ' (6 - 1} 



k x ' Q 2 



with 9q~ = 2m 2 tp + ^ ^+ip + is the fermionic part of the light-front QCD Hamiltonian 
density in 1 + 1 dimension and = — 4g 2 ip + ^T a tp + ■tp + ^T a tp + is the gluonic part. We 
consider the standard one pair (qq) approximation to the meson ground state. Explicit 
evaluations show that 

f\M = J_(_ )G / f dx f m tift^l. (6 .3, 
Jo x Q 2 it Jo Jo (x - y) 2 

Here ip(x) is the ground state wave function for the meson. Thus [|J 

F l ( x ) = 7^ dx ^ ( x ) —R \^\ x ) ~ C f — / d V—( V2 • ( 6 - 4 ) 



ox Q 2 Jo Lx(l — x) it Jo (x — y) 2 
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The bound state equation obeyed by the ground state wave function ip(x) for the meson 
in the truncated space of a quark-antiquark pair is the t'Hooft equation: 

MV(x) = -j^-Mx) - C f 9 - [ (6.5) 

together with the normalization condition Jq 1 dxtjj 2 (x) = 1, we easily verify that the twist 
four longitudinal structure function of the meson obeys the sum rule 

/;^ r ^^l^(0)|P)^f. (6.6) 

In the same model, the contribution to the twist two structure function from the 
fermionic constituents is given by 

F 2{q) (x) = (x + l- x)r(x)ip(x) = i/j*(x)ip(x). (6.7) 

Note that, since there are no partonic gluons or sea in this model, the longitudinal mo- 
mentum of the meson is carried entirely by the valence quark and anti-quark. Thus the 
momentum sum rule is saturated by the fermionic part of the longitudinal momentum 
density. On the other hand, light-front energy density is shared between fermionic and 
gauge bosonic parts and as a consequence the fermions carry only a fraction of the hadron 
mass. This seemingly paradoxical situation further illuminates the difference between the 
physical content of the F 2 and F[ =4 sum rules. 

To get a quantitative picture, next, we explicitly calculate the structure functions 
F 2 ( q )(x) and FL{q ^ for the ground state meson in two-dimensional QCD. We have pa- 
rameterized the ground state wave function as 

tp{x) = A/V(l -x) s . (6.8) 

The factor Af is determined from the normalization condition Jq 1 dxtp*(x)ijj(x) = 1 and 
s is the variational parameter. We have 

The first term in the left hand side is the kinetic energy of the fermions and the other is 
the interaction. Substituting the above wave function and also using the normalization 
condition, we get 

2 (3(2s,2s) is(f3(s,s + l)f-f3(2s,2s)] 

0(2, + i,28 + i) ~ ~ Cf WsTTjsTT) (6 - 10) 
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(m=5, s=4.96) 



(m=l, s=.70) 




Figure 6.1: Fermionic contributions to the structure functions F 2 (x) and — ^— for the 
ground state meson in the 't Hooft model for two different values of m, the quark mass, 
(a): m = 5, s = 4.96. (b) m — 1, s — .70. The parameter s appearing in the wave 

2 

function is determined by a variational calculation. We have set Cf^- = 1. 



where /3(n,m) = J dxx n ~~ 1 (l — x) m_1 . 

We determined the value of s variationally by minimizing M 2 for given values of m 2 
and g 2 . The resulting structure functions are presented in Fig. 1 for two different values 
of m 2 . 

Since both the quark and anti-quark have equal mass in the model, both structure 
functions are symmetric about x = ~. When the fermions are heavy (Fig. 1(a)), the 
system is essentially non-relativistic and the structure functions are significant only near 
the region x = ~. When the fermions become lighter (Fig 1(b)), contribution to the struc- 
ture function from the end-point regions become significant indicating substantial high 
momentum components in the ground state wave function. Note that — ^ measures the 
fermion kinetic energy (in light-front coordinates). The exponent s in the wave function 
is a function of the fermion mass and s decreases as m decreases. In the massless limit, s 
vanishes so that the wave function for the ground state becomes ijj(x) = 8(x)8(l — x). 
This results in a flat F 2 structure function. However, because of the presence of m 2 , FJT^ 
vanishes. Because of an exact cancellation between the self-energy and gluon exchange 
contributions, the gluonic part of the F[ =4 also vanishes. Thus the sum rule is satis- 
fied exactly since, in the zero quark mass limit, the ground state meson is massless in 
two-dimensional QCD. 
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In the next section, we investigate the twist four longitudinal structure function for a 
positronium in 3 + 1 dimension in light-front QED. 



6.2 Twist Four Longitudinal Structure Function for a Positro- 
nium in 3 + 1 Dimensional QED 

We consider a positronium like bound state | P) given by 
I % -v P dk^ d k^ r dk 2 d k 2 

MP I h, <ri\ k 2 , a 2 )^2((2nfP+5 3 (P -k x - k 2 )tf(h, a,)d\k 2 , a 2 ) | 0) 
^ — ^ r dk-^ d k-^ r dk^ d k^ f dk^ d k^ 
+ \h,J y/2(2ir) 3 kt ' p{2nfk+ ' v /2(27r)3fc+ 

C3 



MP I fci, <ri\ h, a 2 ; k 3 , A 3 ) p(2<KfP+5-\P - h - k 2 - k 3 ) 

6t(fc 1 , fl r 1 )rft(fc 2 , r 2 )at(fc3,A 3 ) |0). (6.11) 

Here <p 2 is the probability amplitude to find an electron and positron in the positro- 
nium, 03 is the probability amplitude to find an electron, positron and a photon in the 
positronium. We restrict ourselves to the three particle sector. 

We calculate F£ =A in the weak coupling limit in light-front QED using light front gauge, 
A + = 0. We truncate the Fock state upto the three particle sector. This approximation 
includes the effect of dynamical photon. It can be shown that, for a weak coupling theory 
the results are the same as obtained using a non-relativistic approximation. However, the 
entire calculation is fully relativistic and exact in the leading order in a [[/]]. It can be 
shown that the recently developed similarity renormalization scheme for light-front QED 
gives the same results in the leading order in bound state perturbation theory as in the 
weak coupling treatment and particle number truncation ||. 

The twist-4 part of the fermionic component of the longitudinal structure function is 
given by 

F[ff } (x)=M 1 + M 2 , (6.12) 
Mi = -q 2 A^ J dy~ e^ p+ y- x (P\^(y-)j-m-mi^(y~)\P), (6-13) 
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and 



Mo 



(p 1 ) 2 i x 2 (p+y 



dy e 



\P+y x 



(P+) 2 Q 2 2tt 



(P | ^(r)7 + ^(0) 



(6.14) 



We shall take the mass of the state | P) to be M and the electron and positron mass to 
be m. In the weak coupling (non-relativistic) limit, the helicity dependence of the wave 
function factorizes away, so it is sufficient to consider one helicity sector. Here we shall 
take the two particle state with a\ and a 2 up. 

For a positronium state, | P) given by Eq. ( |6 . 1 1| ) we obtain 



F L(M X \iag 



+ {M 2 ) diag 

4 r 

= — I dxid 2 K,H(K,i-) 2 + m 2 )\ ip 2 I S(x — Xi) + 5(1 — x — x{) 
Q 2 J L 

+ -prz X! / dx\d 2 K^ I dx 2 d 2 n 2 \ ip3 | 2 ((^i - ) 2 + m 2 )5(x — Xi) + 



Q 2 



((k^) 2 + m 2 )5(x -x 2 ) 



The off-diagonal contributions to F[,£(x) comes from Mi alone. 

4 e 2 



where 



M{ 



ff-diag 



1 



Q 2 2(2ny 



dxxd'nf / dyd'n 1 - M" + M\ + M{ + Mf (6.16) 



E( Xl -y) 2 



2(Ki) 2 y 



5(x — xi) 



2(k ± ) 2 x 1 



(6.15) 



8{x-y) \W(xiA) I , (6-17) 



Ml 



E(x — y) 2 



- (ni) 2 6(x - xi) + (n^dix - y) 



(6.18) 



Ml 



I r2(4) 2 y r/i . 2(k ± ) 2 (1-x 1 ) r/ 1 

1 -<J(1 - 2; - zi) - — — — ^(ar-l + y) 



E'( 2 /-x 1 ) 2 L(l-x 1 ) 



|V>2 1<72 (zi,«r) 



(6.19) 



E'{y-x 1 ) 



{k ± ) 2 5{x + y-l)- (4) 2 ^(! - ar - zi) 

(vr iCT2 (^i ) «r)V'2 1<T2 (2/ ) « ± ) + ^-c.) 



(6.20) 
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In this calculation we have taken all operators to be normal ordered. Also, in Eqs. ( |6.17|) - 
( p.2C| ) we have neglected all mass terms in the vertex, since as described in Appendix 
C, these terms are suppressed in the weak coupling limit. The energy denominators are 
given by 

(/t- L ) 2 + m 2 (/t^) 2 + m 2 (ki - k l ) 2 



E = M 2 



y 1-xi xi-y 

(Ki) 2 + m 2 (K- L ) 2 + m 2 (k± - k M2 



E' = M 2 - v 1J ^ — '-— U J-. 6.21 

xi y yy-xx) 

We define the twist four longitudinal photon structure function as in the previous 
chapter, 

1 xP H 



?r=4( 



Q 2 2tt 



J dy e 



[P | (-)F +A (y-)F A -(0) + \g+-F^{y-)FUU) I F) 



2 



(P | F+ x (y-)F+(0) | P) + (y- - 0)1 . (6.22) 



(P+)2 

Now, from the definition of Fj7£(x) and F[^(x) it can be verified explicitly that 

F T lu){-x) = - F uj){ x )i (6-23) 
Fir 9 )(-x) = -F[^(x). (6.24) 

F£7gUx) has both diagonal and off-diagonal parts. As before, we take all operators to be 
normal ordered and we get 



— f dx x d 2 n{ f dyd 2 K L \ ip 3 | 2 — ^— r*(l -x x -x-y) 

Q J J (l — xi — y) 

dx\d 2 K^ I dyd 2 K~^~ip2( x ii 



X diag Q 2 J J (l — xi—y) 

4 4e 2 



Q 2 2(2ttY 



^{V^h ] —^S(xi-x-y). (6.25) 

(a?i - y) 1 



The second term in the right hand side is the contribution of the instantaneous interaction. 
The off-diagonal contribution is 

where 

4 6 f , _n I f , ,n I X 



G 1 = -—^Lj- f dx 1( ?Kt f dyM 



Q 2 2(2tt) 3 7 1 L J tf E( Xl -y) 2 
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Ht - h )- 2(4) 2 2(k- l 

y 



(xi - y) x\ 
2(4 -K^ 2 



+ 



xi-y) (1-2/)' 



(6.27) 



G 9 



4 e 2 



Q 2 2(2vr) 3 



dxid 2 Ki J dyd 2 ^ — 



x 



(y - xi)' 



4(4 - k^) 2 2(k ± ) 2 2(4) 



+ 



(y-xi) 0--y) (1-^0 

2(4 -k 1 ) 2 (k 1 ) 2 (4) 2 



1^(^,4) 



{y-xi) y 



(r/ ia2 (xi, KtWiv, ^) + h.c)]8{x 1 +x-y) (6.28) 



where E and E' are given by Eq. ( |6.21| ). 
From these expressions, we calculate 



where 



and 



L{gY 



-dx 



x 



, dxd K 1p 2 *^2 



+ 



+ Q2^2 J dxd 2 ^ J dyd 2 q ± ip 3 *ip 3 



x 

2 f„±\2 



+ 



1 — x 
-K 1 - q L ) 2 



+ 



x y (1-x-y) 

.j J dxd 2 K^ J dyd 2 q ± [A 1 + ^2 + 51 + B 2 ] 



A e 2 

Q 2 ^) 
4 4e 2 



Q 2 2(2vr) 3 

ip*(x, K ± )^ 2 (i/,g ± ) 



y dxd 2 ^ J dyd 2 q ± - 



yf 



A, 



1 1 

E(x-y) 



-4, 



^ i r 2 ia2 (x, 4) i + ^r iCT2 «- L )^ ioa (y, q a 



i i 

&{y-x) 



y / 2 ^r i<T2 (^,^)^2 1<72 (y ) ^) 

+ V' 1 |^ ffa (z,/c- L )f 



(6.29) 



(6.30) 



(6.31) 



Si 



1 1 

E(x-y) 



+ w j2 (x,4)^r iff2 (2/,g ± ) 

71 



(6.32) 



Bn 



E> (y - x) 



Here Vi, V2, V\ and V'i are given by 



Vx 



(x - yf 



+ K x T V) ( 
(t. — 1A v 



(a; — 2/) ^ x 2 



K 2 = V, 



(^) 2 (l-2x) + (q ± ) 2 (2y-l) 2(«- L - g^) 2 



x(l - x)(x - y) y(l-y)(z-y) (2 - y) 2 
2(k ± - ? ± ) 2 (2-x-y) f (k 1 ) 2 t (q^ 2 



(6.33) 

(6.34) 
(6.35) 
(6.36) 



(x-t/) 2 (l/-x) V(l_ x )2 (1-y) 2 - 

In these expressions we have kept only those terms in the vertex which survive in 
the non-relativistic limit. The helicities a\ and <Ji are both up and we have neglected all 
mass terms in the vertex since in this limit they are suppressed. Also in the weak coupling 
(non-relativistic) limit, we consider only photon exchange interactions and the off-diagonal 
terms proportional to | ^2 1 2 originating from self energy effects can be neglected. 

Now, from the expression of ^3 given in Appendix C, we get 



dxd 2 ^ I dyd 2 q ± \ ip 3 



dxd k 



2(2tt)3 J ■" J * ^E" 1 1 E '~'l (6 ' 37) 
where B\ and B2 are given earlier. Using this, one can write the second term in the right 
hand side of Eq. (|6.29|) as 

f r 1 1 

dxd 2 K L I dyd 2 q ± — Bi H B? 

J \-E E' 



g 2 2(2vr) ; 



+ m 2 (g ± ) 2 + m 2 



+ 



+ 



q 



x y (1-x-y) 

Considering the fact that the total energy is conserved, one can write this as 



Q 2 2(2nY 



dxd 2 ^ J dyd 2 q ± (B 1 + B 2 ). 



(6.38) 



(6.39) 



So we get 



iF[^(x) + F, 



T=4l 

L(g) [ 



X 



X 



dx = — [ dxd 2 K, ± ip2*' l p2\- 
Q 2 J 



K 



X 



1 — X 



Q 2 2(2tt) 3 

4 4e 2 



dxd 2 K^ J dyd 2 q ± [A 1 + A 2 \ 

[ dxd 2 ^ [ dyd 2 q ± - — 

J J yX y) 



Q 2 2(2vr) 3 
i/jz(x, fir 1 , 1 - x, -K^^iy, q X ) 



(6.40) 
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Also, if we denote by M 2 , then in the non-relativistic limit, it can be shown that 



M 2 — Mq ~ 0(e 4 ). 



(6.41) 



So we neglect this difference in the energy denominators and replace the bound state mass 



in E and E' by M 2 = 



(K x ) 2 +m 2 
x{\—x) 



. The energy denominators then become 



E 



x 



(q L ) 2 + m 2 



Q- 



x-y 



,m 



(x-y) x 



) 2 (x-y) 2 + (^-q^} 



(6.42) 



and 



E' 



to 



m 



1 — x 
1 



^ [( _^_ ) ^_ !/f+(K ,_ ?J 



(6.43) 



We get, in this limit, 

\x) 



I 



x 



-dx 



4 

Q~ 2 



dxd 2 K^ \ 1p2 



|2l« 



TO 



x(l — x) 



4 2e 2 



g 2 2(27r)= 

TON 2 



/ 7/i 

V T 



dyd g $J(a;,K ,l-x,-K )ifa(y,q ,l-y,~Q ) 
1 



X J ( K ±- q ±)2 + ^y (x _ y) 2 

/ TO \2 I 



(6.44) 



Here in the weak coupling limit we have omitted the spin indices. We can see that the 
interaction part of / L dx has exactly the same form as in non-relativistic bound 
state equation (see Appendix C) and hence in the leading order it can be related to 
the expectation value of the Coulomb interaction. Since k 1 - and q 1 - are small in the non- 
relativistic limit, all (ft -1 ) 2 and (g -1 ) 2 dependence in the numerator of the interaction terms 
are neglected compared to the to 2 dependent terms. However, the term proportional to 
^bz^yr cannot be neglected because both x and y are almost equal and this term cancels 
the contribution from the instantaneous interaction in the non-relativistic limit. Both 
of these terms originate from Fu g -\ and one can see that only the gauge bosonic part 
of the longitudinal structure function is important for the Coulomb interaction in the 
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weak coupling limit. The m 2 terms in the energy denominators combine with the other 
terms to give the non-vanishing contribution. The fermionic part of the longitudinal 
structure function gives contribution to the kinetic energy of the fermions. Reminding 
oneself that we are working in the light-front gauge and not in the Coulomb gauge, this is 
a manifestation of the gauge invariance of the separation of the Hamiltonian density into 
a fermionic and a gauge bosonic part. 

The Fermionic part of the Hamiltonian density is given by 

1 



ef 



f — m/>7 d ip 



2^ +t la ± . (id 1 - + gA 1 -) +*y m\-j-\ a ± .(id ± + gA 



+ 7°m 



V> + .(6.45) 



The gauge bosonic part of the Hamiltonian density is given by 



or = -f +x f; + \g + -{F Xa f = \(d+A-) 2 + If* 

= (d l A 3 ) 2 + 2ed l A l (J^j 2(^ + )V + 



The fermionic part of the longitudinal momentum density is given by 
The gauge bosonic part of the longitudinal momentum density 



(6.46) 



(6.47) 



0++ = —F +X F^. 



(6.48) 



For a positronium like bound state, we calculate the matrix element of and . 
The matrix elements have both diagonal and off diagonal contribution. The diagonal 
contribution to the matrix element from the fermionic and the gauge bosonic part is 
given by 

rpM2 



(p | e + - 1 p) 



(P+) 



2 / dxid 2 ^^^ 



;(P | e ++ | P) 



+ 



diag 



x-y (1 — ari) 



+ 



2 I dxid k 1 I dx2d k 2 ^3^3 
8e 2 



(4) 2 + (4) + (-4 - 4) 



Xi 



x 2 



2(2tt) 



(1-Xi- x 2 ) 

1 



,, , dx\d 2 n^ j dyd 2 K^i\)\{x\, K,f)^ 2 (y, 



[x-yf 



(6.49) 
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where 0+" = Of~ + 9+~ . 

The off-diagonal part can be written as 

,x ( PM2 



(P | 



(i 3 



+ ^2 



p I e ++ | p) 



= V! + v 2 



off-diag 



where, 



V 1 = ^fdx<P K Ljdy<Pq 



1 1 

E{x-y) 

+V 2 (r2 aia2 (x,K ± )^(y,q ± ) + h.c.) 



(6.50) 



(6.51) 



V2 = ^I^S^^J^x) I 2 ^" 

+n(^ 1<T2 (^, i - ^) + /i.e.) 



(6.52) 



Here the expressions for V 1: V 2 , V\, V' 2 , E and E' are given earlier. As before, we have 
taken the two particle state with both <j\ ,a 2 up. 

Considering only the photon exchange interactions and putting M 2 = Mq in the energy 
denominators as before, one obtains in the non-relativistic limit for a weak coupling theory 



(P | 9 + ~ | P) - ^<P | 6++ | P)]=2f dxd 2 ^\ ^ 2 



(P+) 



(/sr 1 ) 2 + m 2 
x(l — x) 



-2 



2e 2 



2(2tt) ; 



/W^,^i-x,-^(y,^i-y,-rt 



,m 



1 



J (^-^1)2 + (2)2^.^ 



+ (- 



m 



;) 2 



■(1-*)' (^-^) 2 + (( I ^y) 2 (x-^ 
Introducing the three vector p (see Appendix C), this can be written as 

(P^ 2 



(6.53) 



(P | + - | P) 



(p+y 



(P\9 ++ \P) = 2 d 3 p\ <i>(p) | 2 4[(p) 2 + m 2 



d 3 p / d 3 p'(p*(p)(p(p') — 



(6.54) 



2(2tt) 3 7 1 7 x , v '(^_ p ,y 
Multiplying the bound state equation (see Appendix C) by <p*(p) and integrating we get 

M 2 = J d?p\ </>($ \ 2 4[(p} 2 + m 2 ] - [ d*p [ rfV0* ^(pOt^W (6-55) 



2(2tt) 3 
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(p-piy 



Hence from Eqs. ( |6.44| ), ( |6.53 ) and ( 3.55 ) it can be seen that the sum rule is satisfied in 
the lowest order in weak coupling limit for a positronium target in light-front QED and 
it can be written as 



x ~ Q 2 



[P I # H 



P) 



(p+y 



[p l e A 



In the non-relativistic limit for a weak coupling theory, 



M 2 = 4m 2 + 4m.Be 



where B e is the binding energy of positronium. 
From Eq. ( |6.55| ) we obtain 

i) 2 2e 2 



i M 2 
P)\ = 4—. (6.56) 



(6.57) 



B P 



m 



— / d 3 p / d 3 p><P*(p)<P(p>) (6.58) 
tt) J J (r> — n i 



2(27r) 3 J 'J * ' - '{p-p'f 

The first term in the right hand side is the kinetic energy with y being the reduced mass 
of the two body system and the second term is the expectation value of the Coulomb 
interaction. So we see that in the weak coupling limit, the sum rule reduces to a relation 
connecting the kinetic and potential energies to the binding energy. 



6.3 F[ 4 for the Ground State of Positronium 

The bound state equation (see Appendix C) can be analytically solved for QED, 
which is the primary motivation for studying QED. The ground state wave function of 
positronium is given by 



<f>v, Sa , Se ,{p,s, s') = <f) u (p)6 Se>s 5 SehS 



(6.59) 



where s e and s e i label the spin quantum numbers of the electron and positron respectively 
and v denotes all the other quantum numbers, v = n,l,m correspond with the standard 
non-relativistic quantum numbers of hydrogen. The spin part factorizes out and the 
normalization condition is given in Appendix C. The wave function is given by 



4(e n ) 



where 



((en) 2 + W) 



ma 



2n 
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(6.60) 



(6.6i; 



and 



are Hyperspherical harmonics. Here <\m \< I <n — 1. 
For Is state of positronium, we have 

Y 1 

r 1,0,1 



(6.62) 



In terms of x and the Is state wave function can be written as 

4( ei )l 



[ei) 2 — m 2 + 



1 (K ± )' 2 +m 2 
4 s) 



(6.63) 



(6.64) 



which agrees with ||] for non-relativistic x ~ |. 

The leading order contributions to the structure functions F 2 (x) and F£ =4: (x) can be 
directly evaluated using this wave function. 

1 2 



d K \ 1p2(%7 K ) 



Ax 4 (l - x) 4 



(6.65) 



m 2 [(l - 2x) 2 + a 2 x(l - x)] + (^) 2 

where A = 4.78 x 10" 12 (Met;) 6 . The integral is convergent and can be evaluated analyti- 
cally introducing a cutoff A and taking the limit A — > oo in the end. We get 

x 4 (l — x) 4 



Fo(x) = 28.15 x 10 



-ii 



/ l-2x) 2 + a 2 x(l-x] 



3 ' 



(6.66) 



Near x = ~, F 2 (x) ~ 1861.22x(l — x). F 2 (x) is very sharply peaked at x 
The twist four longitudinal structure function is given by 



x 



A(k x ) 2 x 3 (1 - x) J 



i • 



(6.67) 



m 2 [(l - 2x) 2 + a 2 x(l - x)] + (k^ 2 

Here we have considered only the (ft -1 ) 2 dependent part, the integral of which is directly 
connected to the kinetic energy of the electron-positron pair. Evaluating the integral 
analytically, we get 



Fr(x) 



14.7 x 10 



-ii 



x 3 (l — x)'- 



X 



Q 2 



1 - 2x) 2 + ct 2 x(l - x) 



(6.68) 
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Q 



2 F L (x) 



is sharply peaked at x — ^. Near the maximum 



X 



Q 2 ~ 14.7 x KT 11 — 

era; 1 — x 



(6.69) 



6.4 Discussions 

To summarize, in this chapter, we have investigated the twist-four longitudinal struc- 
ture function for a meson in 1 + 1 dimensional QCD and also for a positronium like bound 
state in light-front QED in the weak coupling limit. 

In the previous chapter, we have derived a new sum rule for the twist four part of the 
longitudinal structure function. The validity of the sum rule has been explicitly checked 
in two dimensional QCD ('t Hooft model). To get a qualitative picture of the twist four 
structure function we have computed numerically both F 2 and F L structure functions 
in the 't Hooft model for the ground state wave function calculated using a variational 
ansatz. 

In the weak coupling limit in light-front QED, we get expressions that look similar to 
the familiar non-relativistic expressions, but the entire calculation is fully relativistic in 
the leading order in bound state perturbation theory. We have explicitly verified a sum 
rule for Fl =A that we previously proposed. We have shown that in the weak coupling limit, 
the sum rule reduces to a relation connecting the kinetic and the potential energies to the 
binding energy of positronium. We have also shown that, in this limit, the fermionic part 
of F[ =4 contributes only to the kinetic energy of the fermions and not to the interactions. 

The twist four part of the longitudinal structure function is important since it is the 
leading non-perturbative contribution to Fl. The leading twist contribution to Fl is 
perturbative, in contrast to the case of F 2 . This analysis for a bound state in weak- 
coupling light-front QED is quite interesting since it gives an idea of what goes in such 
a calculation in light-front QCD. The situation in light-front QCD is somewhat different 
because similarity renormalization group technique generates a confining interaction even 
in 0(g 2 ) in the effective Hamiltonian which makes it impossible to solve the bound state 
equation analytically. However, similarity renormalization upto 0(e 2 ) does not produce 
any additional interaction in the effective QED Hamiltonian and in the weak coupling 
limit, one can work with the canonical Hamiltonian. Similarity renormalization results 
agree with the results using particle number truncation in the weak coupling limit for light- 
front QED. Inspite of these differences, the overall computational framework in QCD is 
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the same and this analysis in light-front QED allows an analytic understanding of the 
problem. 
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CHAPTER 7 



Transverse Spin in QCD and 
Transverse Polarized DIS 

The spin structure of the proton is one of the most challenging problems in present 
day particle physics. A large number of theoretical and experimental investigations are 
going on in this subject. The famous EMC result gave rise to the so called 'spin crisis' by 
indicating that only a small fraction of the proton helicity is carried by the quarks and 
antiqurks, and so Ellis- Jaffe sum rule is violated. A much more debated and interesting 
issue is the transverse polarized structure function g?, which can be measured in transverse 
polarized DIS. The contribution to the cross section in this case is proportional to gx 
and this is suppressed by a factor i compared to the leading contribution of g\ in the 
longitudinally polarized case. So, the structure function g T is measured less accurately 
than g\. Also, gT cannot be expressed as an incoherent sum over on-mass- shell partons 
like g\. The partons must be interacting in order to contribute to transverse polarized 
scattering. Due to these reasons, g? is called a higher twist effect. 

The literature on transverse polarized structure function g T or g 2 (gx = 9\ + 92) is vast 
and even contradictory and confusing |], [J . Early works on g 2 include [|3| . The discussion 
of #2 at large Q 2 using OPE has been done in M . Quark mass plays an important role in 
transversely polarized scattering (see chapter 4 for references). The recent polarized DIS 
data has opened up new avenues to explore gT- 

We analyze the transverse polarized structure function g? in light-front Hamiltonian 
QCD. We show that gT is related to the interaction dependent transverse spin operator 
in light-front QCD which is the reason for its complexity. In this chapter, we present our 
analysis of the transverse spin operator and transverse polarized structure function. 

The plan of this chapter is as follows. In Sec. 7.1, first, we briefly review the complex- 
ities associated with the description of the spin of a composite system in a moving frame 



81 



in the conventional equal time quantization. Then in Sec. 7.2 we give the canonical struc- 
ture of light-front Poincare algebra. In Sec. 7.3 and 7.4 we introduce the light-front spin 
operators for massive and massless cases respectively of arbitrary transverse momentum. 
The explicit form of transverse rotation operators in light-front QCD is derived in Sec. 
7.5. The connection of the transverse spin operators with g? is given in Sec. 7.6. Sum- 
mary and conclusions are presented in Sec. 7.7. For the sake of completeness and clarity, 
the explicit form of the kinematical operators and the Hamiltonian in light-front QCD 
starting from the gauge invariant, symmetric, interaction dependent, energy momentum 
tensor is derived in Appendix D. A complete discussion of transverse spin operators in 
free fermion field theory and free massless, spin one boson field theory is presented in 
detail in Appendices E and F. 

7.1 The Problem of Spin in Relativistic Quantum Field Theory 

From the early days of quantum field theory, it has been recognized that the issues 
associated with the spin of a composite system in an arbitrary reference frame are highly 
complex and non-trivial ||. The familiar Pauli-Lubanski operators readily qualify for 
spin operators only in the rest frame of the particle. For a single particle in a moving 
frame it is known |J how to construct the appropriate spin operators starting from the 
Pauli-Lubanski operators. How to construct the spin operators for a composite system 
in an arbitrary reference frame is a nontrivial problem. In equal-time quantization, the 
complexities arise from the facts that for a moving composite object, Pauli-Lubanski 
operators are necessarily interaction dependent and, further, it is quite difficult [^] to 
separate the center of mass and internal variables which is mandatory in the calculation 
of spin. Due to these difficulties there has been rarely any attempt to study the spin of a 
moving composite system in the conventional equal time formulation of even simple field 
theoretic models, let alone Quantum Chromo Dynamics (QCD). 

In order to resolve the above mentioned problems and puzzles, we have undertaken an 
investigation of the spin of a composite system in an arbitrary reference frame in QCD. 
We have compared and contrasted both the instant form and front form formulations. 
We emphasize that the interaction dependence of the spin of a composite system in an 
arbitrary reference frame is not a peculiarity of light-front dynamics, it is a general feature 
in any formulation of quantum field theory. What is peculiar to light-front dynamics is 
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that one can at most go only to the transverse rest frame of the particle. No frame exists 
in which P + = and one is so to speak "always in a moving frame". As a consequence, 
spin measured in any direction other than that of P + cannot be separated into orbital 
and intrinsic parts. This is to be contrasted with the light-front helicity J7 3 which is 
independent of interactions and further can be separated in to orbital and intrinsic parts. 
The situation is quite analogous to that of a light-like particle. In this case it is well known 
that since there is no rest frame, one can uniquely identify the spin of the particle only 
along the direction of motion since only along this direction one can disentangle rotation 
from translation for a massless particle. Also, in any direction other than the direction of 
motion, one cannot separate the angular momentum into orbital and intrinsic parts. 

7.1.1 Intrinsic Spin in Non-Relativistic Quantum Mechanics 

In the non-relativistic case, the transformation from one inertial frame into another 
is obtained by a Galelian transformation, which is generated by ten generators, namely, 
the Hamiltonian H, which generates time translation, three momenta P which produce 
translation in space, three angular momenta J which produce rotation and the Galelian 
boost operators N. Among these, only the Hamiltonian is dynamical, all the remaining 
generators are kinematical. These generators obey the following commutation relations 
among themselves: 

[Jk, Ji] Jmi [Pki Jl] — itfdmPm, [N k , Jl] 

[P k , Pi] = 0, [N k , Ni\ = 0, [Pk, Ni\ = tikiM 
[P k , H] = 0, [J k , H] = 0, [JVfc, H] = -iP k , (7.1) 

where we have used summation convention over repeated indices. Here M is the mass 
operator. From the above commutation relations, it can be seen that the angular momenta 
are not translationally invariant. The angular momenta can be written as 

J = ^(NxP) + I (7.2) 

where w is the inverse of the mass operator and I is called the intrinsic spin ||. The 
intrinsic spin operators are translationally invariant and obey SU(2) commutation rela- 
tions among themselves. The non-relativistic spin operators are simple due to the fact 
that the Galelian boost generators are kinematical and commute within themselves. 
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7.1.2 Relativistic Spin Operators in Arbitrary Frame 

A relativistic dynamical system in equal-time formulation is transformed from one 
frame to another by Poincare transformation, the generators of which are, the relativistic 
Hamiltonian H, three momenta P, three angular momenta J and the boost operators K. 
The commutation relations obeyed by these generators are given below: 

\Jki J{\ ^klmJrm \Pki Jj\ ^^klmPmi \Kki Jl\ ^klm^m 

[P k , Pi] = 0, [K k , Ki] = -ie k i m J m , [P k , Ki] = i5 kl H 

[P k , H] = 0, [J k , H] = 0, [K k , H] = -iP k . (7.3) 

The crucial difference from the non-relativistic case is that the relativistic boost generators 
are dynamical and they do not commute with each other. The angular momentum J is 
not translationally invariant and therefore does not qualify as spin operators. The Casimir 
invariants of the Lie algebra are 

M 2 = P"P M , W 2 = W^Wp, (7.4) 

where 

W» = -)f vpa J vp P c (7.5) 

is the Pauli-Lubanski operator. It is translationally invariant and its components obey 
the commutation relations: 

[W», W v \ = ie^ vXp W x P p . (7.6) 

For a particle at rest, the 4-momentum is, P M = (M, 0,0,0). The commutation relation 
then becomes 

[W\ W j ] = ie ijk W k M (7.7) 

so that if we define three operators, S l = then S l obey SU(2) algebra. Thus, we can 
define spin in terms of Pauli-Lubanski operators. However, in an arbitrary frame, it is 
much more complex, because W p are interaction dependent. 

We consider two classes of representations which are of physical importance: 

• Positive time-like representations: M 2 > H > 

• Positive light-like representations: M 2 = H > 0. 

In either cases we do not demand that the representations be irreducible (this allows us 
to deal with elementary and composite systems simultaneously). 
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A. Positive Time-like Representations 

Starting from the Pauli-Lubanski operator, one can construct an operator S such that 
it is translationally invariant, transforms as a three vector under pure rotations and its 
components obey SU (2) commutation relations among themselves. 

[S\ P"} = 0, [J j , S k ] = ie ]kl S\ [S j , S k ] = ie ]kl S l . (7.8) 

A suitable solution to the above requirements is provided by 



s - JL 

M 



w 



M + H 



- J P ° K x P (J P) P (79) 

where W are the space components of the Pauli-Lubanski operator. H, P are equal time 
Hamiltonian and momentum operators respectively obtained by integrating the energy 
momentum tensor over a spacelike surface and J and K are the equal time rotation and 
boost generators respectively, which are obtained by integrating the angular momentum 
density over a spacelike surface. Since boost K is dynamical, all the three components of 
S are interaction dependent in the equal time quantization. Nevertheless, the component 
of S along P remains kinematical. We shall see later in this section that, this is to be 
compared with light-front quantization where the third component of the light-front spin 
operator J7" 3 is kinematical. This arises from the facts that boost operators are kinematical 
on the light front, the interaction dependence of light-front spin operators J % arises solely 
from the rotation operators, and the third component of the rotation operator J 3 is 
kinematical on the light- front. 

The operators S cease to be defined when M tends to zero. The commutation relations 
among P, S and M are given by 

[P j , S k ] = 0, [S j , S k ] = ie jkl S l , [S j , M] = 0. (7.10) 

Since P and M stand for the momentum and invariant mass of the system, the above 
relations make clear that S should represent 'intrinsic spin' of the system. 
The invariant W 2 can be completely expressed in terms of M and S as 

W 2 = -M 2 S 2 . (7.11) 
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B. Positive Light-like Representations 

In the massless case, the spin operators do not obey SU (2) algebra. Beginning from 
the basic generators P, J and K (here H = |P|) one has to construct operators S, T l 
and T 2 such that they commute with four momentum P M and amongst themselves satisfy 
E(2) commutation relations: 

[s, r 1 ] = tT 2 , [s, r 2 ] = -iT\ [T\ r 2 ] = o. (7.12) 

A suitable solution consistent with the above requirements is: 

W° 



s 



p 



(I p I +^ 3 ) 



r2 jw 3 + w°) 



r = w ~ p JypY^pry ( 7 - 13 ) 

Note that here S is the component of angular momentum in the direction of motion. 
Another point is, S is a scalar under pure spatial rotation, while shows complicated 
behavior under pure boosts. 

C. Comments 

The spin operators for a single particle in an arbitrary frame can be defined easily in 
the above way. The generators for a multi-particle relativistic system have been analyzed 
by several authors J7j in the equal-time formalism. The expressions obtained are too 
complicated to be used in any practical calculations. In order to describe the intrinsic 
spin of a composite system, one should be able to separate the center of mass motion from 
the internal motion. Even in free field theory, this turns out to be quite involved (See Ref. 
|7]] and references therein). The generators cannot be easily separated into the center of 
mass and internal variables. Moreover, the derivations have been done neglecting the field 
theoretical effects such as pair creation and crossing and so are expected to be valid in the 
relatively low energy region where an expansion in - is permissible. Interactions are to 
be incorporatated by introducing an effective potential which vanish sufficiently rapidly 
for large distance. On the other hand, in light-front theory, since transverse boosts are 
simply Galilean boosts, separation of center of mass motion and internal motion is as 
simple as in non-relativistic theory. (See Appendix G). 

Next, we discuss the light-front case. 
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7.2 Light-Front Poincare Generators 



The light-front Poincare generators have been introduced in chapter 3. In terms of 
the gauge invariant, symmetric energy momentum tensor 0^", the four- vector P^ and the 
tensor M^ u are given by 



P» = -J dx-d 2 x ± e + ", (7.14) 

(7.15) 



M» v = - I <l.r-<f 2 .r 



The boost operators are M + ~ = 2K 3 and M +l = E % . The rotation operators are M 12 = 
J 3 and M~ l = F\ The Hamiltonian P~ and the transverse rotation operators F l are 
dynamical (depend on the interaction) while other seven operators are kinematical. The 
rotation operators obey the F(2)-like algebra of two dimensional Euclidean group, namely, 

[F\ F 2 ] = 0, [J 3 , F { ] = ie ij F j (7.16) 

where e* J is the two-dimensional antisymmetric tensor. The remaining non-zero commu- 
tators are: 

[P+, K 3 ] = -iP + , [P + , F 1 } = 2iP\ [P\ E j ] = i5 ij P + 

[P\ J 3 ] = [p\ F j ] = i5 ij p-, [K 3 , E*] = iE i 

[K 3 , F*] = -iF\ [K 3 , P~\ = -iP~, [E\ J 3 ] = ie ij E j 

[E\ F j ] = 2tK 3 5 ij - 2iJ 3 e i \ [E\ P' J ] = -2iP\ (7.17) 

Another important point is that, the transverse boost generators E % commute with each 
other like the Galelian boosts. 

From the early days of light-front field theory, the complications associated with trans- 
verse rotation operators F % have been recognized. They are interaction dependent just 
like the Hamiltonian. Furthermore, together with the third component of the rotation 
operator J 3 , which is kinematical, F l do not obey the angular momentum algebra. In- 
stead they obey the algebra of two dimensional Euclidean group which is appropriate only 
for massless particles. For massive particles, one can define transverse spin operators H] 
which together with the third component (helicity) obey the angular momentum algebra. 
However, they cannot be separated into orbital and spin parts unlike the helicity operator 
| It]]. Most of the studies of the transverse spin operators in light-front field theory, so far, 
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are restricted to free field theory |TT | . Even in this case the operators have a complicated 



structure. However, one can write these operators as a sum of orbital and spin parts, 
which can be achieved via a unitary transformation, the famous Melosh transformation 
12]. In interacting theory, presumably this can be achieved order by order ]13| 



m a 



suitable expansion parameter which is justifiable only in a weakly coupled theory. 

Knowledge about transverse rotation operators and transverse spin operators is im- 
portant for addressing issues concerning Lorentz invariance in light-front theory. Unfor- 
tunately, very little is known |l4j regarding the field theoretic aspects of the interaction 
dependent spin operators, We emphasize that in a moving frame, the spin operators are 
interaction dependent irrespective of whether one considers equal-time field theory or light- 
front field theory. To the best of our knowledge, in gauge field theory, the canonical struc- 
ture of spin operators of a composite system in a moving frame has never been studied. 
In this chapter, we present a systematic investigation of the spin of a composite system 



in a moving frame in QCD. We refer to our papers and preprints, We show 

that, in spite of the complexities, light-front field theory offers a unique opportunity to 
address the issue of relativistic spin operators in an arbitrary reference frame since boost 
is kinematical in this formulation. 

7.3 Transverse Spin Operators 
7.3.1 Massive particle 

The Pauli-Lubanski spin operator 

W» = -~ e WM vp P a (7.18) 
2 

with e + ~ 12 = —2. For a massive particle, the transverse spin operators || J % in light-front 
theory are given in terms of Poincare generators by 

MJ 1 = W 1 - P l J z = -F 2 P + + K 3 P 2 ~~E 2 P- - P l J z , (7.19) 

2 2 

MJ 2 = W 2 - P 2 J Z = -\f 1 P + - K z P l + -E 1 P- - P 2 J 3 . (7.20) 



The first term in Eqs. ( |7.19 ) and ( 7. 20] ) contains both center of mass motion and internal 



motion and the next three terms in these equations serve to remove the center of mass 
motion. 



8s 



The helicity operator is given by 

W + 



J 3 + jh( Elp2 



E 2pU 



(7.21) 



Here, J 3 contain both center of mass motion and internal motion and the other two terms 
serve to remove the center of mass motion. The operators J % obey the angular momentum 
commutation relations 



(7.22) 



They commute with P M . 



7.3.2 Massless Case 



Again, we start from the Pauli-Lubanski spin operator 

1 



M up p a . 



(7.23) 



For the light-like vector usually the collinear choice is made [|17], [HJ, namely, p + 0, 
p 1 - = 0. Then we get, W~ = 0, W+ = J 3 p+, W 1 = \F 2 p + , W 2 = -\F x p+. 

In free field theory, we have explicitly constructed the Poincare generators for a mass- 
less spin one particle in A + = gauge in Appendix F. Consider the single particle state 
| pX) with p x = 0. From the explicit form of the operators, we find that 



J 3 | p\) 
F i | pX) 



X\pX), 

= 0, 1 = 



1,2 



(7.24) 



since p ± = 0. 

For calculations with composite states (dressed partons, for example) we have to con- 
sider light-like particles with arbitrary transverse momenta. Let us try a light like mo- 



mentum P M with P 1 - 7^ 0, but P 
massive particle, 

W+ = 

w 1 = 

w 2 = 
w~ = 



3p+ 



rp 



so that P 2 = 0. Then we get, as in the case of 



E lp2 _ E 2pl^ 
1 



-F 2 P + + K 3 P 2 - -E 2 P 
2 2 



— F 1 P + - K^P 1 -, - 

2 2 
p 2 p l - F X P 2 - J 3 P- 



E 1 P- 



(7.25) 



89 



Thus even though W 1 and W 2 do not annihilate the state, we do get W J, W lt (= \W + W~ + 
\W~W + - (W 1 ) 2 - (W 2 ) 2 ) I kX) = as it should be for a massless particle. 

Just as in the case of massive particle, we have the helicity operator for the massless 
particle, 



J 



^1 = J3 + -L(E 1 P 2 -E 2 P 1 ). 



(7.26) 



In analogy with the transverse spin for massive particles, we define the transverse spin 
operators for massless particles as 



x = W l - P l J 



i 



(7.27) 



They do satisfy 



J 1 | kX) = 0, 

J 3 | kX) = X | kX), 



(7.28) 



where k is an arbitrary momentum. The operators J % and J 3 obey the i?(2)-like algebra 



J 1 , J 2 



J 3 , J 1 =iJ\ J 6 , J 



r3 n-2 



(7.29) 



7.3.3 Comments 



In order to calculate the transverse spin operators, first we need to construct the 
Poincare generators P + , P l , P~ , K 3 , E l , J 3 and F l in light-front QCD. The explicit form 
of the operator J 3 is given Ref. ||10|| . The construction of F % which is algebraically quite 
involved is carried out in the next section. The construction of the rest of the kinematical 
operators is given in Appendix D. In this appendix we also present the Hamiltonian in a 
manifestly Hermitian form. 

In order to have a physical picture of the complicated situation at hand it is instructive 
to calculate the spin operator in free field theory. The case of free massive fermion is 
carried out in Appendix E. In free field theory one can indeed show that (see Appendix 
E) J 1 | A;A) = \ J2x> °a'a I kX'). The case of free massless spin one particle is carried out 
in Appendix F. 
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7.4 The Transverse Rotation Operator in QCD 



In this section we derive the expressions for interaction dependent transverse rotation 
operators in light-front QCD starting from the manifestly gauge invariant energy momen- 
tum tensor. It is extremely interesting to compare and contrast the situation in the equal 
time and light-front case. The angular momentum density 



The rotation operators are J 1 = e^ k M^ k . Thus in a non-gauge theory, all the three 
components of the rotation operators are manifestly interaction independent. However, 
as we have seen in the last section, the spin operators S l for a composite system in 
a moving frame involves, in addition to J\ the boost operators K % = M° l which are 
interaction dependent. Thus all the three components of S l become interaction dependent. 

A gauge invariant separation of the nucleon angular momentum is performed in Ref. 
||19|| . However, as far the spin operator in an arbitrary reference frame is concerned, the 
analysis of this reference is valid only in the rest frame where spin coincides with total 
angular momentum operator and in an arbitrary reference frame the need to project out 
the center of mass motion, which is quite complicated in equal time theory is not empha- 
sized. Moreover, the distinction between the longitudinal and transverse components of 
the spin is never made. It is crucial to make this distinction since physically the longitudi- 
nal and transverse components of the spin carry quite distinct information (as is clear, for 
example, from the spin of a massless particle). Moreover, even for the third component 
of the spin of a composite system in a moving frame, there is crucial difference between 
equal time and light front cases. J7 3 (helicity) is interaction independent whereas S 3 is 
interaction dependent in general except when measured along the direction of P. 

In light-front theory, generalized angular momentum 



(7.30) 



In equal time theory, generalized angular momentum 




(7.31) 




(7.32) 



J 3 which is related to the helicity is given by 




(7.33) 
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and is interaction independent. On the other hand, the transverse rotation operators 
which are related to the transverse spin are given by 



F l = M~* = X - J dx-d 2 x ± [x-Q +i - x i e + - 



They are interaction dependent even in a non-gauge theory since +_ is the Hamiltonian 
density. 

In light-front theory we set the gauge A + = and eliminate the dependent variables 
tp~ and A~ using the equations of constraint. In this paper we restrict to the topologically 
trivial sector of the theory and set the boundary condition A l (x~ ,x l ) — > as x~ ,z — > oo. 
This completely fixes the gauge and put all surface terms to zero. 

The transverse rotation operator 



dx d 2 x~ 



are 



+i 



The symmetric, gauge invariant energy momentum tensor 

1- 



e^ 



i) - F^ a F 
1 



va 
X 



(F Xaa y + ^tD x -m)^ 



where 



(7.34) 



(7.35) 



- iff* +gA», 



jpuXa _ Q^j\ Xa _ Q X j\^ a -\- g f abc A^ b A Xc 

= d v A\ - d\A va + gf abc A ub A c x . 



Fx 



First consider the fermionic part of Q^: 

The coefficient of g^ u vanishes because of the equation of motion. 

Explicitly, the contribution to F 2 from the fermionic part of is given by 



(7.36) 



(7.37) 



dx d x 



1— 1 — 

x~-ip(-f + iD 2 + rfiD + )i) - x 2 -^ + iD~ + -f~iD + )^ 

— — 



where 



F(I) 



dx d x x 



^ +t ^ id 2 ip + + if> +t gA 2 il> + + -ipY id + if) 



1- 

T 



(7.38) 



(7.39) 
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TP 2 



- / d.r (l 2 .v--.v 2 



(7.40) 



We have the equation of constraint 



id + ip~ = [a 1 ■ (id 1 + gA 1 ) + 7 m]f 



(7.41) 



and the equation of motion 



1 



id-ij) + = -gA~^ + + [a 1 ■ {id 1 - + gA 1 ) + l°m} — [a 1 - ■ (id 1 - + gA x ) + 7°m]^ + . (7.42) 



Using the Eqs. ( [7.41] ) and (|7.42[) we arrive at free (g independent) and interaction (g 
dependent) parts of Fp. The free part of Fp is given by 



TP 2 

P F(free) 



- I (/./• d 2 x~ L { .r" 



x 



-(d 1 ) 2 + m' 



id+ 



-{d 



±\2 + m 2 



id+ 



+ 



+ m 



a^d 1 + id 2 



d+ 



id^ 



d 



eve 



We have introduced the two-component field £, 







The interaction dependent part of F F( js is 

F 2 F[I)mt = g J dx-d 2 x x x-£}A 2 £ 



+ \gj dx-d 2 x L [£}^[{-i^A l + A 2 )£] 



d+ 



(7.43) 



(7.44) 



(7.45) 



The interaction dependent part of Fpr U s is 

FklDmt = \9l dx-d 2 X^ [£t J_ [H(7 3 A 1 + A 2 )e] + _L [e t (w 3 A l + 



— <? y da; gPx^x 2 



+ (f^)^ 1 " + ^(^ ■ d^i* 1 ■ A^ 
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d+ 



d+ 



dx d x x 



We have introduced cr 1 = a 2 and cr 2 = — a 1 . 

Next consider the gluonic part of the operator F 2 : 



— dx d x 



x~Q+ 2 -x 2 Q+- 



where 



Q+2 p+\a p2a 



1 



-F +Xa F- x a + -g + -(F Xaa ) 2 . 



Using the constraint equation 

-3+ A- a = 8 l A ia + gf abc -^-(A ib d + A ic ) + 2g-^-^T a £), 



we arrive at 



Z?2 _ p2 _|_ p2 

g g(free) ' g(int) 



where 



(7.46) 



(7.47) 



(7.48) 



(7.49) 



(7.50) 



TP 2 



- J dx~d 2 x ± lx' (A ja d + d j A 2a - A 2a d + d j A ja + A ja d + d 2 A ja ) 



-x 2 (A ka (d 1 ) 2 A ka )^ 



2 / dx~d 2 x ± A 2a d 1 A la . 



The interaction part 
F 2 - i 



dx~d 2 x ± x~ I gf abc d + A ta A 2b A tc 
+ g{f abc -^:(A ib d + A ic ) + 2±@T°tj)d+A 



(7.51) 



l - J dx-d 2 x ± x 2 l2gf abc d i A ja A ib A jc +!jf abc f ade A ib A jc A td A je 
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+ 2gd i A ia -^(f abc A jb d + A jc + 2£}T a £ 
+ g 2 (f abc ^(A lb d + A lc ) + 2±?T a £) 

(f ade ^(A jd d + A je ) + 2^^T a £ 

So the full transverse rotation operator in QCD can be written as 



(7.52) 



F 2 = Fi+Fh + F?n, 



(7.53) 



where 



F? 



Fh 



- J dx~d 2 x\x~Vl ~ ^C^o + V)] 



(7.54) 



- I dx d 2 x 
2 



a 6 d L + id 2 



H — I dx d 2 x ± m 
2 1 



a 



eve 



'in 



+ 2y dx~d 2 x ± g 
J dx~d 2 x ± 2(d 1 A 1 )A 2 
- l -Jdx~d 2 x^g^T^)A 



,(7.55) 



- - J dx-d 2 x ± gf abc -^{A lb d + A ic )A 2a 



(7.56) 



where Pq is the free momentum density, 7i is the free Hamiltonian density and V are 
the interaction terms in the Hamiltonian in manifestly Hermitian form. The operators 
Ffj and Fj n whose integrands do not explicitly depend upon coordinates arise from 
the fermionic and bosonic parts respectively of the gauge invariant, symmetric, energy 
momentum tensor in QCD. From Eq. ( 7.19| ) in Sec. 7.3 it follows that the transverse spin 
operators J\ {i = 1, 2) can also be written as the sum of three parts, J} whose integrand 
has explicit coordinate dependence, J % u which arises from the fermionic part, and J\ u 
which arises from the bosonic part of the energy momentum tensor. 

7.5 Connection with Transverse Polarized Scattering 

From the phenomenological point of view, the issue of transverse spin has become 
very important in high energy physics thanks to recent experimental advances . Since 
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transverse spin for a free massless gluon is identically zero, transverse spin measurements 
for gluonic observables directly probe the long distance, nonperturbative features of QCD. 
Analogous to longitudinally polarized scattering, where quark helicity carries roughly only 
25 % of the proton helicity, one may ask what is the situation in transversely polarized 
scattering. In particular can one relate the operators appearing in the transverse spin to 
the integrals of structure functions appearing in transverse polarized scattering? 

In the previous section, we have shown that, though the transverse spin operators 
cannot be separated into an orbital and a spin part, one can still define a decomposition 
of it into three different parts. In this section, we establish the physical relevance of 
this decomposition by exploring the connection between hadron expectation values of the 
transverse spin operators and the quark and gluon distribution functions that appear in 
transversely polarized deep inelastic scattering. 

It is known that the transverse polarized distribution function in deep inelastic scat- 
tering is given by (we have taken transverse polarization along the x-axis) 

9t(x) = ^jdr 1 e-^x(PS 1 \m( 1 l -^ 1 + ) 1 ,m + 

where P^ and are the four momentum and the polarization vector of the target. Using 
the constraint equation for ip~ , we arrive at 

/+oo r+oo 
dxg T (x) = / dx(g T{I) (x) + g T{II) (x)) (7.58) 
-oo J— oo 



+00 I 

dxg n i)(x) = 7^(PS 



PS 1 ). (7.59) 



+ ±.[$(io*A 1 + A 2 )]t] 



d+ 

Thus the integral of gT(i)(x) is directly proportional to the nucleon expectation value of 
Fjj. Both gr(i) and Ffj depend on the center of mass motion whereas both g T and J 1 
are independent of the center of mass motion. The removal of the center of mass motion 
from gr(i) is achieved by gT(ii)- We have 

r+oo 1 pi 

J dxg T{II) (x) = Mj^iPS 1 I I PS 1 ). (7.60) 
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The integral of gT(ii){x) is directly proportional to the nucleon expectation value of the 
quark intrinsic helicity operator 



J%) = \j dx-d 2 x^a% 



(7.61; 



Consider the polarized gluon distribution function that appears in transversely polar- 
ized scattering (see Ref. 



Gt(x) = ^-^I ter^SaPp J d^iPS 1 \ F + ;( v )F +a v (0) \ PS 1 ). (7.62) 



Here is the polarization vector, 



S" 4 = -uj^u 



(7.63) 



where u is the Dirac spinor. For longitudinal polarization, the normalized spinors are [£2] 



MP) 



V2P 1 



( P+ + M \ 
P l + iP 2 
P+-M 
V P 1 + iP 2 J 



ui(P) 



( -P l +iP 2 \ 
P+ + M 
P 1 - tP 2 
V -P+ + M J 



(7.64) 



These give S + = P + , S 1 - = 0. For transverse polarization (along x direction), the 
normalized spinors are 



MP) 



1 



V2P 1 



( M + P + -iP 2 \ 
-P 1 
P 1 

V —M + P + + iP 2 J 



, n(P) 



( —M + P + - 1P 2 \ 
-P 1 
P 1 

V M + P + + iP 2 J 



. (7.65) 



Using these, we get S + = 0,S 2 = 0, S 1 = M, S~ = 2^M. 

For a transversely polarized nucleon, F + _ = 0. Since for a, /j,,v — —, the contribution 
is automatically zero, (3 = — . Further, let us pick, without loss of generality, the transverse 
polarization along the x axis. Thus a is forced to be 1 or +. Then 



G T {x) = G T{I) (x) + G T (ny 



x 



(7.66) 



where 

G T (i)(x) 

and 



8nxMP^ 



G T (ii){x) 



J d^e-^iPS 1 I F + %{r])F +a + (0) - F +a + {r])F + %{0) \ PS 1 ), (7.67) 



mJp+Mi s ~ I dr]e ~ lvx{PSl 1 F+a ^ F+a ^ 

- F + %{r])F + i(0) I PS 1 ). (7.68) 
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We get 

/+oo I 
^ dxG T(I) (x) = —(PS 1 | A* a (0)-d+A- a (0) | PS 1 }. (7.69) 

From the constraint equation, we explicitly see that the operator structure of the integral 
of Gt(i) is similar to Fj n . 
We also have 

/■oo 1 pi 

^dxGr^Cx) = ^p+iPS 1 I (Ald+A* - A a 2 d + Ai) \ PS 1 ). (7.70) 

Thus the integral of G T (ii)(x) is proportional to the nucleon expectation value of the 
gluon intrinsic helicity operator 

(7.71) 

Thus, provided the interchange of the order of integrations is legal, we have shown that a 
direct relation exists between the coordinate independent part of J % which arises from the 
gauge invariant fermionic and gluonic parts of the symmetric energy momentum tensor 
and the integrals of the quark and gluon distribution functions Qt and Gt that appear in 
polarized scattering. 

7.6 Summary 

In this chapter, we have investigated the transverse spin operator in QCD and showed 
its connection with the transverse polarized structure function g T . In equal time quanti- 
zation, one encounters two major difficulties in the description of the spin of a composite 
system in an arbitrary reference frame. They are 1) the complicated interaction depen- 
dence arising from dynamical boost operators and 2) the difficulty in the separation of 
center of mass motion from the internal motion. Due to these severe difficulties, there 
have been hardly any attempt to study spin operators of a moving composite system in 
the conventional equal time formulation of quantum field theory. 

In light-front theory, on the other hand, the longitudinal spin operator (light-front 
helicity) is interaction independent and the interaction dependence of transverse spin 
operators arises solely from that of transverse rotation operators. Moreover, in this case 
the separation of center of mass motion from internal motion is trivial since light-front 
transverse boosts are simple Galilean boosts. 



T 3 

J 9{i) 



5 



dx d 2 x 



A\ 



A\cr 



A\ 
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We have investigated the case of transverse spin operators for both massive and mass- 
less particles. We have introduced the transverse spin operators for massless particles with 
arbitrary transverse momentum. This is done for the first time in light-front field the- 
ory. To provide physical intuition for transverse spin operators which have a complicated 
structure in interaction theory, we have provided the explicit form of these operators in 
Fock space basis for both free fermion field theory and free massless spin one field theory 
in the appendix. 

In QCD, our starting point is the formula for transverse rotation operators expressed 
as the integral of generalized angular momentum density given in terms of gauge invariant, 
symmetric, energy momentum tensor. We have emphasized the differences between spin 
operators in field theory in equal time and light-front quantization schemes. 

Appropriate to light-front quantization, we choose the light-front gauge. We use the 
constraint equations for ip~ and A~ to eliminate them in favor of dynamical degrees of 
freedom. In this initial study, we restrict to topologically trivial sector of QCD and set 
the requirement that the transverse gauge fields vanish as x~' % — > oo. This eliminates the 
surface terms and completely fixes the gauge. In the gauge fixed theory we found that the 
transverse rotation operators can be decomposed as the sum of three distinct terms: F\ 
which has explicit coordinate dependence in its integrand, and F % n and F\ u which have 
no explicit coordinate dependence in their integrand. Further, F l u and F\ u arise from 
the fermionic and bosonic parts of the energy momentum tensor. Since transverse spin is 
responsible for the helicity flip of the nucleon in light-front theory, we now have identified 
the complete set of operators responsible for the helicity flip of the nucleon. 

Our construction and decomposition of the transverse spin operators in QCD also have 
important phenomenological consequences. We have shown that nucleon expectation 
values of F % u and F\ n are directly related to the integrals of quark and gluon distribu- 
tion functions that appear in transversely polarized deep inelastic scattering. After the 
experimental discovery of the so-called spin crisis, the question of the sharing of nucleon 
helicity among its constituents has become an active research area. On the theoretical 
side, the first step involves the identification of the complete set of operators contributing 
to nucleon helicity. In this work, we have made this identification in the case of transverse 
spin. We have explicitly shown that the operators involved in the case of the helicity and 
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transverse spin are very different. Because of their interaction dependence, operators con- 
tributing to transverse spin are more interesting from the theoretical point of view since 
they provide valuable information on the non-perturbative structure of the hadron. 

It is extremely interesting to contrast the cases of longitudinal and transverse spin 
operators in light-front field theory. In the case of longitudinal spin operator (light-front 
helicity), in the gauge fixed theory, the operator is interaction independent and can be 
separated into orbital and spin parts for quarks and gluons. It is known for a long time 
that the transverse spin operators in light-front field theory cannot be separated into 
orbital and spin parts except in the trivial case of free field theory. We have shown that, 
in spite of the complexities, a physically interesting separation is indeed possible for the 
transverse spin operators which is quite different from the separation into orbital and spin 
parts in the rest frame familiar in the equal time picture. 

In short, in this chapter we have explored in detail the theoretical aspects of spin 
operators in quantum field theory in the context of QCD and their consequences. It is 
interesting to establish a transverse spin sum rule in analogy to the helicity sum rule and 
explore its phenomenological consequences. Since transverse rotational symmetry is not 
manifest in light-front theory a study of these operators is essential for questions regard- 



ing Lorentz invariance in the theory ]23|. An important issue in the case of transverse 
spin operators concerns renormalization. Since they are interaction dependent, they will 
acquire divergences in perturbation theory just like the Hamiltonian. It is of interest to 
find the physical meaning of these divergences and their renormalization. We address 
these issues in the next chapter by computing the expectation value of the transverse spin 
operators in a dressed quark state. 
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CHAPTER 8 



Transverse Spin in Light-Front QCD: 
Radiative Corrections 

In this chapter, we continue our investigation of the transverse spin operators in QCD. 
We have shown in chapter 7 that there exists a physically interesting separation of the 
transverse spin operators into a coordinate dependent part and two coordinate inde- 
pendent parts. We have also shown that the transverse spin operators are manifestly 
interaction dependent. They will acquire divergences in perturbation theory just like the 
Hamiltonian and it is important to find the physical meaning of the divergences and their 
r enormalizat ion . 

In light-front QCD Hamiltonian, quark mass appears as m 2 and m terms, m 2 in the 
free helicity non-flip part of the Hamiltonian and m in the interaction dependent helicity 
flip part of the Hamiltonian. It is known that m 2 and m renormalize differently, m 2 
and m also appear in J 1 . Do these terms require new counterterms in addition to those 
necessary to renormalize the Hamiltonian? 

Besides renormalization, there is another important issue regarding the transverse spin 
operators. Recently it was shown that starting from the manifestly gauge invariant 
symmetric energy momentum tensor, in light-front QCD (the gauge A + = and light- 
front variables), after the elimination of constrained variables, J7 3 becomes explicitly 
interaction independent and can be separated into quark and gluon orbital and spin 
operators. Thus one can write down a helicity sum rule which has a clear physical meaning. 
The orbital and intrinsic parts of the light-front helicity operator have also been analyzed 
recently in ||. Even though J % cannot be separated into orbital and spin parts and they 
are interaction dependent, one can still ask whether one can identify distinct operator 
structures in J % and whether one can propose a physically interesting sum rule. In chapter 
7, we have seen that, in the gauge A + = 0, by eliminating the constrained variables, one 
can decompose J into three parts, J % = J} + J l n + J\ n where only J\ has explicit 
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coordinate (x~,x l ) dependence in its integrand. The operators J\ l and J\ u arise from 
the fermionic and bosonic parts respectively of the gauge invariant energy momentum 
tensor. We have also shown its connection with transversely polarized scattering. One 
has to investigate whether this decomposition is protected by radiative corrections. 

In this chapter, we explore the theoretical consequences of the decomposition of J 1 . 
We compare and contrast the consequences of this decomposition and the corresponding 
decomposition of the helicity operator into orbital and spin parts. Next we address the 
issue of radiative corrections by carrying out the calculation of the transverse spin of 
a dressed quark in pQCD in the old-fashioned Hamiltonian formalism. To the best of 
our knowledge, this is for the first time that such a calculation has been performed in 
quantum field theory. This calculation is facilitated by the fact that boost is kinematical 
in the light-front formalism. Thus we are able to isolate the internal motion which is only 
physically relevant from the spurious center of mass motion. We carry out the calculations 
in a reference frame with arbitrary transverse momentum P 1 - and explicitly verify the 
frame independence of our results. We find that because of cancellation between various 
interaction independent and dependent operator matrix elements, only one counterterm 
is needed. We establish the fact the mass counterterm for the renormalization of J 1 is 
the same mass counterterm required for the linear mass term appearing in the interaction 
dependent helicity flip vertex in QCD. It is important to mention that the divergence 
structure and renormalization in light-front theory is entirely different from the usual 
equal-time theory. If one uses constituent momentum cutoff, one violates boost invariance 
and also encounters non-analytic behavior in the structure of counterterms ||. We have 
done one loop renormalization of the transverse spin operators by imposing cutoff on the 
relative transverse momenta and on the longitudinal momentum fraction. Upto one loop, 
we find that all infrared divergences (in the longitudinal momentum fraction) get canceled 
in the result. 

The plan of the chapter is as follows. In section 8.1, we discuss the physical relevance 
of the decomposition of the transverse spin operator and also compare and contrast it 
with the helicity operator. In section 8.2, we present the calculation of the transverse 
spin for a dressed quark state upto 0(a s ) in perturbation theory. Comparison of this 
calculation with the perturbative calculation of helicity is given in section 8.3. Discussion 
and summary are given in section 8.4. The calculation of transverse spin for a system 
of two free fermions is given in Appendix G. The full evaluation of the transverse spin 
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operator for a dressed quark in an arbitrary reference frame is given in appendix H. There 
we also show the manifest cancellation of all the center of mass momentum dependent 
terms. Some details of the calculation are provided in appendix I. We also refer to our 
original papers, |fj. ||]. 

8.1 Decomposition of Transverse Spin and Comparison with He- 
licity Sum Rule 

The helicity 3 3 and the transverse spin operators J % in light-front theory for a massive 
particle is given in terms of Poincare generators in chapter 7. In fl| it has been shown 
explicitly that the helicity operator 3 3 in the light-front gauge, in terms of the dynamical 
fields in the topologically trivial sector of QCD can be written as 

J 3 = Jfi + Jf + Jgi+ Jg (8.1) 

where 3fi is the fermion intrinsic part, Jj is the fermion orbital part, is the gluon 
intrinsic part and Jg Q is the gluon orbital part. The helicity sum rule is given by, for a 
longitudinally polarized fermion state 

^(PSW I J% + J% + 3% + J% | PSW) = ±\. (8.2) 

In the transverse rest frame (P 1 - = 0), J7 3 coincides with J 3 and the helicity sum rule 
takes the form 

| J% + J% + J% + 4 | P£«) = ±\. (8.3) 

For a boson state, right hand side of the above equation should be replaced with the 
corresponding helicity. Here, M is the normalization constant of the state. 

Unlike the helicity operator, which can be separated into orbital and spin parts, the 
transverse spin operators cannot be written as a sum of orbital and spin contributions. 
Only in the free theory, one can write them as a sum of orbital and spin parts by a 
unitary transformation called Melosh transformation. In chapter 7, we have shown that 
the transverse rotation operators can be separated into three distinct components 

F 2 = Fj + Ff, + Ff n , (8.4) 
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where 



F? 



Fh 



l - J dx-d 2 x ± [x-V$ - x 2 {H + V)] 



•5) 



- I dx d 2 x 
2 



a-'d' + id* 



+ 



dx d x m 



a 



id* 



id+ 



'in 



+ — J dx d 2 x ± g 
-J dx~d 2 x ± 2{d 1 A 1 )A 
Ijdx-d^g^T^A 



1 

|2 



^[(-^A 1 + A 2 )£] + ^(ia'A 1 + A 2 )]£ 



(8.6) 



- / dx-d 2 x ± gf abc ^-(A ib d + A ic )A 2a (8.7) 
2 J o + 



where Vq is the free momentum density, 7i is the free Hamiltonian density and V are 
the interaction terms in the Hamiltonian in full Hermitian form. It follows that the 
transverse spin operators JP, (i — 1, 2) can also be written as the sum of three parts, J\ 
whose integrand has explicit coordinate dependence, J l n which arises from the fermionic 
part, and J\ n which arises from the bosonic part of the energy momentum tensor. Thus 
although the transverse spin operators cannot be separated into an orbital and intrinsic 
parts, we can still talk of a decomposition into an orbital-like part and two intrinsic-like 
parts. 

At this point, we would also like to contrast our work with 0, where a gauge invariant 
decomposition of nucleon spin has been done. The analysis in has been performed in 
the rest frame of the hadron and no distinction is made between helicity and transverse 
spin, whereas, we have worked in the gauge fixed theory in an arbitrary reference frame. 

In analogy with the helicity sum rule, we propose a decomposition of the transverse 
spin, which can be written as 



Af 



{PS 1 | Jj + Jjj + J} H | PS 1 ) = ± 



1 



(8.8) 



for a fermion state polarized in the transverse direction. For a bosonic state, right hand 
side will be replaced with the corresponding transverse component of spin. 

The question that now naturally arises is, does the above decomposition have any 
physical relevance ? In the case of the helicity operator, we know that the hadron ex- 
pectation value of the fermion intrinsic part J% is related to the first moment of the 
quark helicity distribution function g\ measured in longitudinally polarized DIS. Also 
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is related to the first moment of the intrinsic gluon helicity distribution. In light-front 
QCD, helicity is kinematical and independent of interaction, as a result, the helicity sum 
rule has parton interpretation, that is, the different terms can be expressed as parton 
helicity distributions. But the transverse spin is dynamical and the corresponding sum 
rule has no parton interpretation. However, we have shown in chapter 7 that the above 
decomposition has relevance in the context of transversely polarized DIS. There exists a 
direct connection between the hadron expectation value of the fermionic intrinsic-like part 
of the transverse spin operator J l n and the integral of the quark distribution function 
that appear in transversely polarized deep inelastic scattering. Also we can identify 
the operators that are present in the hadron expectation value of J\ u with the operator 
structures that are present in the integral of the gluon distribution function that appear 
in transverse polarized hard scattering. The physical relevance of the decomposition is 
made clear from the identification. Our results show the intimate connection between 
transverse spin in light-front QCD and transverse polarized deep inelastic scattering. As 
far as we know, such connections are not established so far in instant form of field theory 
and this is the first time that the first moment of is related to a conserved quantity. 
Another important point is that in perturbation theory, the helicity flip interactions which 
are proportional to mass play a crucial role both in g? and in the transverse spin operator 
whereas they are not important in the case of the helicity operator. In fact, one cannot 
get a transversely polarized state in the absence of mass. 

Because the transverse spin operators are interaction dependent, they acquire diver- 
gences in perturbation theory. One has to regularize them by imposing momentum cutoffs 
and in the regularized theory the Poincare algebra as well as the commutation relation 
obeyed by the spin operators are violated ||. One has to introduce appropriate countert- 
erms to restore the algebra. In the next section, we perform the renormalization of the 
full transverse spin operator upto 0(a s ) in light-front Hamiltonian perturbation theory by 
evaluating the matrix element for a quark state dressed with one gluon. This calculation 
also verifies the relation ( |8.8| ) upto 0(a s ) in perturbation theory. 

8.2 Transverse Spin of a Dressed Quark in Perturbation Theory 

In this section, we evaluate the expectation value of the transverse spin operator in 
perturbative QCD for a dressed quark state. 
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The dressed quark state with fixed helicity a can be expanded in Fock space as 

\P,a) = <j>W(P,a)\0) 

+ i dkt<fkj r dhtgkj / 2{27v)3p+6 3 {p _ h _ h) 



<j>° aiM (P,\k 1: ;k 2 )tf(k 1: a 1 )a t (k 2 , A 2 ) | 0). (8.9) 

We are considering dressing with one gluon since we shall evaluate the expectation value 
upto 0(g 2 ). The normalization of the state is given by 

(k', X' \k,X) = 2{27r) 3 k + 5 xy 5{k + - k' + )5(k^ - fc ,J -). (8.10) 

The quark target transversely polarized in the x direction can be expressed in terms of 
helicity up and down states by 

I k + , fc-L, S 1 } = k + , ft" 1 , T)± | k + , fc-L, j)) (8.11) 

with s 1 = ±m^, where is the renormalized mass of the quark. 

We introduce the boost invariant amplitudes $± and <3>^ iA2 (x, g- 1 ) respectively by 
A (A;) = $^ and 0^ Aa (fc; k u fc 2 ) = ^r^ lA2 (x, g x ), where a: = |£ and g x = fcj 1 - xP x . 
From the light-front QCD Hamiltonian, to lowest order in perturbative QCD, we have 

x(l-x) 1 



{q ± ) 2 + m 2 (l-x) 2 ^T^ 



x -J=T*x* 2-*— + u —*-a^ - oHm— x^ei)*** (8.12) 



x 



Here m is the quark mass and x is the longitudinal momentum carried by the quark. 
Also, 5" 1 = a 2 and a 2 = —a 1 . It is to be noted that the m dependence in the above 
wave function arises from the helicity flip part of the light-front QCD Hamiltonian. This 
term plays a very important role in the case of transversely polarized target states. Also, 
we have seen in chapter 7 that a transversely polarized dressed quark state cannot be 
obtained when quark mass is zero. 

For simplicity, in this section, we calculate the matrix element of the transverse spin 
operator for a dressed quark state in a frame where the transverse momentum of the quark 
is zero. It can be seen from the expressions of J % in chapter 7 that the sole contribution in 
this case comes from the first term in the right hand side, namely the transverse rotation 
operator. A detailed calculation of the matrix elements of the transverse spin operator in 
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an arbitrary reference frame is given in appendix H where we have explicitly shown that 
all the terms depending on P 1 - get canceled. 

The matrix elements presented below have been evaluated between states of different 
helicities, namely a and a'. Since the transversely polarized state can be expressed in 
terms of the longitudinally polarized (helicity) states by Eq. flS.ll| ), the matrix elements 
of these operators between transversely polarized states can be easily obtained from these 
expressions. 

Here, we have used the manifest Hermitian form of all the operators. It is necessary to 
keep manifest Hermiticity at each intermediate step to cancel terms containing derivative 
of delta function. 

The operator |F 2 P + can be separated into three parts, 

l _ F 2 p+ = 1 F 2p+ + 1 ^ p+ + If? u P+ (8.13) 

where Ff, Fjj and Ff n have been defined earlier. The matrix elements of the different 
parts of these for a dressed quark state are given below. The evaluation of the matrix 
element of |F/P + is quite complicated since it involves derivatives of delta functions. A 
part of this calculation has been given in some detail in appendix I. We first consider the 
operator 

l -F}P+ = ~F/(1)P+ ~ \fK2)P+ - if?(3)P+. (8.14) 

The first term contains the momentum density, the second and the third terms contain 
the free and the interaction parts of the Hamiltonian density respectively. The matrix 
elements are given by 

(P, a | ^F 7 2 (1)P+ | P, a') = (P, a\\j dxd? ''q 1 x~ P 2 V | P, a') 



-- £ / dxdWK^x-^ + h.c (8.15) 

spin 



(P,a | l -Fj{2)P+ | P,a'> = (P, a | ~ J dxd 2 q^ x 2 P Q l -P+ \ P,a') 



spin " ^ 



I ^.1 ' /V " ' \ x l-x 



+ % - £ / W 2 ?V^ A1 ^ + h.c. (8.16) 



spin 



dq 2 
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In the above two equations, both the single particle and two particle diagonal matrix 
elements contribute. Here, h.c. is the Hermitian conjugate, J2 sp in is t ne summation over 
ax, cr'i, Ai, A'i. Pq is the free part of the Hamiltonian density. The interaction part: 

(P,a | \f?(3)P + I P,a') = (P,a \ \ J dxd 2 q L x 2 PrJ-P+ \ P,o') 

= -7=1= E / dxdY^—( - Urxll^ ■ e 1 ) 

+ i&± ^ kA + fr.c.). (8.17) 

Pj~ t is the interaction part of the light-front QCD Hamiltonian density. The manifestly 
Hermitian form of it is given in Appendix D. The interaction part of the Hamiltonian 
contains qqg, ggg, qqgg, qqqq and gggg terms. Only the qqg part of it contributes to the 
dressed quark matrix element. 

The operator ^Fj I P + which originates from the fermionic part of the energy momen- 
tum tensor, can be separated into three parts, 



-FhP- 



l.tt 



where ^F^ nII P + is the explicit mass dependent part of the operator, ^F 2 ±II P + is the part 
containing derivatives with respect to x x and ^F 2 jjP + is the interaction part. The matrix 
elements are given by 

(P, * I \PlnP+ I P, O = f*r< + y E / dxd 2 q ^Zxxlv\^\, (8.19) 

spin 



(P,a\-F 2 ±II P + \P,a') 



\ Y: f dxd^KlxxUVx^x'l, (8-20) 

spin 



(P, a | -F 2 P + | P, a') 



1 <7 



4 -/2(2tt) 3 



E / ^ 2 ^7r^=(^r 

spin J V 1 J- 



xl(-^ 3 4 + eDx^i - -xl(^ 3 4 + ) X& i 



CTl A 



f8.21^ 



In Eqs. (|8.19|) and ( |8.20|) , contributions come from only diagonal matrix elements whereas 



Eq. (|8.21| ) contain only off-diagonal matrix elements. The matrix element of \F 2 U P 
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which comes from the gluonic part, is given by 

{P,a\ l -F 2 n P+ | P,a>) = —fJ— E / dxd^^=L 
2 V2(27r) 3 t&J vl^ 

{*?Z<*W=A + h.c) - J dxd^-^ £ A$- A < A . (8.22) 



x 



A,<71,<7j 

The first term in the right hand side is the off-diagonal contribution which comes from the 
interaction dependent part of the operator. The second term is the diagonal contribution 
coming from the free part. 

The expectation value of the transverse spin operator between transversely polarized 
states is given by 

(P, S 1 | MJ l | P, S 1 ) = (P, S 1 | \f 2 P + + K 3 P 2 - ^E 2 P- - P X J Z | P, S 1 ). (8.23) 

Since we are in the reference frame with zero P x , only the first term in the right hand 
side, i.e. the |F 2 P + term will contribute, as mentioned earlier. We substitute for $^ lA 
using Eq. ( |S.12|) . The final forms of the matrix elements are given by 

ma s ^ , rl-e 
~~'J 2 



(P, S 1 | Mj}{\) | P, S 1 ) = -^C f ln% jf ~ e dx(l + x), (8.24) 
(P, S 1 | Mj}{2) | P, S 1 ) = ^C/ln^ J''" dx(l - 2x), (8.25) 

(P, S l | MJj l (3) | P, S 1 ) = -^C f \n% dx(l - x) (8.26) 

where MJ}{l),MJ}{2) and MJ}{$) are related respectively to F 7 2 (l), Ff(2) and F 2 {3) 
defined earlier, [i is the hadronic factorization scale for separating the 'hard' and 'soft' 
dynamics of QCD, i. e. we have set a hadronic scale such that | q L f >> jj 2 >> m 2 . e is 
a small cutoff on the longitudinal momentum fraction which can be safely taken to zero 
at the end of the calculation. 

So we obtain, from the above three expressions, using Eq. ( fc>.14| ), 

(P, S 1 | MJ} | P, S 1 ) = -^CAn%. (8.27) 

Air /r 

The contribution to the matrix element of MJ} T entirely comes from Fjj. The various 
parts of this matrix element are given by 

1 I |2 , , Q_ 2 f 1 '" j 1 

> 2 



(P, S 1 | MJl u | P, S 1 ) = ~m| | 2 + ^C/ln^ /' ' dar-J-, (8.28) 



in 



(P, S 1 I MJ\ U | P, S 1 ) = -^CAn% f 1 6 dx(l - x), (8.29) 

AlT \X Je 

(P, S 1 I MJg H | P, S 1 ) = ^ch^l (8.30) 

where MJ'^jj, MJ^ ±II and MJ g x u are related respectively to F^ n , F\ n and F 2 n . In 
Eq. (|8.28|) we have to use the normalization condition 



|2 1 , Q 2 r^-t \+x 2 



$? =1" — C/hi-V / • (8-31) 

111 2ir 1 /j, 2 Je l-x ( J 

Upto 0(a s ), the normalization condition will contribute only in the first term of Eq. 
( ggg ). We get, from Eq. ( ggg) , 

1 , ma SA1 U Q 2 f l ~ e ^( 2 1 + x r 
__ 

It is clear that the singularity at x = 1 is canceled due to the contribution from the 
normalization condition. The overall contribution coming from MJ} T is given by 



(P, S 1 I MJl n | P, S 1 ) = \m + ^C f \n% f 1 " dx f-A_ _ I + f!} . ( 8 .32) 

2 47T /i z Je V 1 — X 1 — X I 



(P, | MJh | P, S 1 ) = y (l + ) , (8-33) 



(P, S 1 | Mj/ Zi | P, S 1 ) = ^^C/ln^ jT 1 '(I - x)dx. (8.34) 



which does not involve any x divergence. The matrix element of MJ~j n is given by 

2ma s ^ , Q_ 2 /•!-«, 

It is to be noted that all the contributing matrix elements are proportional to the quark 
mass. Among the different parts of the operator, only J^n an d a P ar t of the interaction 
terms in J\ are proportional to the quark mass m. These mass dependent terms flip the 
quark helicity. It is also to be noted that the terms proportional to m 2 do not flip the 
helicity. In all the other terms, though the operators do not depend on m explicitly, the 
contributions to the matrix elements arise from the interference of the m terms in the 
wave function of Eq. ( |8.12| ), with the non-m dependent terms through the different parts 
of the transverse spin operator. Since in light-front formulation, helicity and chirality are 
the same, these linear in m terms are explicit chiral symmetry breaking terms. From Eq. 
( F27D and Eq. ( ggg ) we find that 

ma s ^ 1 Q 2 
"> 2 



(P, S 1 I Mj} + Mj} n I P, S 1 ) = —C f \n% f 1 \l - 2x)dx = (8.35) 

47T U 2 Je 



112 



which means that the entire contribution to the matrix element of the transverse spin 
operator is given by 

<P, | Mji | P, S 1 ) = y (l + 3 ~^Cf^) ■ (8-36) 

This contribution entirely comes from MJ}^ Contribution from the orbital-like part 
{MJ}) exactly cancels the contribution from the gluon intrinsic-like part (MJ} n ). 

The renormalized mass tur of the quark is given in terms of the bare mass upto order 
a s in light-front Hamiltonian perturbation theory by ||, 

m R = m (l + ^C>ln^ . (8.37) 

In the light-front formulation of QCD, there are two mass terms in the Hamiltonian, one 
is quadratic in m which is present in the free part and does not break chiral symmetry, 
the other is linear in m which we discuss here and which explicitly cause chiral symmetry 
breaking. An important feature of light-front QCD is that, these two mass scales are 
renormalized differently even in the perturbative region. The renormalization of m 2 is 
different from the result stated above. 

Adding all the parts, for a dressed quark in perturbation theory upto 0(g 2 ), the 
expectation value of the transverse spin operator is given by 

(P, S 1 | MJ 1 | P, S 1 ) = (P, S 1 | Mj} + MJh + Mj} n | P, S 1 ) 

= ^- (8.38) 

It is important to mention that here we are calculating the expectation value of the 
operator MJ % . In order to extract the eigenvalue of J % one has to know the eigenvalue 
of M. Both MJ 1 and M are dynamical operators. However, in this case, the mass M in 
the left hand side in the renormalized theory is nothing but the renormalized mass of the 
quark, which therefore gets canceled from the above equation, and we get 

(P, S 1 \J 1 \ P, S 1 ) = (P, S 1 I J} + J}j + Jjjj I P, S 1 ) 

= \ (8-39) 

The identification of J with spin, therefore, requires knowledge of the mass eigenvalue, 
independently of the boost invariance properties of the light-front dynamics. 

We can explicitly verify the relation between the integral of gx and the expectation 
value of the fermion intrinsic-like part of the transverse spin operator to order a s in 
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perturbative QCD. The transverse polarized structure function for a dressed quark is 
given ]TU| by, 

e 2 m ( a s Q 2 vl + 2x-x 2 f 1 , ,1 + x' 2 



+ -*)]}, (8.40) 



x 



so we get 



j\ T {x)dx = ^(P,S l | MJh | P.S 1 ) (8.41) 



which explicitly shows the connection between the integral of the transverse polarized 
structure function and the matrix element of the fermion intrinsic-like part of the trans- 
verse spin operator. It can be seen that the operator structure of MJ^jj as well as the 



corresponding part of is related to the quark transversity distribution h\ [11 



8.3 Comparison with the Perturbative Calculation of Helicity 

J 3 

It is quite instructive to compare our calculation of the transverse spin of the dressed 
quark with the helicity of the dressed quark in perturbative QCD. All the operators 
contributing to helicity are kinematical (interaction independent) and hence all of them 
give rise to only diagonal contributions. Further, in this calculation mass of the quark can 
be completely ignored since they give rise to only power-suppressed contribution. In the 
massless limit, helicity is conserved at the quark gluon vertex. This means that the quark 
in the quark-gluon state has the same helicity as the parent quark. Since the transverse 
gluon carry helicity ±1, we get a non- vanishing contribution from the gluon intrinsic 
helicity operator. However, both the quark and the gluon in the quark-gluon state have 
non- vanishing orbital angular momentum due to transverse motion. Conservation of total 
helicity implies that orbital contribution has to cancel gluon intrinsic helicity contribution. 
This is precisely what happens [Q and we find that the total quark plus gluon orbital 
part exactly canceled the intrinsic gluon contribution and the overall contribution to 
the helicity is ±|, which entirely comes from the intrinsic part of the fermionic helicity 
operator. 

In contrast, in the case of transverse spin operator, it has both interaction independent 
and interaction dependent parts. The latter gives rise to off-diagonal matrix elements and 
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they play a very important role. Of special interest is the gluon intrinsic-like transverse 
spin operator. This operator gives vanishing matrix elements for a free gluon. However, 
since gluon in the quark-gluon state has intrinsic transverse momentum, both diagonal 
and off-diagonal terms give rise to non-vanishing contributions and we get a net non- 
vanishing matrix element for the gluon intrinsic-like transverse spin operator. However, 
we find that contribution from this matrix element is completely canceled by that from 
the matrix elements of orbital-like transverse spin operators. This is analogous to what 
happens in the helicity case. 

In this section, the calculation of the matrix elements has been done in the frame with 
P 1 - = 0. The complete calculation of the matrix element of the transverse spin operator in 
an arbitrary reference frame is given in appendix H. It is clear from the expressions there 
that all the terms explicitly dependent on P 1 - get canceled in the expectation value of 
MJ 1 . The parts that remain after the cancellation of the P x dependent terms are those 
given above. In the above expressions, we have used the manifest Hermitian form of the 
operators. We again stress the fact that this manifest cancellation of contributions from 
center of mass motion is typical in light-front field theory because the transverse boost 
operators are kinematical. The situation in the equal time relativistic case is completely 
different and there one cannot separate out the center of mass motion from the internal 



motion in a straightforward way even in the free theory case |12| because of the compli- 
cated boost generators. Due to the manifest cancellation of the center of mass momenta, 
J % can truly be identified as the transverse spin operator. 

8.4 Summary and Discussions 

In chapter 7, we initiated the investigation of the transverse spin operators in QCD. 
In this chapter, we have continued the analysis and discussed two very important aspects 
of it. In analogy with the helicity sum rule, we have proposed a decomposition for the 
transverse spin. Earlier we have shown the relationship between nucleon matrix elements 
of J l n and J\ n and the first moments of quark and gluon structure functions respectively, 
appearing in transverse polarized hard scattering. This is the first time that the integral of 
gr is related to a conserved quantity, namely the transverse spin operator. It is important 
to mention here that the proposed decomposition of the transverse spin operator will 
not be affected if one adds a total derivative term to the angular momentum density. 
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Such a term can at most produce a surface term which we are neglecting since we have 
restricted ourselves to the topologically trivial sector of the theory. We have started 
with the angular momentum density defined in terms of the symmetric gauge invariant 
stress-energy tensor, which is obtained from the Noether's stress-energy tensor by adding 
a total derivative term. Even though the angular momentum density differs from the 
Noether angular momentum density by a total derivative term, both give rise to the 
same generators. Another point worth mentioning is that we have worked in the gauge 
fixed theory. In the light-front gauge, A + = 0, the transverse spin operator can be 
separated into three parts, and J\j is related to the first moment of Qt measured in 
transverse polarized scattering, which is a gauge invariant object. This is similar to the 
helicity case, where only in the light-front gauge and using light-front quantization, the 
intrinsic fermionic helicity is related to the gauge invariant first moment of g\ measured in 
longitudinally polarized scattering. The corresponding gluon intrinsic helicity cannot be 
measured directly in polarized deep inelastic lepton-nucleon scattering but in some other 
process like polarized hadron-hadron scattering. A similar situation holds in the case of 
transverse spin. 

We have also discussed another important issue, namely, the renormalization of the 
transverse spin operators. In light-front theory, in addition to the Hamiltonian, transverse 
spin operators also contain interactions and have a complicated structure. They are 
dynamical operators. Since transverse rotational symmetry is not manifest in light-front 
theory, a study of these operators is essential for questions regarding Lorentz invariance 
in the theory ||. Since they are interaction dependent, they will acquire divergences in 
perturbation theory just like the Hamiltonian. It is of interest to find the physical meaning 
of these divergences and their renormalization. The renormalization of only the intrinsic- 
like fermion part of the transverse spin operator has been discussed in the literature so far. 
In this chapter, we have carried out the renormalization of the full transverse spin operator 
for the first time upto 0(a s ) in light-front Hamiltonian perturbation theory by evaluating 
the matrix elements for a dressed quark target. We have shown that the entire contribution 
to the matrix element comes from the fermion intrinsic-like part of the transverse spin 
operator and is equal to |. The contributions from J\ and J\ n exactly get canceled. 
Also, the mass of the quark is very crucial in this case, since the helicity flip interactions 
which are proportional to the quark mass play a very important role. However, the terms 
proportional to m? do not flip the helicity and do not contribute. Since helicity flip is 
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involved, we do not encounter any quadratic divergence unlike the case of renormalization 
of the light-front Hamiltonian. Further, we have compared and contrasted the calculations 
of transverse spin and helicity of a dressed quark in perturbation theory. 

In this chapter, we have also verified the frame independence of our results. We have 
explicitly shown that, in an arbitrary reference frame, all the terms depending on the 
center of mass momenta manifestly get canceled in the matrix element. The cancellation 
is as simple as in non-relativistic theory since boost is kinematical on the light-front. For 
future studies, it is an interesting problem to evaluate non-pertubatively |L3| the matrix 
element of the transverse spin operator in light-front QCD. Also, in this work, we have 
used cutoff on the relative transverse momenta and the small x divergence gets canceled 
in the one loop result. It is interesting to study the renormalization of the transverse spin 
operators using similarity renormalization technique ||. 
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CHAPTER 9 

Summary and Conclusions 



In this thesis, we have investigated the higher twist structure functions, namely the 
twist four longitudinal structure function Fl and the transverse polarized structure func- 
tion g T in the recently developed method based on light-front Hamiltonian QCD. Deep 
inelastic scattering is a light cone dominated process and the general features of DIS are 
most suitably realized in terms of parton distributions as provided by Feynman's parton 
model. The field theoretic realization of the parton model is possible in light-front Hamil- 
tonian QCD approach. It is known that the constituent quark model (CQM) where the 
hadrons are thought of as composites of only valence quarks and antiquarks, is very suc- 
cessful in describing the non-relativistic properties of hadrons. In QCD, this description 
is not possible because of the complicated structure of the vacuum. The simplicity of 
light-front vacuum gives us hope of reconciling CQM with QCD. Our approach is more 
intuitive than the conventional approach based on OPE. Furthermore, here we deal with 
the structure functions themselves, instead of the moments, which naturally arise in the 
OPE method. It is to be remembered that it is the structure functions which are actually 
measured in DIS experiments. 

Because our formulation is entirely different, we have obtained various new and inter- 
esting results. Also, it is possible to address the issues which are not easy to address in the 
conventional approach. Both the structure functions twist four Fl and gx contain non- 
trivial interaction dependence in the operator structure and therefore involve quark-gluon 
dynamics. In chapter 5, we have derived a sum rule which connects the twist four part of 
Fl to the mass square of the target hadron. This is a new sum rule and it has important 
consequences regarding the understanding of how the mass of the hadron is distributed 
among its quark and gluon constituents. A very important aspect of our approach is 
that we can analyze both the perturbative and non-perturbative contents of the struc- 
ture functions. The sum rule has both perturbative and non-perturbative consequences. 
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The Fock space expansion of the target state and particle number truncation allow us to 
express the structure functions in terms of light-front multiparton wave functions. The 
structure functions can be calculated once these wave functions are known. Major inves- 
tigations are under way to calculate these wave functions by solving the light-front bound 
state equation. In chapter 6, we have investigated the twist four longitudinal structure 
function for bound states like a meson in 1 + 1 dimensional QCD and also for positron- 
ium in 3 + 1 dimensional QED in weak coupling limit. These calculations give analytic 
understanding of the underlying non-perturbative dynamics. We have also verified the 
sum rule non-perturbatively in both these cases. 

Another most interesting issue that we have addressed is the transverse polarized 
stucture function g T . The operator here also involves non-trivial quark-gluon dynamics. 
In chapter 7 we have shown that qt is related to a part of the light-front transverse spin 
operator. The problem of relativistic spin operators for a composite system in an arbitrary 
reference frame is non-trivial in equal time theory because of the interaction dependence 
of the Pauli-Lubansky operators. One cannot separate the center of mass motion easily 
because equal time boost operators are dynamical. We have shown in chapters 7 and 8 
that the light-front formulation is a much more straightforward way to address the issue 
of spin since light-front boost operators are kinematical and one can separate the center of 
mass motion from the relative motion. In chapters 7 and 8 we have shown that though the 
light-front transverse spin operators cannot be separated into orbital and spin parts, there 
exists a physically interesting decomposition of them into three parts, one orbital-like and 
two intrinsic-like. We have worked in the light-front gauge. One of the intrinsic-like part 
is related to the fermionic part of the gauge invariant symmetric QCD energy momentum 
tensor, and this in turn is connected with gx- We have compared and contrasted this 
decomposition of the transverse spin with the corresponding decomposition of helicity for 
the nucleon. 

In perturbative calculations using field theory one encounters divergences. The actual 
results are, of course, observables and are finite. Therefore one has to renormalize the 
divergent results. Renormalization in light-front Hamiltonian theory is entirely different 
from the usual equal time case because it is manifestly non-covariant and also a gauge 
fixed theory. In chapters 5 and 8, we have renormalized the twist four longitudinal struc- 
ture function and the full transverse spin operator respectively in light-front Hamiltonian 
perturbation theory. Replacing the target by a dressed quark, we have shown that the 
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counterterm required to renormalize the twist four part of Fl has the same structure as 
the mass counterterm and that needed for the renormalization of transverse spin operator 
is related to the linear mass counterterm in light-front Hamiltonian perturbation theory. 
It is to be remembered that m 2 and m are renormalized differently in light-front theory. 
Thus, our results show that the regularization and renormalization of the higher twist 
structure functions Fl and are dictated by the renormalization of the light-front QCD 
Hamiltonian. 

Our alternative approach gives a new insight to DIS. We have shown that the twist 
four longitudinal structure function is related to the light-front QCD Hamiltonian density 
and the transverse polarized structure function is related to the transverse spin operator. 
Earlier, it has been shown that F 2 is related to the longitudinal momentum density and 
gi is related to the helicity. Thus it gives a unified picture of DIS and connects the 
electroproduction structure functions to the Poincare generators in light-front theory. It 
is important to mention that, these relations are non-perturbative in general. 

There are various new and interesting problems along this line that needs to be in- 
vestigated in future. One of the most challenging problems is to calculate the structure 
functions for a bound state like a nucleon or a meson in 3 + 1 dimensional QCD using the 
recently developed similarity renormalization technique [TJ. In our work, we have used 
large ultraviolet transverse momentum cutoff and a small infrared longitudinal momen- 
tum cutoff to regularize the theory. It is interesting to address the issue of renormalization 
of the transverse spin operator using similarity renormalization. Another possible way of 
a non-perturbative calculation of the hadronic matrix elements is the transverse lattice 
formalism [Q. 
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APPENDIX A 

Notations and Conventions 



Light front variables are defined as, 



x+ = x° + x 3 , 



x = x° - x 3 . 



(A.l) 



We denote the four vector x^ by, 



(A.2) 



Scalar product of two four-vectors, 

1 



x-y = -x^y + -x y^ - x ■ y 



The metric tensor is, 



9/j.f = 



( \ \ 
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0-10 

V o o o -i J 
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-1 





V 











-1 / 



Partial derivatives are given by, 



d + = 2<9_ = 2 



d 
dx~ 



(A.3) 



(A.4) 



(A.5) 



d- = 2d + = 2^. (A.6) 

Four dimensional volume element: d A x = ^dx + dx~d 2 x ± . Lorentz invariant volume ele- 
ment in momentum space 

We define the integral operator, 

^p/OO = \J dy~e{x- - y-)f(y-). (A.8) 
123 



jU/OO = 1 1 d y~ I - y~ I /GO- ( A -9) 



Unless otherwise specified, we choose the following convention for the gamma matrices: 



7 ± = 7 °±7 3 . (A. 10) 



where a are the Pauli spin matrices. 

75 = «7o7i7273, a = 7°7- (A. 13) 

Projection operators 

A ± = ^yV- (A.14) 
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APPENDIX B 



Fj^i Comparison with Other 
Conventions 



The hadron tensor relevant to unpolarized electron-hadron deep inelastic scattering is 
given by 

W>" = (-g^ + ^)w l (x,Q 2 ) + (P^-^)(P v -^)w 2 (x,Q 2 ).(B.l) 

This is valid for any Q 2 and here we have not made any approximations. Here W\ is 
dimensionless and W 2 has dimension The dimensionless functions 

F L (x,Q 2 ) =2[-H/ 1 + [M 2 -^|^]H/ 2 l, (B.2) 

and 

F 2 (x,Q 2 ) =uW 2 (x,Q 2 ) (B.3) 

are the unpolarized structure functions. 

Our notations are different from the definition of Close (l| where both W\ and W 2 are 
dimensionless. The longitudinal photo absorption cross section is given by, 

47r 2 a 

a L = —W L (B.4) 
K 

where K is the incident flux of photons and Wl = ^f = — W\ — ^W 2 + M 2 W 2 . 

A q 

2 

In the limit M = 0, we get Fl = (—Wi — ^W 2 ). Our definition of Fl is also different 



from that of Reya 



P'rwr = %Y x F ^ va - (R5) 



In terms of W\ and W 2 , 



v 2 

pReya = 2x \ Wl _ W 2 + M 2 W 2 ] - M 2 \W X ~ M* W 2 ] . (B.6) 

q 
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In the limit M = 0, we get, F L eya = xF^. Thus, P^P V projects out Fl from only in 
the massless limit. 

Next, we compare our conventions with those of Ellis, Furmanski and Petronzio ||. 
We write in terms of Fl and Ft where Ft = —■ 

1 a^a v 

2 q z 

P^P V 1 1 

+ \x + PV + P v <f) - -gnF T 

1 P.q 2P.q K H H ' 2 y 1 T 

+ M 2 W 2 {-g^ + ^). (B.7) 
q z 

EFP has used the same expression but M = 0. Towards the end, they have considered 
the mass effects separately. In our work, we have taken both the quark and the target to 
be massive from the beginning. 
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APPENDIX C 



Bound State Equation for 
Positronium in the Weak Coupling 

Limit 



The Fock space expansion of the positronium state is given in chapter 6. The rela- 
tivists LF version of the Schroedinger equation is, 

M 2 + (P 1 ) 2 

p- 1 p) = y j | P) (C.l) 

where P~ is the LFQED Hamiltonian. Substituting the Fock space expansion of the 
state in this expression and taking projection with the three particle state, 
6 t (A;i,o-i)rf t (A; 2 ,o- 2 )a t (A;3, A 3 ) | 0) we get 



ilj% ia2X3 (x,K 1 ;x 2 ,K 2 ;l-x-x 2 ,K 3 ) = Mi + M 2 , (C.2) 
where the amplitudes are given by 



M 1 = -(-) r—— ~ h -x 2 ,-4;x 2 ,4) (C.3) 
E ^2(2^)3 Vl - x - x 2 

and 

m 2 = I— ^=-=^_ w 2 ^2 l4 (^4;i-^-4) (c.4) 

with the energy denominator 

E=[M 2_ rn 2 + (nj) 2 _ m 2 + (4) 2 _ (4) 2 , (C 5) 
x x 2 1 — rr — x 2 
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The vertices are, 



W 



[p .Ki — vm) 



22 X 



(T + (J 



(a .«2 — im) 

1-22 



1X^.(6^)* (C.6) 



and 



w 2 = xU E [t- 

Here, 



<7~ 



[a .k 2 — im) [o .Ky—im) 



2 — 22 



•?'2 



1 — 2 



^X-aUeir- (C7) 



Kg K^ K 2 ■ 



(C.{ 



Also, taking projection with a two particle state b\k\, a\)d\k 2 , cr 2 ) | 0) we get 



M 2 (4) 2 + ^» 2 (4) 2 + ™ 2 - 



2i 



1 — 2i 



V>2 lCT2 (^i ; 4, 1 - zi, -4) =Mi + M 2 + I. (C.9) 



Here 



^= v r e dyfd 2 K x - 1 

[/ 1 ^3 1<J2A (?/, K" 1 , 1 - 2i, -«|, X\ -y,Ki- K" 1 ) 



(CIO) 



2 ^=£ /' ' e *W< 



2(27r) 3 <(T ; 



2 



1 



VI - xi - y 



U 2 ip3 1<T2X {x 1 , nf,y, fir 1 , 1 - 2i - y, -KTf - fir 1 ). 



The vertices are 
Ux = xl 



-2((4-^).e^) (a^Ka^) (o^Ko^) 

1 1 Xa[ 



2i -y 



y 



and 



U2 = xU 



1 — 2i — y 1 — 2i 

/ is the instantaneous interaction given by 

4e 2 



y 



X — CT2 



dyd k ip 2 (y,K ,l-y,-K ) 



(x - y) 2 ' 



(C.ll) 



(C.12) 



(C.13) 



(C.14) 



2(2vr) 3 

Substituting the expressions for ip 3 in terms of i[) 2 in Eq. ( |C.9| ) we get the bound state 
equation for ip 2 . 
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The bound state equation is 



M 



) 2 + m 2 



x 



x 



(-— 

V F Or - 



2(2tt) 



dy / d 2 q A 



-E (x - y) 
1 1 



V\ifo{x, ^,1-x, -K- 1 ) + V 2 ip2{y, q ± , l-y, -q A 



E> (y 



X) 



«-\ 1 - a;, -^) + 7' 2 ^ 2 (y, 1 - y, -? ± )]) 

4g2 y i 
y dyd 2 q L i) 2 (y, q X , l-y, -q ) . _ , 2 - 



(C.15) 



2(2vr) 3 

Since we are interested in the non-relativistic limit of the bound state equation, we have 
kept only those terms of the vertex which survive in this limit. Also, it has been shown 
that in this limit the helicity dependence of the wave function goes away, so it is sufficient 
to look into one helicity sector. We are looking into the helicity up-up sector and the 
linear in m terms in the vertex, where m is the electron mass, that cause helicity flip are 
not considered. However, both the linear m terms and the terms quadratic in m in the 
bound state equation coming from the vertex are suppressed in the non-relativistic limit. 
The last term in the right hand side of Eq.( |C15| )) is the instantaneous interaction where 
we have neglected annihilation effects in this limit. The energy denominators are 



E = M 



E' = M 2 



(q J 



) 2 + m 2 



) 2 + m 2 



y 

(k 1 ) 2 + m 2 



1 — x 

(q L ) 2 + m 2 



\K - q-^ 
x — y 

(k 1 ~ q L ] 



(C.16) 



x y (y-x) ' 

It is to be noted that the entire mass dependence comes from these energy denominators 
in the non-relativistic limit. 

Substituting these expressions in the bound state equation and considering only the 
photon exchange interactions in the non-relativistic limit, we obtain the non-relativistic 
bound state equation 



[M 2 - 



(fir 1 ) 2 + m 2 



x 1 



x 



}lp 2 (x, K" 1 , 1 — X, -K ± ) 



I dyd 2 q-My,q J 



+ 



711 



y,-q 

2 



- q J 



+ - 

V x 



(C.17) 



We introduce a new three vector 
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(C.18) 



where k z is defined through a coordinate transformation from x G [0, 1] to k z G [—00, 00] 
by 

x=\ + . Kz (C.19) 

In the non-relativistic limit, 

x ~J-(rn + Kg ), (C.20) 
2m 

(k 2 — k ) , , 

rr — y ~ — , (C.21) 

y 2m ' v ; 



x(l — x) 
Then 

2 



(/ " )J " / " 2 =4((p) 2 + m 2 ). (C.22) 



/;/ \ - 1 / m \ 2 



x y ( K ±_ g ±)2 + ( ^ )2(a; _ y)2 ^(i-^y ( K ±_ g ±)2 + ( _^_) 2(;r _ y)i 



We introduce the bound state wave function <p(p) which is normalized as 

J dxd 2 K ± \ ^{x,^) f = J d 3 pJ{p)0*{p)4){p) = 1 (C.24) 

where in the non-relativistic limit we have suppressed the spin which factorizes out. The 
Jacobian of the coordinate change can be written as 

J ® = trL- < a25 > 

The tilde on the wave function is removed by suitable normalization, 

J d 3 pJ{p)0*{p)4)(p) = J d 3 p<f)*{p)<t){p) = 1. (C.26) 
The bound state equation then becomes 

M 2 - 4(($ 2 + m 2 )]0(p) = -^-y 3 J d 3 p^J(mP')mj^-^- 2 - (C27) 

This can be written as 

r i 2p 2 f -t Am 

M 2 - 4((p) 2 + m 2 )]0(p) = J d 3 Mp')j^y 2 - (C28) 

The instantaneous interaction and the photon exchange interactions combine to give the 
Coulomb interaction in the non-relativistic limit. 
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APPENDIX D 



Poincare Generators in Light-Front 

QCD 



In this appendix we derive the manifestly hermitian kinematical Poincare generators 
(except J 3 ) and the Hamiltonian in light-front QCD starting from the gauge invariant 
symmetric energy momentum tensor 0^. To begin with, is interaction dependent. 
In the gauge fixed theory we find that the seven kinematical generators are manifestly 
independent of the interaction. 

We shall work in the gauge A + = and ignore all surface terms. Thus we are working 
in the completely gauge fixed sector of the theory The explicit form of the operator 
J 3 in this case is given in Ref. [|j] which is manifestly free of interaction at the operator 
level. The rotation operators are given in chapter 7. 

At x + = 0, the operators K 3 and E 1 depend only on the density 
ward calculation leads to 

= id + ip + + d + A i d + A i . 



0" 



Then, longitudinal momentum operator, 



A straightfor- 
(D.l) 



- I dx d 2 x J 
2 



id + ip + + d + A j d + A j 



(D.2) 



Generator of longitudinal scaling 

4 



dx d 2 x x x 



dx d x x~ 



id + ip + + d + A j d + A j 



(D.3) 



Transverse boost generators 



E' 



- dx d x x l Q 
2 

1 



id + ip + + d + A>d + A> 
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(D.4) 



The transverse momentum operator 

pi = I J dx~d 2 x ± Q +i 



(D.5) 



which appears to have explicit interaction dependence. Using the constraint equations for 
and A~, we have 



e +i = 



3 + A j d i A j - d + A j d j A i + d+A l &A 3 - 2^ +t Aty H 



Thus 



-5/ 



(D.6) 
(D.7) 



(D.8) 



Thus we indeed verify that all the kinematical operators are explicitly independent of 
interactions. 

Lastly, the Hamiltonian operator can be written in the manifestly Hermitian form as 



P~ = \ J dx~d 2 X ± Q + ~ = ^J dx~d 2 X^{n + Hint) 

where Hq is the free part given by 



no^-AUdfAi+e 

The interaction terms are given by 



id+ 



-(g ± ) 2 + m% t 



id+ 



'Hint — 7~Lqqg + ^ggg + H-qqgg + ^qqqq + ^gggg^ 



(D.9) 



(D.10) 



:d.iu 



where 



^<?<?9 — 



+ S^)^ • + ' ^)(^ • A^ 

- mg^^ ■ + mg{}{^ • A^£ 

+ mg{^){^ ■ A^ - mg^^Ka 1 - ■ A^}, (D.12) 



H 999 = 2gr bc [d l AiAlAi + {d l A i a ) — (A{d + A{) 



(D.13) 
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nj 2 
ri qqgg {J 



+ 4^(r bc Ai ) d + Al)^T a 0], 



(D.14) 



(D.15) 



2 



+ 2^(^-4)^(49+^)]. 



(D.16) 



133 



APPENDIX E 



Transverse Spin in Free Fermion 

Field Theory 



E.l Poincare Generators: Operator Forms 



The symmetric energy momentum tensor 



The momentum operators are given by 

1 r 



I /■ I 

P + = - dx~d 2 x ± ^ + - id + ip 



St +A id + }^ 



The Hamiltonian operator is 

1 r 1 1 

P~ = -J dx~d 2 x ± ilj^y~- id + +7 + ^ id~ i> 

= \ J dx-£x^-^[ml - (d^) 2 ] - (^K " (^W) 
The longitudinal scaling operator (at x + = 0) is 



— / dx d 2 x J ~x~ 
2 



— 1 ~ " 



(E.l) 



(E.2) 



(E.3) 



(E.4) 



(E.5) 
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The transverse boost operators are 

E l = ~j dx~d 2 x ± x i 



— 1 ~ 

■^7 + - id + ip 



The generators of rotations are 

J 3 = l - I dx-d^lx 1 



and 



r • ' i /7/ J + 7 2i 



7 ^ + }^ 



— x 



= J dx-d'x^i'-ix 1 ^ 2 - x 2 ^\ - ['-(x 1 ^ 2 - x 2 T 



-II 



dx d x <x~ 



- r 



1 1 " 



= l -Jdx~d 2 x ± £} 



x \m 2 -{d^) 2 )^-x-^- 



_l J ie 3 <j 

9+ l dx i dxi 



- I dx d 2 x 
2 J 

d_ 

d + ^ dx 



d+ 
+ e lj ma j 
1 



x\m 2 - (c^) 2 )-- - x- ^~ 



d+ 



dx i 



1 r 







OXl > 



E.2 Fock Representation 



Free spin-half field operator is 

dk + d 2 k 



a 7 2(2tt)VA; + 
In terms of Fock space operators 



P + 



r dh + d 2 h L 



135 



Hk + H 2 k ± 

2(2^+** E [ fet ( fc ' A )^> A ) + ^> " A )^> "*)]■ 



(E.ll) 



p ~ = / 2 XT 3 ^ m ^ +(fc±)2 ? [6t( *' A)6( ^ A) + rft(A; ' ~ AWA; ' ~ A)] - (K12) 



« /■ dk + d 2 k^ + ^ / 
"2 7 2(27r¥k+ h 



2(2vr) 3 fc+ 



*(M)^ A) + ^-A) d(l% _ A)] 



<9A;+ 



dk+ 



_ [6 » ( ,, A) ^M) +(i t (4> _ A) ^|^) ]) 



<9A;+ 



(E.13) 



[& t (A:,A)^|^ + rft(A; ! -A) M(A; '- A) 



(E.14) 



J 



[kl w* ~ k2 w m A) + ( [kl ^ ~ e w ]d1 ^ - x)] ) d{k > " A) 

-d\k,-X)[k^-k^]d(k,-\)] 

1 r dk+d 2 ^ 

+ =: I Q ? E A t fet (^ A)6(/fc, A) + d f (A:, X)d(k, A)] (E.15) 



2 7 2(2vr) 3 A; 



with A = ±1. 

F l 



dk+d 2 k L 
2(2vr) 3 fc+ 



A 

[&'(*, A) 



^(*,A) 6(4A) + c^-A) d(|fci _ A) 



g6(fc. A) 
<9A;+ 



+ </*(*;, -A) 



dk+ 
dd(k, -A) 



+ 



« /• d/rc 
2 J 2(27 



2(2tt) 3 A;+ £;+ 

db{k, A) 



a&t(fc, a) 



dk+ 
b(k,\) + 



) 



dd^(k,-X) 



dk l 



d(k,-X) 



[b*(k,\)- 



dfr 



+ S(k,-X) 



dd(k, -A) 



f^T j^' E ^av A)6(fc, A') - d\k, -X')d(k, -A) 



dk+dtk 1 - eg 
2(27r) 3 A;+fc+ 



m F ]T ct aa' [ fet (^> A)6(fc, A') + dt(Jfe, -X)d(k, -A')] . (E.16) 

AA' 
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E.3 Transverse Spin of a Single Fermion 

For a single fermion of mass m and momenta (fc + , A;- 1 ), we have 



P+ 
P~ 

K 3 
F 1 



fcA) = k + | kX), P 1 | kX) = k 1 | kX), P 2 \ kX) = k 2 \ kX), 
kX) = {k±) l + ™ 2 | fcA), J 3 | ifeA) = ^A | fcA), 

r\ r\ r\ 

kX) = -ik+— | kX), E 1 | fcA) = ik + — | fcA), E 2 \ kX) = ik + — \ kX), 



kX) 



dk+ 1 " 1 ' dk 1 1 " 1 ' dk 2 



We arrive at 



mj 1 | fcA) = {^F 2 P + + K 3 P 2 - 1 -E 2 P- - P 1 ^ 3 ) | kX) 

= m^^flkX'), (E.18) 

A' Z 

mj 2 | fcA) = [~F X P + - K 3 P l + ^^P- - P 2 J 3 ^j | fcA) 



mE^I^')" (E.19) 

A' Z 
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APPENDIX F 



Transverse Spin in Free Massless 
Spin One Field Theory 

F.l Poincare Generators: Operator Forms 

The symmetric gauge invariant energy momentum tensor 

Ql*» = p>* F V _ (F.l) 

where the Lagrangian density 



C = -\f^F, u (F.2) 



with 

= d v A» - d»A v . (F.3) 

We choose A + = gauge. Only the transverse fields A 1 are dynamical variables. The 
momentum operators are given by 

p+ = I j dx-d 2 x ± d + A j d + A j , (F.4) 

pi = I J dx~d 2 x ± (A j d + d j A i - A l d + d j A 3 - A j d + d l A 3 ) . (F.5) 

The Hamiltonian operator is 

P" = 1 y dx~d 2 x ± \^{d + A-f + ^F^Fij 

= \j dx~d 2 x ± d i A j d i A j = ~J dx~d 2 x L A\d i ) 2 AK (F.6) 

The longitudinal scale generator (at x + = 0) is 

1 



K 3 - — 
2 



J dx~d 2 x ± x'd + A 3 d + A 3 . (F.7) 
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The transverse boost generators are 

E i = ~ J dx~d 2 x ± x i d + A j d + A :i . (F.8) 
The generators of rotations are 

J 3 = ^Jdx-d 2 x ± (x 1 [d + A 2 d i A i + d + A 1 (d 2 A 1 -d 1 A 2 )} 
-x 2 \d + A l d l A l + d + A 2 (-d 2 A 1 + d l A 2 )}) 
= l -J dx-d 2 x L (x 1 [d + A 1 d 2 A 1 + d + A 2 d 2 A 2 ] - x^d+A'd'A 1 + d + A 2 d 1 A 2 }) 

+ X - J dx~d 2 x ± [A 1 d + A 2 - A 2 d + A 1 ]. (F.9) 

and 

F l = X - J dx-d 2 x ± (x~ {A^d+WA 1 - A i d + d j A 1 - A j d + d i A j ) 

-x^A^&YA^ - 2 J dx-d 2 x ± A i r t ij d j A j , no summation over i, j ,(F.10) 

with i] 12 = i] 21 = 1, f] 11 = rf 2 = 0. 
F.2 Fock Representation 



The dynamical components of the free massless spin field operator in A + = gauge 



are 



2 r Jh+^uL 



In terms of Fock space operators, we have 

dk+d 2 ^ 



(F.ll) 



(F.12) 



J 



(F.13) 



(F.14) 



K 



i r dk + d 2 k^ , „ r 
~2J 2(2n) 3 k+ 4^ 



2 J 2{2nfk 



da]{kA) .)a(k,X)-a\kA) da{kA) 



dk+ 



dk+ 



(F.15) 
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d rf*V E [( aat jf.' A) )a(fe, A) - a\k, \) da{ *\ X) )}. (F.16) 
2 7 2(27r) 3 fc+ Y ^fc* 



j3 - \s^^ k w- k w>^y^ 



2(27r) 3 fc+ 

Introduce creation and annihilation operators 



+ l / ^7^ZT3T+ ( flt (^' 2)a(fc, 1) - a^fc, l)a(fc, 2)). (F.17) 



a(fc, T) = ^[a(fc, 1) " ™(*, 2)], a(fc, |) = -^=[a(fc, 1) + ia(fc, 2)]. (F.18) 



Then 



J3 = U^^w-^w^)^ 

-aHk,\)(k^-k^)a(k,\)} 



where A now denotes circular polarization. 



i /■ dk + d 2 k 1 - (fc -1 ) 2 ^ .daUk,X) ^ + ^.da(k,X). 

+2 7 2(27FF^ ? (-at;-* A ) - at ^> ) dfc^ 



/■ rlk+^k 1 - k'j 

F.3 Transverse Spin 

Using the explicit form of the operators, we get for a state of momentum k(k + ,k ± ) 
and helicity A, 

JT 3 | fcA) = | fcA) = A | fcA), 
| kX) = k l X | fcA), 
W 2 | fcA) = fc 2 A | kX), 

W~ | fcA) = ^^A | fcA). (F.21) 
fc + 

T | fcA) = 0. (F.22) 
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APPENDIX G 



Transverse Spin for a System of Two 
Non-interacting Fermions 



In order to show the non-triviality of the transverse spin operators even in the free 
theory and the manifest cancellation of the center of mass motion in this case, here we 
evaluate the transverse spin for a composite system of two free fermions. The mani- 
fest cancellation of the center of mass motion for the interacting theory is much more 
complicated and is given in appendix H. 

Let the mass of each fermion be m and momenta {kf , k^~), i — 1, 2. We take the state 
to be | P) = Si)tf(k 2: s 2 ) | 0), where Si and s 2 are the helicities. 

MJ 1 | P) = (^F 2 P + + K 3 P 2 -^E 2 P~ - P 1 ^ \ P). (G.l) 

We introduce Jacobi momenta, (x^g -1 ) defined as 

k± = q- 1 + x 1 P ± , k£ = -q 1 - + x 2 P ± kf = Xi P + (G.2) 

with J2^i = 1. 

Here M is the mass of the composite system and (P + , P^) are the momenta of the 
center of mass. 

The partial derivatives with respect to the particle momenta can be expressed in terms 
of these variables as 



and 



d 


d 


d 




d 




d 


d 


dk\ 


X<2 dq i 


+ dP v 




dk\ 


-X! 


dq i 






d 


x 2 d 


d 


P ± 


d 








dkf 


P+dx ± 


+ dP+ 


~ X2 pT 


dq- 


_ ' 






d 


x\ d 


d 




d 








dkt 


P+dx 2 


+ dP+ 




dq A 







(G.3) 



(G.4) 



(G.5) 
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Then we have 



(G.6) 



- \e 2 P~ \P)=[- \{{P^f + M 2 )^ + *P 2 ] | P> , 



(G.7) 



PV 3 1 p) = V|r - s 1 ^) + P l8 i + plS i\ I ^> 



(G.8) 



-F 2 P A 
2 



Substituting 



we get 



d 



d 



si q 



' dxi 

2 A *1 



• s - ? ! '/ ; ' 2 + (g ± ) 2 )(^ 2 - 



2x 2 + 2 (m ' 



d 

%2 dq 2 



a 



+ iP + P' 



^) - *g 2 P 

d i 
dP+ ++ 2 

d P 1 



dq- 

+ +km 2 + (q ± ) 2 ) 



9 .„ 2 / 9 



+ 



A, 

1 9 9 
h -(P ) 

dP 2 2 1 ; <9P 2 



M 2 



(m 2 + (g^) 2 ) 



MJ l \P) = tq 2 (x 2 - x x 



d 



d_ 



' 2 + (? ± ) 2 )(- 

Xi To na z 



+ -in i 



xi. d 
x 2 dq 2 



+ 



2 l a;i x 2 j 2 &i ^2 



Explicitly we see that MJ 1 does not depend on the center of mass momenta. 



(G.9) 



(G.10) 



(G.11) 
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APPENDIX H 



Transverse Spin of a Dressed Quark 
in an Arbitrary Reference Frame 

We introduce a wave packet state 

\^) = lfdP + d 2 P ± f(P)\P,a) (H.l) 

which is normalized as 

(Vv|Vv) = <W- (H.2) 

Here f(P) is a function of P, the exact form of which is not important. Using the 
normalization condition of the state we get 

X - J dP + d 2 P ± f*(P)f(P)(2n) 3 P + = 1. (H.3) 

The expectation values of the various operators involved in the definition of MJ % are given 
below. It is to be noted that we have done the calculation in an arbitrary reference frame, 
in order to show that the dependence on the total center of mass momenta (P + ,P ± ) 
actually gets canceled in the expectation value of MJ % . 

The matrix elements presented below have been evaluated between wave packet states 
of different helicities, namely a and a'. Since the transversely polarized state can be 
expressed in terms of the longitudinally polarized (helicity) states (see chapter 8), the 
matrix elements of these operators between transversely polarized states can be easily 
obtained from these expressions. We introduce 

Vi = f(PW, C x a = f{P)K^ (H.4) 

The matrix elements are given by 

<^|K 3 P 2 |W> = ljdPU 2 P\2nfP + (^r^P + P 2 

+ Wj dxd^P^ x d -§^P + + h.c). (H.5) 

spin 

143 



Wo 



-E 2 P 



free 



4E/^^ <A ' rm2+(9±+xP±)2 

spin 



OP 2 



X 



+ 



{-q ± + (l-x)P ± f 



1-x 



+ h.c). (H.6) 



. g - / dP + d 2 P ± (2n) 3 P + V / cfcdV-^= 



(.,a* d ^\ t 



1 — x 



X 



1. (l-g) ,. ± J 

im cr -e 

2 x v 



+ h.c.y 



(H.7) 



Here /i.e. is the Hermitian conjugate, X) sp j n is summation over a±, a[, A, A'. -Pj ree is the 
free part and P^ nt is the interaction part of the light-front QCD Hamiltonian density. 

| P X J Z | i>o') = - J dP + d 2 P ± (2n) 3 P + [j2 J dxd 2 q ± 

~ ' spin 



+ 



A,(T2,cr; 



(H.8) 



The first term in the above expression is the quark-gluon orbital part, the second and the 
third terms are the intrinsic helicities of the quark and gluon respectively. Finally, the 
operator ^F 2 P + can be separated into three parts, 



-F 2 P + = -FfP + + -Ff r P + + -Fl T P + 
2 2 1 2 11 2 111 



(H.9) 



where F 2 ,F 2 j and F] H have been defined earlier. The matrix elements of the different 
parts of these operators for a dressed quark state in an arbitrary reference frame are given 
below. A part of this calculation has been given in some detail in appendix I. We have 



\f 2 p+ = \f 2 {i)p + - \fKi)p + - \fK^)p + - 

The matrix elements of these three parts are 

<W | \Fl{l)P + I VV> = bl>°\\j dxd 2 q^x~P 2l -P + | VV> 
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(H.10) 



dP + d 2 P ± (2n) 3 P + 



p+p2 



2^ dP+ 



+ 



dib a ' 



spin 

i 

~ 2 



E / dx^^^-p 

""in 



spin 

i 
2 



j j2 _L 2 / *ct ^cr'i-V 



+ /i.e. 



(H.11) 



k | ip 2 (2)P+ | W> = | ± y dxd 2 q^x 2 P l -P + | ^! 

= ^^w ( 2.)3p + ^r^^^p + 

spin 



-) 2 | (-g ± + (l-x)P i 



^2, 



X 



1 — rr 



spin " 

spin ^ 



spin " 
+ 



1 -x y 



i E / + H • (H - 12) 

spin " 

In the above two equations, both the single particle and two particle diagonal matrix 
elements contribute. 

ip|(3)P+ | W) = (VV I i / dxd^x 2 P int 1 -P+ | W> 



^2(2tt) 3 ^ 



i| dP + d 2 P ± (2ir) 3 P + j dxd 2 q L -j= 

I \ - / * I I \ / ~ I l\ / 



dP 2 



(g ± -e ± ) 1 (5- x - e x )(o- x • g- 
%m 

2 x 



) 1. (l-x) ,_ ± ±1 
2 x J 

:' lX + h.c). (H.13) 

Only the off-diagonal matrix elements contribute in the above equation. The matrix 
elements of the three different parts of \F] I P + are given by 



1 — x 

T: Jb 



^JdP + d 2 P ± (2n) 3 P+ 



in 



145 



m 

2 ^ 

spin 



, 1 



dxd'q^xxWXa^y-l , (H.14) 



(H.15) 



spin 



+ -jdxd 2 q ± £ APty&V& a 



| -F 2 gII P + | W> 



4 -/2(2tt) 3 



£ \ j dP + <fP ± (2<ir) 3 P + J dxd 2 q A 



\fl=x~ 



spin 



a; 



+ /i.e. 



H.16) 



In Eqs. (H.14) and ( [H.15 ), contributions come from only diagonal matrix elements 
whereas Eq. (|H.16|) contain only off-diagonal matrix elements. The matrix element of 



1 Ff II P + is given by 



<W | -F 2 H P + | W) 



-.- J dP + d 2 P ± (2ir) 3 P + [j2 J dxd 2 q ± y=L= 

^ spin V ^ 



2(2tt) 3 ^^ vi - s 



A,(Tl,CT 1 

A,e7-i,cr'i 



Finally, the expectation value of the transverse spin operator is given by 
(W | MJ 1 | = <VV I ^ 2 ^ + + K 3 P 2 - X -E 2 P- - P X J 3 | Vv>- 



(H.17) 



(H.18) 



From the above expressions it is clear that all the explicit P x dependent terms get canceled 
in the final expression. To be specific, it can be easily seen that all the terms in the 
expectation value of K 3 P 2 — \E 2 Pj ree — P 1 <7orbitai are P 1 ' dependent and they exactly 
cancel the P x dependent terms in \Fj^ ree ^P + \ the two P 1 dependent terms in the intrinsic 
part of P l J 3 exactly cancel the two similar terms in the expectation value of ]^F 2 I P + + 
^FfjjP + and the expectation value of \E 2 P^ nt completely cancel all the P 1 - dependent 
terms in the expectation value of \F 2 ^ int ^P + . 
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APPENDIX I 



Details of the Calculation 



Here, we explicitly show the evaluation of one of the matrix elements of the interaction 
part of Ff. Consider the operator, 

P+ 



O g = - I dx d 2 x ± x 2 g 



This can be written in Fock space as 



(a ■ d + m) 



d+ 



(1.1) 



° 9 = § E l(dk 1 )J(dk 2 )J[dk 3 \(b\k 1 ,s 1 )a(h,X)b(k 2 ,s 2 ) 



Sl,S2,A 



f {dX ' k l t ™V ■ e^xs^^mh -k 2 - h) 

+ tf(ki, Sl )a\k 3 , X)b(k 2 , s 2 ) X \ } aL ' im) ■ ef) X . 



k\ 



S'2 



(1.2) 



where (dk) 



dk+cPk 1 



and [dk] 



dk+cPk 1 - 



. We evaluate the expectation value of this 



2(27r)3vT L """ J 2(27r)3fc+' 

operator for the dressed quark state given by Eq.( |H.l| ). Only the off-diagonal parts of the 
matrix element will give non-zero contribution. The matrix element is given by 



(Wiojvv) = f £ J(dpyJ{d Pl } J{d V2 }p{2^fp + 5\p- Pl -p 2 ) 



r r~Z + (v • P — im) 

ttiyptxl- — p+ — - 



<9P 2 



1 



2(2tt) 3 <5 3 (P - Pi - p 2 ) -P+ + h.c.} 



(1.3) 



<Tl,A 
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ig 1 



^ J(dPy J dxd 2 q^y= 

(71, A v 



— Xa(cr P -im)(a -e A )x 



<9P 2 



+ /i.e. 



(71, A 



^ / (dP)' / rfg 



(TlA 



X^V^ • e A)Xa' + h.c. 



;i-4) 



where {dp} 



dp+d 2 p A 



v /i(2^V 
evaluated in a similar method. 



and (dP)' = \dP + d 2 p L 2{2'KfP + . The other terms can also be 
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